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Abstract 

We study different aspects of integrable boundary quantum field theories, focussing 
mostly on tlie "boundary sine-Gordon model" and its applications to condensed matter 
physics. 

The first part of the review deals with formal problems. We analyze the classical limit and 
perform semi-classical quantization. We show that the non-relativistic limit corresponds to 
the Calogero-Moser model with a boundary potential. We construct a lattice regularization 
of the problem via the XXZ chain. We classify boundary bound states. We generalize the 
Destri de Vega method to compute the ground state energy of the theory on a finite interval. 

The second part deals with some applications to condensed matter physics. We show how 
to compute analytically time and space dependent correlations in one- dimensional quantum 
integrable systems with an impurity. Our approach is based on a description of these systems 
in terms of massless scattering of quasiparticles. Correlators follow then from matrix elements 
of local operators between multiparticle states - the massless form-factors. Although, in 
general an infinite sum of these form-factors has to be considered, we find that for the current, 
spin and energy operators only a few (two or three) are necessary to obtain an accuracy 
of more than 1%. Our results hold for arbitrary impurity strength, in contrary to the 
perturbative expansions in the coupling constants. As an example, we compute the frequency 
dependent condunctance, at zero temperature, in a Luttinger liquid with an impurity, and 
also discuss the succeptibility in the Kondo model and the time-dependent properties of the 
two-state problem with dissipation. 
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Chapter 
Introduction 



In one space dimension, a quantum field theory can be defined either on a circle, or on 
an open interval with certain boundary conditions. In Hamiltonian formulation, boundary 
conditions amount to the presence of additional interaction term on the boundary 

H = Hhuik + V{{<j)B}), 

where {4>b} denotes the set of boundary degrees of freedom - fluctuating flelds at the bound- 
ary. The reasons for studying the theories with boundary conditions seem natural since in 
practice one has to deal often with a bounded system having some interface with the external 
world, as in the problem of polymer adsorption p|. Also, the problems on the semi-infinite 
interval in ID sometimes appear as the reductions of 3D problems to the s-wave dependence 
as e.g. in the case of the monopole-catalysed baryon decay 0. 

One of the most appreciated and rich at the present applications of the boundary field 
theory is to the impurity problems of condensed matter physics [|, ^ |70|, |76|, provided 
that the scattering on impurity can be mapped onto the scattering off the boundary with 
some boundary potential. An incomplete list of examples includes the Kondo model llO] , the 



dissipative quantum mechanics |^ and the quantum Hall liquids with constriction [^. All 
of the mentioned above three models look alike from the point of view of the bulk part, which 
is a free massless boson, but differ by the boundary interaction. The role of integrability here 
is two-folded. First, it allows to find the static properties (e.g. via the Bethe ansatz) of the 
model from the bare Hamiltonian, including the mass spectrum, scattering matrices, or the 
free energy. Second, it allows to take advantage of working with the physical excitations in 
the formalism of the FS theory. Namely, integrability suggests a convenient basis of massless 
particle states which are particular combinations of plane waves that scatter diagonally off 
the boundary. Working with these massless particles at first sight adds some complexity, 
but it is paid off by the final simple and manageable results. Let us stress that we deal only 
with such boundary interactions that preserve the integrability of the bulk part. Fortunately, 

^E.g., for the Ising model such an interface can be modeled by a boundary magnetic field 
^ Corresponding classical solutions are presented in 
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there are quite a few of such models in the real world. 

Different integrable boundary theories have been introduced and solved throughout the 
history of integrable models. A noticable contribution in this field has been made by Cardy P, 
^ , who studied the critical surface behavior withing the framework of conformal field theory. 
Another substantial advancement has been done by the work of Ghoshal and Zamolodchikov 
The authors have succeeded in formulating and (partially) solving the set of general 
constraints for the factorizable boundary field theory, thus generating a powerful method 
of obtaining the integrable boundary models from the integrable bulk field theories. The 
equations of [|I| resolve the boundary integrability (boundary Yang-Baxter equation) together 
with the physical constarints of unitarity and crossing symmetry and possess the restrictive 
power to determine the scattering matrices up to "CDD ambiguity". In a more simple 
words, given a FS theory in the bulk with its two-particle scattering matrices, we can derive 
all possible integrable boundary models compatible with this bulk theory. Of course, what 
we get is merely the FS description, leaving the identification of the Hamiltonian structure 
to our intuition. But other methods, developed in [§] and [^, allow us to obtain the exact 
free energy and the correlation functions correspondingly. Thus, the FS approach turns out 
to be very productive. 

The paper [|I| in such a way sets basis for a more systematic study of the boundary 
integrable models, marking them out as a subfield of integrable models. We undertake in 
this dissertation a close study of some particularly interesting boundary models, such as the 
boundary sine-Gordon model (chapters 2,3,4,6), as well as develop some general techniques 
(chapters 5,7). We assume that all the properties of the bulk models of interest are known 
to us, so we can focuse on the peculiar boundary phenomena. One example of the boundary 
phenomena are the boundary bound states (chapter 4). The existence of such states, localized 
at the boundary of the crystal lattice, was first pointed out by I.E.Tamm. Finally, in chapter 
7 we describe in some detail the physics behind the impurity and dissipative two-state models 
and set up the technique for calculating the physical observables (correlation functions) in 
these models by making heavy use of the boundary factorized scattering. Quite remarkable, 
such characteristics as the conductance can be expressed in terms of the boundary scattering 
matrices! 

In chapter 1 we review the Bethe ansatz technology. We start with the traditional Hamil- 
tonian approach and show how to extract physical observables (mass spectrum, scattering 
matrices etc) from a bare Hamiltonian. Although we do not present any new contributions 
with respect to the existing extensive literature on the subject, this material is necessary to 
understand further chapters. We focuse our discussion on the Thirring model, which is the 
fermionic analog of the sine-Gordon model. 

In chapter 2 we consider the sine-Gordon model on a half-line, with an additional poten- 
tial term of the form —M cos f (</? — v^o) at the boundary. We construct the classical solutions 
by using the bulk sine-Gordon theory and the "generalized method of images." From the 

* The term "integrable models" in our context amounts to the model having an infinite number of 
conserved charges in involution (i.e. mutually commuting), both on quantum and classical levels. The 
integrability results to the factorized scattering property by making use of the argument that S-matrix must 
commute with the infinite number of charges 0] . 
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classical solutions in hand we extract the time delay (|2.18|). From the time delay we recon- 



struct the semi-classcial phase-shift using the method of Jackiw and Woo [Q. We establish 
the agreement with the semi-classical limit /3 — ^ of the exact boundary reflection matrix, 
( p.34| ). The exact expressions for the boundary reflection matrices are known up to CDD 
ambiguities . They were obtained as a "minimal" solution to the general set of constraints 
for the integrable boundary field theory. 

The purpose of chapter 3 is to investigate the non-relativistic, 6 —>■ 0, limit of the bound- 
ary sine-Gordon model and to determine the quantum-mechanical potential induced by the 
presence of boundary. We show that the generalized Calogero-Moser model with boundary 
potential of the Poschl- Teller type describes the non-relativistic limit in question. 

In chapter 4 we address the exact quantum field theory solution of the boundary sine- 
Gordon model, which we obtain by means of the Bethe ansatz technique. Among other 
things, this solution allows to re-derive the boundary reflection matrices of ( |3.1| )-( p^ , 
and to relate them to the physical parameters in the Hamiltonian The present chapter 
includes a complete study of boundary bound states and related boundary S-matrices for the 
sine-Gordon model with Dirichlet boundary condition. Our analysis is based on the solution 
of the boundary bootstrap equations, representing the integrability constraints, together 
with the explicit Bethe ansatz solution of the inhomogeneous XXZ model in a boundary 
magnetic field - a lattice regularization of the boundary sine-Gordon model. We identify 
boundary bound states with new boundary strings in the Bethe ansatz. 

The main purpose of chapter 5 is to study the ground state energy of 1+1 integrable 
relativistic quantum field theories with boundaries. This involves several questions. One 
is the energy associated with a boundary for an infinite system, the other is the way the 
energy of the theory on an interval varies with its length - the "genuine" Casimir effect. 
The elegant method of Destri and de Vega for the periodic systems leads directly to the 
expression for the ground state energy from which the infinite size contribution and the finite 
size correction can be easily extracted. The heart of the DDV method is a non-linear integral 
equation ( |5.45D being derived from the Bethe equations. We generalize the Destri-de-Vega 
method to the systems with boundaries and apply it to compute the ground state energy for 
the boundary sine-Gordon model. 

In chapter 6 We study an open XXZ chain in the regime A > 1 with a boundary 
magnetic field h and discuss some of its peculiar features due to the presence of boundary. 
In the Bethe ansatz formalism, boundary bound states are represented by the "boundary 
strings" as described in chapter 4. We find that for certain values of h the ground state wave 
function contains boundary strings, and from this infer the existence of two "critical" fields 
in agreement with An expression for the vacuum surface energy in the thermodynamic 



limit is derived and found to be an analytic function of h. We argue that boundary excitations 
appear only in pairs with "bulk" excitations or with boundary excitations at the other end 
of the chain. The case where the magnetic fields at the left and the right boundaries are 
antiparallel has non-trivial differences with the case of the parallel fields. The Ising (A = oo) 
and isotropic (A = 1) limits are discussed thoroughly and found helpful for the intuitive 
understanding of the behavior of the boundary XXZ chain at arbitrary A. 

In chapter 7 we show how to compute analytically time and space dependent correlations 
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in one-dimensional quantum integrable systems with an impurity. Our approach is based 
on a description of these systems in terms of massless scattering of quasiparticles |^ . Cor- 
relators follow then from matrix elements of local operators between multiparticle states - 
the massless form-factors. Although, in general an infinite sum of these form-factors has to 
be considered, we find that for the current, spin and energy operators only a few (two or 
three) are necessary to obtain an accuracy of more than 1%. Our results hold for arbitrary 
impurity strength, in contrary to the perturbative expansions in the coupling constants. 
As an example, we compute the frequency dependent condunctance, at zero temperature, in 
a Luttinger liquid with an impurity, and also discuss the succeptibility in the Kondo model 
and the time-dependent properties of the two-state problem with dissipation. 

This review is based mostly on the published papers [|^, [^, [^, ||70| . 
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Chapter 1 

Introduction to the Bethe Ansatz 



One of the efficient approaches to quantization of interacting fields is based on the conformal 
field theory (CFT), while yet another on the factorized scattering theory (FST). In both 
cases the knowledge of the Hamiltonian is not required. Rather, the data is encoded into the 
particular representation of the Virasoro algebra on the space of fields and the dimension of 
the perturbing operator, or in the exact scattering matrices as in the case of FST. We start, 
however, with the traditional Hamiltonian approach and show in this chapter how to extract 
physical observables (mass spectrum, scattering matrices etc) from a bare Hamiltonian us- 
ing the Bethe ansatz technology. Although we do not present any new contributions with 
respect to the existing extensive literature on the subject |T0[, this material is necessary 
to understand further chapters. 

This chapter is by no means the complete review of the Bethe ansatz. We focuse our 
discussion on the Thirring model, which is the fermionic analog of the sine-Gordon model. 
The relation between the two models becomes an exact mapping of one onto another due to 
the bosonization technique of Coleman and Mandelstam. 



1.1 Hamiltonian formulation 

The massive Thirring model is defined by the Hamiltonian 

Ht = J dx[-i{^tdx^i - ^td^^2) + mo{^t^2 + ^t^i) + 2g^t^+^2H- (1-1) 
It can be rewritten in terms of the creation and annihilation operators of the Fock fermionic 



space and diagonalized [|rT|. The vacuum of the Fock space |0) is a particular eigenstate 
annihilated by ipi and il'2 and is called bare vacuum. The physical vacuum is the state with 
the Dirac sea ffiled which has an infinite negative energy. Thus, to perform calculations a 
high-energy cutoff is required. Physical excitations are obtained by removing pseudoparticles 
from the Dirac sea and placing them above it, allowing in general some complex combinations 
called bound states. 

The model possesses the conserved charge N = J dx{ijjfi/ji + il^2'^2) (total number of 
pseudoparticles), as well as an infinite family of other local commuting charges, and there- 



fore it is integrable. The latter can be shown by passing to the discrete Thirring model |T2|, 
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which is known as lattice XYZ chain.0 The presence of conserved charges puts tremendous 
constraints on the dynamics of the model and, in fact, makes it analagous to the free Dirac 
fermions. The non-trivial difference comes from the mutual pairwise interaction of pseu- 
doparticles in the vacuum. Loosely speaking, the Dirac sea is sensitive to the removal of one 
of its pseudoparticles, i.e. a polarization of vacuum occurs. In yet another sense the properly 
regularized Thirring model is analogous to the quantum-mechanical N-body problem with 
the pairwise interaction potential V ~ 6{x — y). 

N-particle bare wave-functions are constructed by glueing up free solutions to the Dirac 
equation at the boundaries of the domains Xi^ < < . . . < Xij^ , similar to the way one 
proceeds in the quantum-mechanical problem with the delta-function potential. The solution 
to the free Dirac equation in two dimensions can be written in the form 

^ = M(/3)e^P^-*^*, (1.2) 

where 

« = ^ ( 1^/2 ) , E^-p' = ml 
it is convenient to parametrize the energy and momentum in terms of rapidity j3: 

E = mo cosh (3 , p = mosinh/5. (1.3) 



Matching solutions ( |1.2|) leads to the pairwise phase shifts in the wave-function 

, , sinhi(2m-/5) 1 , , 

arising at the hyperplanes Xi = xj, with the bare S-matrix being 

S = e"^, 5(0) = 1. (1.5) 

The domain fi < 7r/2 is called a repulsive regime, while the domain /i > 7r/2 is called an 
attractive regime. The value /i = 7r/2 is the free point. 



1.2 Bethe equations, thermodynamic limit and 
solution for the density of states 

The above heuristic discussion can be made more rigorous. Let us consider first the model 
( |1.1| ) compactified on the circle of length L and then take the limit L ^ oo. The wave- 
functions of the state^ 

N 

\(3i,f]2,...,f3N)= |[l^''^^'^"'d^^ n [l-^X{P^-f3J)e{x,-x,)]■ (1.6) 

i=l i<j<N 

■A+{P,,x,)---A+{Pn,xn)\0), 

^There exist many lattice models that have Thirring model as the continuum limit. XYZ chain is, however, 
the most studied one. See Q for another example. 

^ Notice the peculiar to integrable models factorized structure of the N-particle states: the wave-function 
consists of the product of two-particle terms. 
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where 6{x) is the step function and 



are subject to the periodic boundary conditions, called Bethe equations: 

Taking the logarithm of Eq. ( |1.7|) we get 

- p{f3i)L = ^ <I>(A - + 27ili. (1.8) 

j 

Such equations should be written for each particle's index i and form together a complex 
set of coupled transcendental equations. Different solutions are obtained for different sets of 
integers /j. For example, the Dirac sea corresponds to the dense set of integers from — / to 
J, while excited states correspond to the sets with u of /j's missing. Such vacancies will be 
called holes. Thermodymanic limit is the limit L oo, N oo. Class of states for which 



u is fixed and finite when iV — oo is called scaling states |]T4 . 

Besides the real solutions to (|1.8|) and the Dirac sea solutions with Im/3 = vr, there exist 
various other solutions with some of complex. In the thermodynamic limit it is possible 



to show that all permissible complex roots form strings. From the point of view of 
the wave function, strings are bound states. They contain the rapidities with common real 
part Re/3 and equally spaced imaginary parts and are located symmetrically with respect 
to Im/5 = or Im/5 = vr lines, (see Figure 1.1). The easiest way to determine the spacing 
along imaginary axis is to look at the poles of the bare S-matrix ( |1.5| ). The latter are given 
by Ap = 27iil — 2i^. Since by periodicity —it <Im/3 < vr, it is enough to consider / = 0, 1. 
We find that A(3 = — 2i/i in the repulsive regime, and Af3 = —2iuj in the attractive regime, 
where 

u; = vr — yU. 

More rigorous classification of strings based on the analysis of Bethe equations is given in 
Denote by 

m—l n—1 

*m,n= E E*(/5 + ^/^('^-2p-^ + 2/)) (1.9) 

p=0 «=0 

the scattering phase of m-string on n-string, and by 

= Y.^{(3 + m + tf^{n - 1 - 2p)) = <l>i,„(/5 + m) (1.10) 

p=0 

the scattering phase of the Dirac sea pseudoparticle on the n-string. For technical simplicity 
we restrict /z and uo to certain "rational" values equal to vr/t, where t = 2,3,.... Then 
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lm(3 





vr Dirac sea 




• 




• 








(a) 




*(b) 








*(c) 





Figure 1.1: Various root configurations: (a) 1-string; (b) 2-string; (c) 3-string. 

only the strings of tlie lengtli 1,2, . . . ,t — 1 exist. | Note tliat (t — l)-string is qualitatively 
different from all the shorter strings and is similar to a hole in the Dirac sea as far as the 
bare energy, momentum and scattering phase with other objects are concerned. 

In the thermodynamic limit the distribution of roots of Bethe equations is described by 
continuous positive functions Pj{P) and Pj{P), where Pj{P) is a density of j-strings, p_i is a 
density of the Dirac sea pseudoparticles, and pj is a density of holes. Summing the Bethe 
equations and rearranging, we get the equations for strings: 

Ps{f3,)L = - E E - P,) + 2vr/,„ (1.11) 

I j 

where 

s-l 

= E MA + - 1 - 2a)), 

a=0 

p_i(/5) =p(/5 + ^7r), 

and Pi are now the real numbers. Introduce =sign(/j+i — Ii)s- Since the density of roots 
must be positive, then 

%^-^.(P + P).. (1-12) 



^At first glance it seems that the t-string should be allowed, too. However, it is easy to check that it has 
vanishing energy, momentum and scattering phase. So, it is a "ghost." 



12 



Eq. ( |1.12D is the definition of tfie density functions. After tlie tliermodynamic limit is taken, 
the Bethe equations assume the form: 



Introducing 



we get 



I 



/vr 



P 

Ztx , 



;i.i3) 

;i.i4) 
;i.i5) 



where the involution operator * is defined as: 

A * p{x) = J A{x - y)p{y)dy. 
The choice of sign of as is determined by the behavior of the function 

namely, whether ys{(3) increases or decreases. Thus, we are interested in the sign of y's{[3): 

Alas, we cannot compute the infinite sum of the terms whose values depend on unknown yet 
(3j. However, if p'Xl^) and Yl, have the same sign, then the conclusion regarding the sign 
of y'g can be made without the evaluation of sum. For example, for the attractive regime 
/X > 7r/2 we have 

^ _cosh/3 < 0, < 

p'm > 0, <l>'_i,,(/3) > 0, t>l. 

Hence, 

(J_i = -1, CTi, . . . , (Tt-i = 1. 

Let us compute the density of roots in the ground state in the attractive regime. We 
should choose appropriate regularization of the momentum, since Eqs. ( |1.15|) make no sense 
with p ~ sinh /5. In |T^, for example, the sharp momentum cutoff is chosen, i.e. p = 
for \P\ > A with some large A. Such a choice makes the function of momentum to be non- 
analytic, and disregards large momentum pseudoparticles at all. Although it works in the 
attractive regime fine, it gives some problems in the repulsive regime for /i < tt/S. We shall 
choose here a smooth cutoff, which comes naturally from the lattice regularization of the 
Thirring model (see Figure 1.2): 



p{P) = —imo In 



sinh !(/? + / — iuj) 
sinh |(/? + / + icu) 



— irriQ In 



sinh \ — f — ioj) 
sinh \ {l3 — f + iu}) 



;i.i6) 
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Figure 1.2: Smoothly regularized momentum 
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p'{[3) = rriQ smuj 



+ 



1.17) 



cosh(/5 + /) — cosu; cosh(/3 — /) — cos a; 

We reproduce the usual relativistic expression in the hmit f —>■ oo: 

p{P) ArriQ sin uje~^ sinh (3 = mi sinh /5, / oo. 

Denote the first term in ( |1.16| ) by /+(/?) and the second one by f-{P)- Then the regu- 
larized energy is 

hif3) = -Ui/3) + f4(3) (1.18) 
As expected, h{P) mi cosh/3 as / — oo. The Fourier image of (|1.17|) is 



piiik) 



e"'^p'_i{l3)dl3 = I e*'^V(/3 + in)d(3 



-Airm, 



sinh uk cos fk 



sinh nk 

We need to solve the Bethe equation for the ground state density p_i: 

This can be easily done by passing to the Fourrier transformed equation: 

^-1-1 • P-i, 



;i.i9) 



p-i 

271 



where 

Thus, we get 



Eventually, 



sinh(7r — 2uj)k 
sinh nk ' 



p-i 



^-1,-1 



P'-i 
27ri_i,_i 



2 sinhtufc cosh(7r — uj)k 
sinh tt/c 



mo 



cos fk 



cosh(7r — uo)k' 



;i.2o) 



j- I e-'^^p^i{k)dk 



mo 
4^ 



+ 



_cosh^(/3 + /) cosh|^(/?-/) 
Upon the cutoff removal, / — ^ oo, the density ( p..21| ) becomes 



P-iiP) 



m^e 



cosh — . 

/i 2fi 



(1.21) 



(1.22) 



We see that what we got looks like the density of free pseudoparticles in the Dirac sea, with 
the bare mass and bare rapidity renormalized as a result of interactions. Note that mass 
renormalization depends upon the regularization procedure chosen. 
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1.3 Thermodynamic Bethe ansatz 



The thermodynamic Bethe ansatz (TBA) method allows one to obtian the excitation energies 
and to compute the exact free energy. We employ it here for the Thirring model and discuss 
briefly its main concepts. The idea of TBA dates back to the work of C.N.Yang and C.P.Yang 
T5[| . They introduced the temperature into the Bethe ansatz technique and used Eqs. ( |1.15|) 



to realize the statistics of states at temperature T. 

The basic equation of the TBA method can be written in the form: 

£(p) = h{l3) -A+ * ln(l + e-=/^). (1.23) 

ZTT 

By virtue of its definition, the pseudoenergy of excitations 6{l3) is 

e^/^ = i (1.24) 
P 

Eq. ( [1.23| ) is the result of minimization of the free energy F = E — TS over p under the 
condition / p =const. Thus, the Lagrange multiplier A is just the chemical potential. The 



entropy of the system S is given in , and E is the total energy of all pseudoparticles. We 
note here that the expression for the entropy differs for the "fermionic" and "bosonic" cases. 
We work with the fermionic case where for each quantum number Ij corresponds only one 
pseudoparticle. 

Eq. ( |1.23[ ) is a non-linear integral equation that cannot be solved analyticaly in general. 
In the limit T — > Eq. ( |1.23|) reduces to 



e{p) = h{P)-A-^^'*e. (1.25) 



Integration in the latter formula is over all rapidities where e{j3) is negative. From ( |1.24 ) 



also follows that for such P we have p(/5) = 0, and, respectively, for the values of (3 where 
e{P) is positive, p{P) = 0. Let us represent e(/3) in the form 

£{(3) =£+ + £-= £H{e) + eH{-e), 
where H{e) is the Heaviside step function. Then from ( |1.25| ) follows: 



£+(/?) = h{(3) - A - (^Ij- + 5^ * e-. (1.26) 

In other words, > corresponds to the particle excitation energy, and to the hole 
excitation energy. e~ coincides with e where e is negative, i.e. in the Dirac sea. For example, 
if we take a pseudoparticle out of the Dirac sea and put it above the sea, the excitation energy 
would be 

Eext = eif3p) - eiPh) = £+(/5i) - ^-(/Js) > 0. 
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The total energy 

E = f hpdf3= f p{e + A+—<^' *e-)df3= (1.27) 
= J[ep + Ap+ —6- ■ ($' *p)]= J {e+p - E-p) - ^ 1 p'_^e- + A J p, 

where we used ( |1.25| ), ( |1.15| ) and the permutation property of the convolution operator. Eq. 
( 1.27 ) has clear physical interpretation: the energy of an arbitrary state at T = can be 
expanded as the sum E = E^ac + -Eezc+ const, where 

E^ac = -77- / P'~i^' < 



27r 

is vacuum energy, and 

Eexc = J {e^P - e'p) > 

is excitation energy. 

TBA equations for T = can be obtained also directly from the variation of total energy 



1% . Under small perturbations p, p, + Spj of the ground state densities p) the total 



energy changes by 

3 3 

Substituting 5pj from the Bethe equations ( |1.15| ) we obtian the basic spectral equations 

hj = + Ajk* cTkC'^. (1.28) 

Similarly, we deduce 

Pj = TC^ + Ajk * (Tk-Kj: . (1.29) 



From ( |1.28| ) it follows for the Thirring model 

iZi = (1.30) 

with 

/^^(fc) = -4v^.mo ""^t"^^^ • (1-31) 
smh tt/c 

We obtain 

27rmo sin/fc 



k cosh pk 

In the limit / ^ 00 we get for the energy of hole in the Thirring model: 

eZi -4moe-"^/2M ^03^ (1.32) 

The hole of the Thirring model corresponds to a soliton of the sine-Gordon model. For the 
rest of the excitations we have 



^ ^3 ^i.-i^-i ^ ^3 ^ h_i, (1.33) 
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where we have substituted ( p.. 301 ) to get the latter. Using the Fourier transforms Ajk together 
with 

sinh(7r — juj)ksm fk 



h'Ak) = AmrriQ 



sinh nk 



we obtain 



and 



sm 



(tt - /i)7rj 
2/i 



cosh 



2^' 



j<t-l 



eUiP) = -e-,{f3). 

t — 1 string in the Thirring model correspond to antisoliton of the sine-Gordon model, while 
the other j-strings correspond to the breathers. 



1.4 Scattering matrices 

The method to compute elastic S-matrices from the Bethe ansatz equations was developed 



m 



We describe the basic idea briefly, taking as an example scattering of holes (solitons 
in the sine-Gordon model). 

By definition, the phase shift for scattering of two holes is given by 

5u{l3i -132) = - log S = vr- y^a, (1.34) 

where ipi is the phase gained by a hole when going around the system and the same phase 
but in the presence of another hole. The 2 are composed of a sum of two-particle bare 
phase shifts; for example 

y.2 = i:p-i(/9i) + E*i,i(A-/3,) 

Here the sum is taken over all the solutions (3j of the Bethe equations ( |1.8D with two holes 
at positions (3i and /?2- The result of subtraction Lpi — ip2 is proportional to the backflow 
function of vacuum: 6h = 27rF{f32\Pi)- The latter is defined using the difference of the two 
solutions of the Bethe equations obtained with and without the hole at Pq, 

F{(3o\f3)^{(3-P)Lp{(3). 

One can show using ( p.. 81) that F satisfies the following integral equation: 

- /5o) = $ * F + 27rF (1.35) 

Equation (|1.35|) describes the backflow caused by a hole at (3 = Pq on lm(3 = it axis. Taking 
a derivative with respect to (3 and applying the Fourier transform to both sides of ( |1.35|) we 
arrive at the following solution (in Fourier space): 

F'ik) = f^/""'^^^ (1.36) 
2 cosh fj,k smh(7r — p)k 
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From (|1.34|) we obtain: 



li- log sm = r Tl^y (1-37) 

t dO J-oo 2 cosh fik smh[7i — fi)k 

where 6 = (3i — f32- Note that this method does not fix the constant normahzation factor in 
the matrix element, which should be fixed by other constraints (e.g. unitarity). Expression 
( |1.37| ) is in agreement with the results of , which enables us to identify coupling constants 



in the Thirring model with those of the SG model: /i = vr — Psq/S. 

The conserved charge equals to the total number of pseudoparticles Q = J pd(3. The 
physical (renormalized) charge is obtained after the subtraction of the charge of vacuum. It 
can be easily calculated from the Bethe equations using the backfiow functions. E.g., for a 
hole Qhoie = — 1 + -^'(0), while for n-string Q„ = n + -F^(O). Thus, we obtain the following 
values: Qhoie = — vr/2a; for the hole, Qt-i = vr/2cc; for the t — 1 string, and Qn = for the 
other strings. 



1.5 Remarks 

There are many interesting issues of the Bethe ansatz technology left beyond the scope of 
this chapter. We want merely to mention some of them here. 

Somewhat more accurate analysis of the permissible solutions of the Bethe equations for 
the XXZ chain, based on the counting arguments, shows that there are certain constraints 
on the space of physical states. For example, the holes can exists only in pairs in this model 
T^ . Another interesting example is so-called imaginary mass Thirring model discussed in 



19], where the pseudoparticles are paired in the Dirac sea. 



The important questions are the norm of the Bethe states ||20[, and their completeness 
|21[] . The latter issue for the XXZ chain is equivalent to showing that the number of the 
physical states obtained from the Bethe ansatz equations equals to 2^, where N is the 
number of sites on the chain. 

Finally, it is worth to mention that the thermodynamic Bethe ansatz technique can be 
applied also to the physical excitations, and in this context it allows to find the exact free 
energy as a function of scaling parameter mL, as well as the central charges of the fixed 
point theories [§. 



19 



Chapter 2 



Classical and semi-classical analysis of 
the boundary sine-Gordon model 

We consider the sine-Gordon model on a half-line, with an additional potential term of the 
form — Mcos|(v9 — cpo) at the boundary. We construct the classical solutions in the next 
section by using the bulk sine-Gordon theory and the "generalized method of images." From 
the classical solutions in hand we extract the time delay ( |2.18D as follows. We send a soliton 
(anti-soliton) which lives in the semi-infinite world governed by the boundary sine-Gordon 
model from some large position xq at time to and measure the time ti it takes to bounce 
off the boundary and come back to xq. At the same time to in some other, infinite, world 
governed by the bulk sine-Gordon model we send a soliton with the same speed but in the 
opposite direction from the position —xq and measure the time t2 it takes to arrive at Xq. 
Now, At = ti — ^2 is our time delay. From the time delay one can reconstruct the semi- 
classcial phase-shift using the method of Jackiw and Woo ||22|. We establish the agreement 
with the semi-classical limit — *• of the exact boundary reflection matrix, (|2.34|) . The 



exact expressions for the boundary reflection matrices are known up to CDD ambiguities 
|l|]. They were obtained as a "minimal" solution to the general set of constraints for the 



integrable boundary field theory. This chapter is based on 23 



2.1 The boundary sine-Gordon model 

The Lagrangian of the boundary sine-Gordon model is given by: 

dx + Mcos^{(p{x = 0)-ipo) (2.1) 

Vanishing of the variation of Lagrangian under a small perturbation (f ^ (f + 6if is 
equivalent to two conditions for the field (f. One is the sine-Gordon equation on the interval 

[0,oo) 

(f>tt - (f^xx = -sin{(t)), (2.2) 



SG 



1 r+oo 

2 Jo 



{dt'^f - {d.^'^f + ^ cos/3(^ 
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where (3ip = 0, whereas another is the boundary condition 

5x0|x=o = M sin ^ (0 - 0o) L=o • (2-3) 

The condition ( p.3|) is the most generic boundary condition preserving the integrabihty of 
the bulk problem It has two arbitrary paramenters, boundary mass M and a phase 0o- 



2.2 The classical solutions 



The initial-boundary value problems compatible with integrabihty for the classical non- 
linear equations is a field itself withing the non-linear science. Such problems were studied 
extensively by mathematicians. One of the first results on the initial-boundary value problem 
for the classical sine-Gordon equation were obtained by Sklyanin In ||2^ the particular 
form of the most generic boundary condition was reported, dxcf) = Msin(0/2) at x = 0. 
Later, this condition was widened in |25| to include the Dirichlet (j){x = 0) = const and dxip = 



M cos(0/2) cases. Some solutions to the posed boundary value problems were obtained using 
the inverse scattering method and the Backlund transformations |2^, |2^. Still, the most 
generic boundary condition was missed in |2^. The condition ( p.3| ) first appeared in the 
physics litearure in the seminal work of Ghoshal and Zamolodchikov 

In this section we construct explicitly the solution to (^.1D-(P75D that replicates the soliton 
(antisoliton) scattering from a boundary, as well as the solution localized at the boundary, 
called boundary breather. The idea that we use to construct such solutions is a remniscent of 
the method of images (and therefore we refer to it as the "generalized method of images." ) We 
employ the known solutions on the full line - kinks and anti-kinks. To provide the reflection 
of the kink from a boundary (and thus the correct assymptotic behavior at infinity), we send 
another kink of the same kind towards it - the "mirror image" . Such a solution automatically 
satisfies (p.l|) . Besides, in order to satisfy also (|2.3|) , we place the third, stationary kink at 
the origin. We check by a direct calculation that condition (p.3|) is satisfied. To obtain the 
three-kink solution, we need the knowledge of the r-function for sine-Gordon equation. 



2.2.1 The r-functions 

On the infinite interval, (—00,00), the classical multi-soliton solution to the sine-Gordon 
equation is well known It is usually expressed as: 



0(x, t) = 4 arg(r) 



flmir 
4 art an I 



7^e( 



r 



where the r-function for the A^-soliton solution is: 



r 



Mj=0,l 



exp 



N 



kj 



+ 2 /^*/^.- log (y-ry 

l<i<j<N V^i -r i^j 



(2.4) 



(2.5) 
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The parameters kj, aj and ej have the following interpretations. The velocity of the j^^ 
soliton is given by: 

- (tn) ■ 

(Note that vj is positive for a left-moving soliton.) The aj represent the initial positions 
of each of the solitons, and ej = +1 if the j^^ soliton is a kink, while ej = —1 if it is an 
anti-kink. The rapidity, 6, of the soliton is defined hj k = e^, and we have normalized the 
soliton masses to unity (in further discussion, the words "soliton" and "kink" will be used 
synonymously) . 

It is fairly obvious how to get a single soliton solution on [0, oo) with either 4'\x=o = or 



dx(l)\x=o = 0. One exploits the symmetry of (|2.2| ) under (f) —(f) and x — > — x, and simply 



takes a two soliton solution on (— cxd, oo) where one soliton is a mirror image of the other 



through X = |^. If one does this with a double-kink solution then it satifies the foregoing 
Dirichlet condition, while the kink-anti-kink solution satisfies the Neumann condition. It is 
also not hard to guess how one can go beyond this solution: For M = oo, the boundary 
condition ( p.3|) reduces to 0|a;=o = 0o- The only way that this can be obtained from a 
multi-soliton solution on (— oo, oo) is to put a third, stationary soliton at the origin. 

We therefore consider the three soliton solution with ki = k, k2 = 1/k and k^ = 1. That 
is, we consider: 



k-1 



2 




r = 1 - et;2e-^(^ +1)^-" - eo e-^^'^+i)^"" F(t) 

where we have introduced the shorthand: 

e = €162 , eo = €3 , a = ai + 02 , b = . (2.8) 

The function F{t) is defined by: 

F{t) = eie-^('='-^)*-"^ + eae^^'^'-^)*-''^ . (2.9) 

This solution has = at x = 00, and for > 1 it describes a left-moving soliton moving 
from X = 00 with a right-moving "image" starting at x = —00. There is a stationary soliton 
with center located at x = —b. Viewing this as scattering off a boundary at x = one 
can easily see that a is the phase delay of the returned soliton. To make this more explicit, 
observe that r has the following asymptotic behaviour: 



r(x,t) ~ 1 + ,e e-*[(^'+^)"+(*^'-^)*l-'^i x, -t ^ 00 with 



X 



t F + l 



r(x,t) ~ 1 + ,e e-^[(' x, t ^ 00, with - = 4 



P-1 



t k^ + 1 ' 
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The problem now is to first show that the r-function given by ( |2.7| ) provides a solution to 
the boundary value problem on [0, oo) defined by ( |2.2| ) and ( |2.3| ). Our second purpose is to 
relate the parameters a and b of ( p.7|) to the parameters M and 0o of ( p.3|) , thereby obtaining 
the classical phase delay, a, in terms of M and 0o- 

2.2.2 The classical phase delay for M = oo case 



In the center of mass reference frame the solution to ( p.2| ) obtained by means of the r-function 
method reads: 

2e^'^^^~"'^c]l[ts]l9) lb e^~"^ =P ga;(l+2ch0)-2ai-a3+21og(u2ii) 
= =F4 arctan 2gx(l+che)-ai-a3+21og(«)p}^(^^gJ^5)-) g2xch0-2ai+21og(i>) :jp I' (2.10) 

where m = tanh(|), v = tanh(^), and we have set Oi = 02, which means that the solution is 
invariant under the transformation t — ^ — t. The upper (resp. lower) sign corresponds to the 
situation when a stationary soliton (resp. anti-soliton) is used to adjust the value of field at 
the boundary. Solution ( |2.10| ) refers to the case when the incoming and outgoing particle is 
the soliton with asymptotic value = 27r at x = +00. 

Let us represent ( ^.101 ) in the form of rational function of variable ch.{tsh.6): 



rs + S2ch(tsh^) 

0(x, t) = 4 arctan ■ . 

ri + sich[tshu) 

The condition (f){x = 0,t) = 0o implies that 

r2 S2 




tan — , (2.11) 

ri si 



from which follows immediately 



03 = log |^=F'W^ tan —j . 

For the argument of logarithm to be positive one should take 0o > with the stationary 
anti-soliton and 0o < with the stationary soliton. This is illustrated in figure 2.1. Further, 
we obtain from (|2.11|) 



2ai = log 



■ 2 2 M' + tan2(^^) 



(2.12) 



Note that the time delay, obtained by ( 2.12 ), is in fact always a time advance in both the 
attractive and repulsive cases. For the same value of 0o the time delay for the soliton that 
lives on the left half-line a; < is not the same as that of the right half-line soliton (except 
for 00 = ±7r). The position of the "left" soliton is not the exact mirror image of the "right" 
soliton for generic 0o- 
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2.2.3 The classical phase delay for generic case 

To summarize the computation briefly, one substitutes ( p.7| ) into ( p.3| ), and obtains the 
constraint: 



\Re{T)dncXm{T) - 7^e(r)9^Tm(r)]^^o = 

1 1 
M 2cos(-0o) ^e(r)Jm(r) - sin(-0o) (ne{r)^ - Im{T) 



(2.13) 



x=0 



One can solve this by brute force substitution for the real and imaginary parts of r, but it 
is somewhat simpler to find constants a, /3, 7 and b such that: 



\a ne{T) + (3 Im{T)] , 
[7 7^e(r) + 5 Xm(r)] 1^^^ . 



(2.14) 



One then finds that ( p.l3|) can be satisfied if and only if one has a = 6, which indeed turns 
out to be true. Using this one arrives at: 



M cos(-0o) 



(P + 1) 



M sin(-0o) 



1 + e v'^e-'' 



-2b 



'k-1 

k + 1 



1 + ev' 



'k-r 

k+ 1 



(2.15) 



where 



A 



+ e 



-2b 



'k-V 



e V 



(2.16) 



It is algebraically very tedious to invert this relationship. One proceeds by eliminating 
e~^, and then solving for a. It is very convenient to introduce a new parametrization of M 
and 4>n: 



jj, = M cos(-0o) 



u = M sin(-0o) 



2 cosh(C) cos(r7) 
2 sinh((^) sin(r7) 



(2.17) 



where < ( < 00 and —n < r] < n. (In the (/i, i/)-plane the curves of constant ( are ellipses, 
while the curves of constant rj are hyperbolae whose asymptotes make an angle of |0o with 
the /X-axis.) 
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Figure 2.1: a solution of the classical Sine- Gordon equation with the fixed boundary con- 
ditions. For (j)\x=o — <I>Q — ^ solution is constructed out of a left-moving soliton, its right- 
moving image at x < and the stationary antisoliton in the middle (upper graph). For 
4>\x=o = 4>Q = —Y' solution is built out of three solitons (lower graph). 
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We then find that the phase delay is given by: 



a = Iog<^ -e tanh2(e/2)tanh2^ 



tanh 1(6' + 277) tanh ^{9 — irj) 



tanh |( 



+ tanhi(0-C) 
There are several things to note about this formula. 



(2.18) 



(i) The ambiguity of the ±1 power comes from solving a quadratic equation, and is a direct 

reflection of the fact that ( p.3|) is not invariant under — >■ + 27r (whereas ( p.2|) is 
invariant under this shift). Equivalently, one can flip between the + and — powers by 
sending 0o 0o + 27r. 

(ii) The argument of the log is always real and positive. The discrete parameter, e = £162, 

must be chosen to arrange this. Hence: 



e = +1 for - C < ^ < C 



-1 for 1^1 > C 



(2.19) 



This means that a kink reflects into kink for — C < 9 < and reflects into an anti- 
kink for 1^1 > (. This is consistent with the fact that Dirichlet boundary conditions 
(M = 00) cause a kink to reflect as a kink, whereas Neumann boundary conditions 
(M = 0) cause a kink to reflect as an anti-kink. Note that these two domains of 
parameter space (in which a kink reflects differently) are separated from one another 
by a logarithmic singularity in the classical phase delay. 



(iii) The choice of the power ±1 in ( p.l8| ) is correlated with the parameter rj and whether 
there is a kink, or anti-kink at the origin. Specifically, we have: 



± sign{tan{-)) 



(2.20) 



where the ± is the same as that of ( p.l8| ). 
(iv) In the M — > 00, or Dirichlet, limit we see that: 

1 , 



C ~ log(M) 



V 



and 



(2.21) 



and the phase delay collapses to ( p.l2| ): 



log { tanh^ l^j tanh^(6') 



tanh ■ 



+ i- 



tanh ■ 



— t- 



±1' 



(2.22) 



In this limit (p.3|) enforces Dirichlet boundary conditions. It is, however, important 
to note that there are two possible distinct boundary values: (f)\x=o = (f>o and (f)\x=o = 
00 + 271. Since the boundary potential is — Mcos(|(0 — ^o)), one sees that (f)\x=Q = 4>o 
is stable, while (j)\x=o = 0o + 27r is unstable. 
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From now on we consider only the stable boundary value that corresponds to the positive 
sign in (|2.22| ). Then one has from (|2.2(]| ) 

eo = - sign(tan{^(l)\^=o)) . (2.23) 

It is essential to observe that we have taken (f)x=oo = 0, ab initio. For different boundary 
conditions at x = oo, one should replace 0o by (j)\x=o — 0|x=oo- This physical parameter is 
defined mod Att. 



2.2.4 Boundary breather solutions 

To fully understand the semi-classical scattering computation one also needs another class 
of classical solutions, which we call here boundary breathers. It is well-known that breathers 
represent bound states in the soliton-anti-soliton channel in the bulk sine-Gordon theory. In 
the same way that the classical bulk breathers can be obtained from the appropriate solution 
by analytic continuation of 6 to imaginary axis, one might expect that the same procedure 
would give boundary breathers on the half-line. To see this, we set 6 = 0} {0 < i} < |) in 
the formula (^A0|) . Next, we impose the following conditions for a solution to be boundary 
breather: 



a) it should be a real function, 

b) it should have finite energy and 

c) the asymptotic value at x = +oo must be fixed and equal to 2?™ {n - integer number). 

The three-soliton (resp. soliton-anti-soliton-soliton) configuration satisfies the first con- 
dition provided that ^9 < — ^ (resp. < ^). However, the other conditions are satisfied 
by the soliton-anti-soliton-soliton configuration only, which one could have foreseen from the 
analogy with the bulk theory. The boundary breather solution has the form (see figure 2.2): 



= 2tt- (2.24) 
—4 arctan ■ 



2cot^9cot|/^cot^e^+^"°'^''cos(tsin^9) + e2^"°"''/s:cot2| + 1 



2 cot ^ cot f VK ™^ ^ cos(t sin ^) + cot f (e^^ ™^ ^iT cot^ | ) 
where 

ir = cot(— --)cot(— + -). 

So, we have a continuum of classical boundary bound states when < < ^ < | and 
for other (pQ according to the 27r-periodicity. In the quantum theory this continuum shrinks 
into a discrete set of bound states (see next section). Note that in the limit -i? — > ^ the 



boundary breather ( |2.24| ) reduces to the ground state, figure 2.3, and the phase delay has 
singularities at 6* = ±|«0o- An analogous picture of bound states occurs for the anti-soliton 
scattering with fixed boundary conditions. 
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2.2.5 General solutions, integrability and 
Backlund transformations 

Thus far we have only apphed the method of images to obtain certain special classical 
solutions of the boundary sine-Gordon problem It is natural to suggest that general 

solutions of ( p. ID can be obtained by similar methods. This in turn would establish the 
classical integrability the boundary problem ( |2.1|) . It is, in fact, rather straightforward to 
show that both of these conclusions are true, at least for the problem (|2.1| ) with (pQ = 0. The 
method we will employ can almost certainly be generalized to problems with 0o 7^ 0, and 
also has the virtue that it can be used to construct the integrable boundary potentials for 
the more general Toda models. A related approach has been followed in |2^, |3y]. The basic 
idea is to use the fact that any integrable hierarchy has Backlund transformations: that is, 
solutions can be mapped into one another by non-trivial gauge transformations constructed 
from the affine Lie algebra action on the corresponding LAX system. For the sine-Gordon 
equation, the requisite Backlund transformation can be cast in the following form: 



where u = x — t, v = x + t, and ( is an arbitrary constant parameter. The point is that 
satisfies the sine- Gordon equation (|2.2|) if and only if ip does so as well. Suppose that is 
a solution to sine-Gordon on [0, oo) satisfying a Dirichlet boundary condition: ip\x=o = 2^7, 
where ?7 is a constant. It follows immediately from ( p.25|) that is a solution to sine-Gordon 
satisfying ( |2.3| ), where M and 0o are given in terms of ( and rj by ( p.l7| ). Thus, if one 
can solve the Dirichlet problem, one can solve the more general problem by a Backlund 
transformation. 

Observe that if = 0, or equivalently 0o = 0, then the Dirichlet problem can be solved 
trivially by method of images: one gets the solution on the half-line by extending it as an 
odd function on the full line. Thus, the Backlund transformation essentially defines the 
generalized method of images. It is also by no means an accident that the parameters 
entering into the Backlund transformation [t] and () are precisely the rapidity parameters 
that turn up in the phase delay (|2.18|) . One should also note that the form of the integrable 
boundary potential is given directly by the Backlund transformation. This fact should easily 
generalize to Toda systems. 

Backlund transformations, in general, are invertible transformations on the solution space 
of an integrable hierarchy. The simplest forms of them modify the constants of the motion 
of solution, and possibily add or subtract a soliton. One can certainly find a Backlund 
transformation that will introduce a stationary soliton into the soliton-soliton solution of 
sine-Gordon. As a result, the general three soliton solution employed above can be obtained 
from the two soliton solution that is appropriate for the "trivial" Dirichlet problem with 
00 = 0. We therefore expect that any solution of the trivial Dirichlet problem can be 
mapped onto the generic problem thus establishing the classical integrability. 




(2.25) 
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Figure 2.2: a boundary bound state (boundary breather) for the values 6 — ^, (f)o — ^ Sit 
four distinct instants of time. 
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Figure 2.3: two possible ground state configurations for 0o = The configuration with 
asymptotic behaviour (^6 — > at infinity has lower energy than the other one. 
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2.3 The semi-classical analysis 

2.3.1 The exact boundary reflection matrix 

Let us recall some of the results of The exact (quantum) boundary S matrix elements 
are 

s+{0) = p+{e) = cos[^-{t-i)ie]R{e) 

sz{0) = P-{e) = cos[^ + {t-i)ie]R{e), (2.26) 

and 

s- = Q+{e) = s+ = Q4e) = ^ sin [2{t - i)te] R{e) , (2.27) 

where 6 is the rapidity of an incoming particle. The parameter t is defined by: 

Sir 



The function R{9) decomposes as 



R(9) = Ro{9)Ri{9), (2.28) 



where Ro is a normalization factor ensuring unitarity and crossing symmetry that does not 
depend on the boundary conditions. The dependence upon boundary conditions appears in 
Ri, which reads 

Ri{9) = ^a{r],9)a{te,9). (2.29) 
cos 4 

Two of the four parameters k, ^, rj, B are independent, and we have the relations 

cosr] cosh 9 = — cos^, cos^ + cosh^ 6 = 1 + — . (2.30) 

These parameters are related with M and ipo in an unknown way. An expression for the 
functions Rq and cr involving infinite products of F-functions is given below (see next chapter). 
There are also simple integral representations: 

^(^.^g) ^"^^ cxpfr - ^'""Ht - 1 + ^ait-DoA (2.31) 

cos[^ — i{t — 1)9] X 2cosh(t — l)xsinha; J' 



Ro{e)=e.J-r ^ sinh|^-l)xsi^h(|-l)x ^^.32) 



and 

RJft\ = PVT. J - 

X sinh I sinh 2(t — l)x 

One has cr{C,,9) = a{—^,9). The only difference between the scattering of solitons and 
anti-solitons arises therefore from the pre-factor cos[^ =F i{t — 1)^] in 
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For simplicity we consider only the limit in which M — > oo. One can identify the 
corresponding values k = and Q = oo easily since, at these values, the topological charge 



f3 f°° 
Q = — d^ipdx 
2tt Jo 



is conserved and therefore the amplitudes Q± must vanish. One has also 1] = ^ so 

Rm = ^a{^,9). (2.33) 
cos 

Consider now the quantum boundary S matrix at the leading order in ^ as — > 0. The 
computation is most easily done by using the integral representation given above, and eval- 
uating the integrals explicitly by the residue theorem. This provides convergent expressions 
where the /3 — > limit can be investigated term by term. To get non-trivial results one must 
scale ^ as 



0^ 



We find then 



where 




S{e; y)=exp{^ I dv In tanh^ ^^^^^ ) , (2.35) 



5(6*; 0) being the semi-classical approximation to the bulk soliton-soliton S-matrix ||22| (in 
the following we denote it also as S{d)), and /t = sign^ (in the following we assume that 
e > 0). 



Before discussing the relation between formula ( 2.34 ) and the classical computations of 
the preceding section, it is useful to comment on the bound states of the quantum theory. 
Poles of the Rq term are associated with breathers and do not correspond to the boundary 
(new) bound states. The latter correspond rather to poles of the a{^, 6) term which are 
located in the physical strip Im^ e [0,7r/2]. By inspection of the F-product expression 
one finds two families of poles of a: 



t-l\ 2/ t 

where n, I are integers. Let us restrict to ^ > 0. Then the only physical poles are those that 
correspond to + sign in 9^^^ and — sign in 9^'^\ The first pole that enters the physical strip 
(from the bottom) is ^r C ^ f • The number of poles of a in the physical strip increases 
monotonically with ^ for ^ small enough, and as t gets large it becomes simply of the form 
^/tt. These poles are cancelled by the cosine term in P_ and therefore appear only in the P+ 
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amplitude. As C, reaches the value ^ = ^ the poles densely fill the interval < Im^ < | and 
a pole at 6 = in /2 appears corresponding to the emission of a zero momentum soliton. As 
argued in this corresponds to a change in the ground state of the system. For < v?o < ^ 
the ground state is ^ at infinity but for | < (/?o < 3| it is — ^ at infinity (see figure 
2.3). Therefore the value ^ = ^ corresponds to (fo = ^ and it is the upper physical value for 

^: as (po varies arbitrarily, ^ never gets larger than ^ and the set of poles 6^^] never enters 
the physical strip. Observe that there are bound states in the quantum theory when in the 
classical case the kink sitting in the middle and the incoming one are of opposite topological 
charges in complete agreement with the discussion of boundary breathers in section 2.1.4. 



2.3.2 Classical time delay and quantum phase shift 



To establish the relation between (|2.34| ) and ( p.22| ), first recall |2^ the general relation 
between the quasi-classical scattering phase shift, S{9), and the classical phase shift, a{9) : 



6{6) = constant H / aiifjdrj, 

2 JO 



Here, m is the classical mass of the particles involved in the scattering. Using the 
classical Levinson theorem" to determine the constant in this formula we deduce, 
(|]22]), the relation 



where 



26{e) = 2nB'n + -| / 



drj log 



tanh^ rf tanh^ ^ 



tanhi(r/ + i^)tanhi(r/ 



(2.37) 

"semi- 
using 

(2.38) 
(2.39) 



and is the number of bound states. Since according to the preceeding discussion ub = 
^ = || for < ,^ < ^ as /3 ^ we see that formula ( p.34| ) is in complete agreement with 

the classical phase shift (|2.22|) and that ^ is proportional to 0o- In the region — ^ < ,^ < 
there are no physical poles for P+ and ra^ = 0. A similar discussion can be carried out for 
< and P_. As /5 ^ the number of physical poles of P_ varies hb = = and 
again we have agreement with the semi-classical Levinson theorem. 

From comparison of ( |2.34| ) and the classical phase shift we see that and (^q are related 
linearly as ^ = ^v^q. This leads correctly to the emission of a zero momentum soliton 
with the ground state degeneracy as discussed in This linear relation can hold only for 
\i\ < Beyond that value the ground state changes. To compare the quantum result 
with the classical phase shift we must then correlate appropriately the value of at infinity 
in the latter. The net effect is to replace 0o by 0o ~ 27r. Eventually, the variation of ^ with 
ifQ is therefore 




/3v?o 
27r 



1 

+ 2 



(2.40) 



as illustrated in figure 2.4. It is very likely that (|2.40|) is exactly true for finite /3 as well. 
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Figure 2.4: Variation of as a function of (po. 
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We can finally recover ( |2.34| ) without appealing to our knowledge of quantum boundary 
bound states. For this one has to evaluate the action for the classical configuration. Following 



the discussion in [22| we have, 



or/m ^ ^ , tanh^ r/ tanh^ 2 



if \P^o\<7i, (2.41) 



where we used the soliton mass m = and 6 satisfies the differential equation 



r+oo / r+oo \ 

6{9) - tanh9 6'{9) = J dtij rfxyi^ - 8 sinh ^ tanh ^ j . (2.42) 

We restrict to the situation where we send in a kink, and (po is positive and so there is 
an anti-kink at the origin (see figure 4). In this case there are quantum boundary bound 
states. It is difficult to perform the double integral for the three-soliton solution explicitly 
because of the cumbersome expression for the integrand ( p.lO|) . One might hope that in 
order to find the ^-independent piece of the phase shift it is sufficient to evaluate both 
sides of ( p.42| ) in the limit as ^ — > -|-oo or ^ — 0. However both of these limits do not 
seem to be helpful, because due to the non-uniform convergence of the right hand side of 
( p.42| ) it is not allowed to interchange such a limit with the integration. We evaluated the 
right hand side of ( p.42|) (with fixed 9) numerically for several different values of 00 using 
Mathematica. Combining ( |2.41| ) and ( |2.42 ), we obtain an estimate for C((y9o); in agreement 
with C(v?o) = yV^o with accuracy 0.1%. This result was checked for different values of 9. 
We therefore obtain agreement with the semi-classical Levinson theorem. 



2.4 Remarks 

Although the method of images works nicely in the classical theory, it does not seem to 
extend to the quantum case: we have not been able to recast the boundary S matrix of 
as a product of bulk S-matrix elements. A related phenomenon is the non-trivial structure 
of the boundary S-matrix with one-loop corrections (the semi-classical expressions being the 
tree approximation). Recall that in the bulk one has 



Si{9;y) = exp ( — j (if Intanh^ 



V + iy 



where 



87r-/32 t-l 
Scaling 
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one finds the following next-to-leading correction to the boundary S-matrix in the Dirichlet 
problem: 



P±{0) 



exp ±K ^.^ + K, 



1/2 



72 



12 



1/2 



In addition to the usual replacement of by (3'"^ we see the appearance of an entirely new 
factor involving the square-root of a hyperbolic tangent. 

The massless limit mo = of the boundary sine-Gordon model ( p^.lD was studied in pi 
where, in particular, the massless classical solutions were obtained. 



35 



Chapter 3 

Non-relativistic limit of the quantum 
sine-Gordon model with Dirichlet 
boundary condition 



In the previous chapter we studied the quantum sine-Gordon model on a half-hne ( |2.1|) 
in the semi-classical limit (3 —>■ 0. The purpose of this chapter is to investigate the non- 
relativistic, 6—^0, limit of this model and to determine the quantum-mechanical potential 
induced by the presence of boundary. We show that the generalized Calogero-Moser model 
with boundary potential of the Poschl- Teller type describes the non-relativistic limit of (|2.1|). 



The discussion is based on 32 



3.1 Introduction 

When 6^0, the energy and momentum of relativistic particles - in our context solitons 
and anti-solitons - degenerates to 

E = cosh e 



p = Ms sinh e Ms9. 

So, 6* is a speed of particles measured in the units of c (the speed of light). 

Since such properties as integrability and factorized scattering survive in the above limit, 
we expect to get some quantum-mechanical integrable model on a half-line. Moreover, since 



a large class of such models based on the Lie-algrebraic classification is known |Q, it is 
natural to look first for the appropriate candidate withing this class. Indeed, we show 
that the boundary Calogero-Moser model with boundary potential of the Poschl- Teller type, 



whose integrability was established in p4| , p5| , describes the non-relativistic limit of (|2T 



Let us give an idea how one can come up with the boundary Calogero-Moser model by 
means of the heuristic arguments. It was known since long ago |36[ that the solitons 
and anti-solitons of the bulk sine-Gordon model interact via sinh~^(x) (particles of the same 
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kind) and cosh^^(x) (particles of different kinds) potentials in the non-relativistic limit. Q 
Now, the question is to find out the form of the potential induced by the boundary. But 
from the semi-classical picture of the previous chapter we infer that the "boundary" can 
be represented by a stationary kink or an anti-kink and a moving "mirror image". So, 
we obtain the superposition of cosh~^(x) and sinh~^(x) potentials as a boundary potential, 
i.e. the Poschl- Teller potential The sinh~^(a;) term is a hard-core reflecting potential. 



whereas the term — cosh (x) is necessary to provide the existence of the boundary bound 
states. 



3.2 The exact quantum field theory solution 

In this chapter we consider the sine-Gordon model on a half-line ( p.l|) , with the fixed value 
of field at the boundary: ip{x = 0,t) = (po, or M = cx) in (|2.1|). In [|^], the quantum 
integrability and the exact S-matrix were conjectured for ( |2.1| ). The boundary scattering 
matrix is diagonal and, according to [IH, the reflection amplitude of the sohton P+ (resp. P_ 
for anti-soliton) reads: 

P^(e) = cos{^ ±Xu)R{u,^) = cos{^±Xu)Ro{u)Ri{u,^), (3.1) 

where 9 = iu is the rapidity, ^ = ^V'o and A = || — 1; 

r(l-^)r(A + ^) ^ r(4AA;-^) 

"^'^"^ = r(i + i)r(A-i)nr(4AA: + i) 



r 


(l + 4AA; - ^] 


r 


(A(4A;+1) + ^] 


r(l + A(4fc-l) + ^) 


r 


[l + 4AA; + ^] 


r 


^A(4A; + 1) - ^) 


r (l + A(4A; - 1) - ^) 



(3.2) 



^ ifi r(| + 2/A + ^)r(i + 2/A + ^) 
^ l\ r (I + 2/A + A + ^) r (i + 2/A + A + Si 
r (I + 2/A + A + ^) r (i + 2/A + A - ^) 



r (i + 2/A + 2A + ^) r (i + 2/A + 2A 



(3.3) 



The poles of P± located in the physical domain 0<M<7r/2atMn = ±^ — ^fx^^r correspond 
to the "boundary" bound states of the theory. The latter exist in the soliton (resp. anti- 
soliton) scattering channel if ^ > (resp. ^ < 0), and their energy is En = MgCosUn- 
Note that the "physical" values of the parameter ^ are bounded: |^| < 47r^//5^ (see previous 
chapter). In the semi-classical limit of the quantum field theory ( p.l|) , /3 ^ 0, the principal 
("tree") approximation to the amplitudes ( |3.1| ) has the form ( p.34| ). 

^The quantum breathers correspond to the bound states in the — cosh^^(a;) potential. 
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3.3 The Calogero-Moser model on a half-line 



Our purpose is to show that the non-relativistic dynamics of quantum sine-Gordon sohtons 
in the presence of a boundary is described by the generahzed Calogero-Moser Hamiltonian: 



H 



1 ^ (P 



M 



N 



^Y.-^ + Y.^^AA{Xi-Xi:)+VAA{Xi + x'i)) 



i<i' 

M N M 

+ i^AA iVj - Vf) + Vaa iVj + Vj')) + YY i^AA {Xi - Vj) + Vaa {Xi + Vj)) 

j<j' «=i i=i 

N M 
i=l j=l 



(3.4) 



Here Vaa and Vaa are bulk nonrelativistic potentials obtained long ago in ^ 



Vaa (x) 



"o pip 



Ms sinh aox ' 



Vaa {x) 



Mg cosh^ uqX ' 



with 



P = lrl 



Stt 

J' 
mo 
2 ' 



(3.5) 

(3.6) 
(3.7) 



and Wa and Wa are boundary potentials of the Poschl- Teller |^ type 



Wa{x) 
Wa{x) 



at, 



/i(/i-l) z/(z/-l) 



2Ms ysinh aQX cosh a^x ^ 



an 



2Ms \ sinh oqx cosh oqx 



/i > 1, z/ > 1. 



(3.8) 



Let us comment on the properties of the one-particle Schrodinger equation with the 
Poschl- Teller potential Wa{x). The energy of the bound states, which appear when u > 

2 

is given by En = —^^iy — /i — 1 — 2n)^, where n = 0, 1, 2... For a fixed value of z/ — /i there 



are in total 



obtained to be equal to 



2 



bound states. The reflection coefficient, which is a pure phase, can be 



SA{k) 



V ( ih.\ F f 1 _|_ /'^'^ tfc \ p f fj.+u _ ik \ 

^ \ao) ^ 2 2ao) \ 2 2ao ) 

-p f ik \ "P / 1 _|_ M— _|_ ik \ -p / pL+u _|_ ifc \ 

^ V ao/ V2 2 2ao/ ^ V 2 20:0/ 



(3.9) 



This expression has "physical" poles on the upper imaginary half-axis in the complex mo- 
mentum plane which correspond to the bound states. Besides, it has poles at the points 
kn = (1 -|- n)ao that come from the first F- function in the numerator of (|3.9| ). The latter set 
of poles is infinite and does not correspond to any bound states of the theory. The S'-matrix 
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for the potential can be obtained from ( |3.9[ ) by the substitution /i z/. One can see 
that Sa and satisfy the boundary Yang-Baxter equation of 

Sa cos ^"^(^^ ~ lA ~ = cos {^{^ ~ /^) + -^""^ ■ (3.10) 

3.4 Relativistic vs non-relativistic 
integrable models 

To estabhsh the equivalence we will show in particular that the S'-matrices of the quantum 
sine-Gordon theory and the model (|3.4| ) coincide in the appropriate limit. The system ( p. 41) 



is integrable both at the classical and quantum levels To see this one takes the 

hyperbohc-type Calogero-Moser Hamiltonian for A^ + M particles based on the BCn+m root 
system |Q and shifts the coordinates of the particles + 1, ...,N + M by m 12. The result 



is ( |3.4|) . Integrability means that the system admits a Lax representation and has N + M 



integrals in involution. Moreover, since as t — > ±cxd these integrals reduce asymptotically 



to symmetric polynomials in particles' momenta, one can use the standard argument |^ to 
show that the S'-matrix is factorized. A small modification arises due to the presence of the 
boundary; namely, one can consider the particles' collisions both very far from the boundary 
where the problem is reduced to the bulk one, and near the boundary where the colliding 
particles have enough time to reflect and go to x = +oo. In the first case factorization gives 
the nonrelativistic Yang-Baxter equation for the bulk S'-matrix , while in the second case 
we get exactly the boundary Yang-Baxter equation of [0], which in the Dirichlet case allows 
to express the boundary S'-matrix of the anti-kink through that of the kink ( p.lOj ). The 
unitarity requires the latter to be a pure phase, but otherwise leaves it undetermined. So 
in both theories the S-matrix is factorized and fully determined by the bulk two-particle 
S-matrix and the boundary S-matrix. Thus to establish the equivalence it is sufficient to 
show that these S'-matrices coincide when the nonrelativistic limit is taken. 

Note that the translational invariance of the Hamiltonian ( p.4|) is broken not only by the 
boundary potentials, but also by the interaction of particles with their mirror images. This 
is very natural from the point of view of the underlying sine-Gordon theory, since, as was 
shown above, that the one soliton problem on a half-line is equivalent to the three-soliton bulk 
problem, with one of the particles staying at x = 0, and the other two being "generalized 
mirror images" of each other. The analogy becomes exact if we take y^o = 0. Then the 
ordinary method of images works, and it is obvious that the system of + M solitons on 
a half-line is equivalent to the system of 2(A^ -|- M) solitons on a line with symmetric initial 
conditions. Hence the corresponding nonrelativistic Hamiltonian can be obtained from the 
known |^, |36| nonrelativistic bulk Hamiltonian. One can easily see that the result is just 



with /i = z/ , /i(yU — 1) = p(p — l)/4. 
We will show below that yU and v are related to the parameter ^ of the sine-Gordon model 
ni) as follows: 

— ^ = ^- (3-11) 
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3.5 Taking the non-relativist ic limit 



The nonrelativistic limit of ( |3.1j ) corresponds to the values ^ 1. Simultaneously we must 
take the limit P <^ 1, so that Mg = ^ mo (otherwise the S-matrix becomes 1 and we 
do not get anything interesting.) In general, these two limits are to be taken without any 
further assumptions on the relative magnitude of 6 with respect to /3. However, it is worth 
to note that the region -C 6* ^ 1 corresponds to the quasiclassical limit ^ 1, where 

k = Ms9 (3.12) 

According to ( p. 7] ) ^ = which in general is not necessarily large. In what follows 
we assume that ^ is positive and P^po ~ 1, so that ^ scales as Then P+ has poles 

corresponding to the boundary bound states, and the energy of the bound states lying 
close to the edge of the continuous spectrum in the theory ( p^.l] ) becomes En — Mg — 

^^siE — 1 — 2nj . P„ does not have poles in the physical region. One can easily see then 
that (p.8|) , (|3.9|) describe correctly the spectrum of the boundary bound states provided that 
equations (|3J|),( pn| ) are fulfilled. To complete the identification of the boundary 5- matrices 
we have to compare the phase shifts. Since ^ scales as 1//?^, by virtue of (|3.11|) the expression 
( |3.9| ) can be rewritten as SNuik, v — \i)f{Q\ where 



^ Vaoy' ^ V 2 2ao 

■p I ik \ -p I II— f _j_ ik \ 

M ao) \ 2 ^ 



2aoJ 



is meromorphic and contains the poles located in the arbitrarily small neighbourhood of ^ = 
as /3 — s> 0, whereas the factor f{0) can be expanded into the power series f{6) = l + O'fi'' 
with all the coefficients and a radius of convergence ~ 1 as /5 ^ 0. In the same limit P± can 
be factorized analogously: P± = Sj\fji{k,±2^/TT)f±{9) with f± admitting expansions of the 
form f±{0) = 1 + I^i^i with all the coefficients and the radius of convergence ~ 1 in 
the limit (3^0. Therefore the boundary S'-matrices agree when ^ <^ 1, and Snr represents 
the nonrelativistic limit of the boundary S'-matrix of the sine-Gordon theory. One can check 
that the same statements are true for the bulk two-particle sine-Gordon S'-matrix of 
and the S'-matrix of the particles interacting via the potentials (|3.5| ), provided that (|3.6| ), 
( ^.71 ) are satisfied (this result was first established in in the quasiclassical approximation 
and later confirmed in note that our approach allows to give an exact sense to the 
statement that the nonrelativistic limit of the bulk sine-Gordon theory is the hyperbolic- 
type Calogero-Moser model.) Thus the equivalence of (|3.4|) and the nonrelativistic limit of 



(2.1) is estabhshed. 



3.6 Remarks 

It is instructive to compare also the combined nonrelativistic/quasiclassical limit of the S- 
matrix of ( p.l|) with the Poschl- Teller S-matrix in the regime (3"^ 9 1 (i.e. k ^ mo) 
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without appealing to the exact formula ( |3.1| ), similar to how it was first done in for 
the bulk soliton scattering. This means that one should first take the limit (3 —>■ and 
then, using (|2.34|) , pass to the limit ^ 0. Expanding the integrals in (|2.34|) and using the 



asymptotic formulas for the F-functions in ( |3.9| ) we get for ( |2.34| ) and ( p.9| ) in the principal 
order the following result: 

P+{k) = e "'o "0 , p_[k) = e'^o "^o , (3.14) 

once again confirming the equivalence of the two theories. Note that it is impossible to 
determine /i and v separately, since in our limit the boundary S'-matrices in both theories 
depend only on the difference n — v. 

The identification of the nonrelativistic limit in the case of the most general integrable 
boundary condition (|2.3| ) requires a nontrivial generalization of the Calogero-Moser Hamil- 
tonians. Indeed, if in (|2.1|) M < oo, then the boundary S'-matrix does not conserve the 
topological charge, and we are not aware of any integrable nonrelativistic model which al- 
lows such a process. 
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Chapter 4 

Boundary bound states and boundary 
bootstrap 

In the previous chapters we discussed classical, semi-classical and non-relativistic limits of 
the boundary sine-Gordon model. Now we address the exact quantum field theory solution 
of this model, which we obtain by means of the Bethe ansatz technique. Among other things, 
this solution allows to re-derive the boundary refiection matrices of ( |3.1| )-( P^ ), and to 
relate them to the physical parameters in the Hamiltonian [^]. The present chapter includes 
a complete study of boundary bound states and related boundary S-matrices for the sine- 
Gordon model with Dirichlet boundary condition. Our analysis is based on the solution of 
the boundary bootstrap equations, representing the integrability constraints, together with 
the explicit Bethe ansatz solution of the inhomogeneous XXZ model in a boundary magnetic 
field - a lattice regularization of the boundary sine-Gordon model. We identify boundary 
bound states with new boundary strings in the Bethe ansatz. This chapter is based on [Q. 



4.1 Introduction 



In the seminal work |I[ it appeared clearly that the boundary sine-Gordon model presents 
an extremely rich structure of boundary bound states, which was further explored in |^ . 
Our first purpose here is to study this structure thoroughly in the particular case of Dirichlet 
boundary conditions, that is the model 



SG 



1 

2 JO 



rrin 



(dtifY -{d,ipy + ^ cos (3ip 



dx 



(4.1) 



with a fixed value of the field at the boundary: ip{x = 0,t) = ipQ. 

Also, the consideration of boundary problems poses interesting challenges from the point 
of view of lattice models, here lattice regularizations of ( [4.11) . In and also in [52| it was 
shown in particular how to derive the S-matrices of [|l| from the Bethe ansatz. Our second 
purpose is to complete these studies by investigating which new types of strings correspond 
to boundary bound states, and by deriving as well the set of S-matrices necessary to close 
the bootstrap. Observe that lattice regularizations are useful to define what one means by 
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creating a bound state at the boundary. Indeed, some bound states (showing up as the pols 
of S-matrices) have no straightforward interpretation, and although they are easy to study 
formally using the Yang Baxter equation and the bootstrap, their meaning in the field theory 
is unclear. 

In section 4.2 we consider the bootstrap problem directly in the continuum theory. We 
identify boundary bound states and we compute the related boundary S matrices. In sec- 
tion 4.3 we write the Bethe ansatz equations for the inhomogeneous six-vertex model with 
boundary magnetic field, which is believed |3^ to be a regularization of ( |4.1| ). We show 
that these equations are also the bare equations for the Thirring model with f/(l)-preserving 
boundary interaction, which is the fermionized version of ( [4.1[) . In section 4.4 we discuss 
in details new solutions ("boundary strings") to the Bethe ansatz equations made possible 
by the appearance of boundary terms. In section 4.5 we study the physical properties of 
the model, in particular the masses and S-matrices corresponding to these boundary strings, 
and we partially complete the identification with the bootstrap results of section 4.2. Several 
remarks, in particular formula for the boundary energy of the boundary sine-Gordon model, 
are collected in sections 4.6-4.7. 



4.2 Boundary bootstrap 

4.2.1 Solving the boundary bootstrap equations 

The S-matrices for the scattering of a soliton (P^) and an anti-soliton (P^) on the ground 
state |0)b of the sine-Gordon model with Dirichlet boundary conditions ([4.1|) read according 
to i: 

P±(^) = cos{^ ± Xu)Ro{u)Ri{u,0, (4.2) 

where 6' = m is the rapidity, ^ = Atclpq/P and A = Svr//?^ — 1. The explicit form of Ro,Ri 
is given in ( p.2| )-( p73D . Since the theory is invariant under the simultaneous transformations 
^ — > and soliton— i>anti-soliton, we choose hereafter ^ to be a generic number in the 
interval < ^ < 47r^//5^ (see the discussion in chapter 2 about the value of the upper 
bound). 

The function Rq contains poles in the physical strip < Im^^ < 7t/2 located at m = 
n7r/2A,n = 1,2,... < A. These poles come from the corresponding breather pole in the 
soliton-antisoliton bulk scattering, and should not be interpreted as boundary bound states 

i- 

When ^ > tt/2, the function P^{0) has additional poles in the physical strip, located at 
u = Vn with 

f 2n + l TT , , , 

{n = 0,1,2,.. .) corresponding to a first set of boundary bound states which we denote by 
\Pn), with masses 

/^^ 2n + l \ ^. 
mn = m cos Vn = rn cos I — ^ — vr I , (4.4) 
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Ialpha> lalpha> 



Figure 4.1: A diagrammatic representation of the bootstrap equations (^4.5|) . 

where m is the sohton mass. These bound states are easy to interpret [Q, For < ipo < 
tt/P the ground state of the theory is characterized by the asymptotic behaviour — > as 
X oo, but other states, whose energy differs from the ground state by a boundary term 
only, can be obtained with ip ^ { a. multiple of 27r//5} as x — oo. Since the Pn appear as 
bound states for soliton scattering, they all have the same topological charge as the soliton, 
which we take equal to unity by convention, so they are all associated with the same classical 
solution, a soliton sitting next to the boundary and performing a motion periodic in time 
("breathing"), with ip{x = 0) = ipo and ip ^tt/P as x ^ oo (see section 2.2.4). 

To deduce the scattering matrices on the boundary bound states we use the boundary 
bootstrap equations as given in |^. We assume that these S-matrices are diagonal, which 
is true if all the boundary bound states have different energies. In this case the bootstrap 
equations read: 

RW) = E RMsf,{e + ivt)st{e - <J. (4.5) 

These equations allow us to find the scattering matrix of any particle h on the boundary 
bound state /5 provided that the latter appears as a virtual state in the scattering of the 
particle a on the boundary state a. 

The masses of the corresponding boundary states are related through 

mp = ma + ma cos (4.6) 

where iv^^ denotes the position of the pole, corresponding to the bound state /3. 

Let Pn stand for the n-th boundary bound state corresponding to the pole Vn in P"*" (|4.3| ). 
Then ( |4.5| ) gives: 

pue) = p+{e)a{e-iv^)a{e + ivn), (4.7) 
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Pp^O) = KO - ^^nW + tVn)P-ie) + - tvMO + ^Vn)P^ie), (4.8) 

where the well known bulk S-matrix elements a{6) = S^^ = SZZ (kink-kink scattering), 
b{6) = S^Z = SZ^ (kink-anti-kink transmission) and c{6) = S^l = Sl^ (kink-anti-kink 
reflection) can be found in [^j. 

It is easy to check that the matrix elements (|4.7] )- (|4.8| ) satisfy general requirements for 



the boundary S-matrices, such as boundary unitarity and boundary crossing-symmetry con- 
ditions IQ, e.g. 

Finally we obtain from (f4.7|) - (^4.8|) by direct calculation: 

cos(^ — Att — i\9) ^_ 
cos(^ — Att -|- iX9) 



Hence the boundary Yang Baxter equation is satisfied since the ratio of the above two 
amplitudes has a form similar to ( [4.2|) with ^ — ^ ^ — Avr, ^ being a free parameter. 

The analytic structure of P^^O) is as follows. The function P^^O) has simple poles in 
the physical strip located at n = | + ^^^^ tt , = 0, 1,2..., and at u = t>„. It has double 
poles at M = ivn + ^ ~ 1, 2, ...n. The function Pp^O) possesses in the physical strip the 
same singularities as Pp^O) plus the set of simple poles at u = wj^ with 

Interpreting these poles in terms of boundary bound states requires some care. First, due 
to the relation ( |4.5| ), one sees that if /? appears as a boundary bound state for scattering of 
a on a, then the poles of the amplitude for scattering of 6 on a are also in general poles of 
the amplitude for scattering of b on (3. It seems unlikely that these poles correspond to new 
bound states, although in our case they would have a natural physical meaning, for example 
one could try to associate them with classical solutions where y?— >-47r//3asx— s>oo. Indeed 
there are strong constraints coming from statistics that we should not forget. For instance 
at the free fermion point = 47r, there is a bound state but although has again a 
pole at Pi, the state of mass 2m/3j is not allowed from Pauli exclusion principle, as can easily 
be checked on the direct solution of the model (see below section 4.3.3). Therefore we take 
the point of view that the poles already present in the scattering on an "empty boundary" 
are "redundant" . The only poles we interpret as new boundary bound states are (|4.1CI|) (the 
additional poles in are related to them by crossing). We denote these boundary bound 
states |5n,7v), and their masses, according to (|4.6| ) and (^^), are given by 



2n + l \ , 2iV-l 

fTT'n,N = m[COSVn + COSWn) = mCOS [ — TT — m COS -- H TT 



A 2A J \X 2A 

mUn sin I ^ + — TT 1 , (4.11) 
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where = 2m sin f is the mass of the p-th breather, p = 1, 2, . . . < A. 

To understand the physical meaning of these new boundary bound states it is helpful 
to consider the semi-classical limit A — oo of the sine-Gordon model. As discussed above, 
the boundary bound states /?„, corresponding to (^I3|), are associated with solutions where a 
soliton is sitting next to the boundary and "breathing" . An incoming anti-soliton can couple 
to this soliton, and together they form a breather sitting next to the boundary and performing 
again some (rather complicated) motion periodic in timeQ. The WKB quantization of this 
solution shoud lead to \6n,N)- The topological charge of the states |5„,Ar) is equal to in 
our units, or, equivalently, to the charge of a free breather in the theory (|4.1|) . 

One can in principle continue to solve the bootstrap equations recursively. For 

example, for the scattering of solitons or antisolitons on the boundary bound states \Sn,N) 
( [4.10|) one obtains the following S-matrices: 

Ps-je) = PpS(^)a{e - iwM)a{e + tw^), (4.12) 

P+ - cos(e-^Ag) 

Pf^ ^ has only one simple pole in the physical strip at u = w^, while Pg^ ^ has also simple 
poles at M = ffc, k = n + l,n + 2, [| — |]. According to the discussion below ( [4.10|) , we 
do not consider these poles as associated with new boundary bound states. Therefore, the 
boundary bootstrap is closed for solitons and antisolitons in the sense that further recursion 
will not generate new boundary bound states. 

So far we have obtained two sets of boundary bound states (|4.4| ) and (|4.11| ) by considering 



all the poles in the physical strip of amplitudes for scattering a soliton and anti-soliton on 
a boundary with or without a boundary bound state. Of course we should also consider 
the scattering of breathers off the boundary. The scattering of breathers on the "empty" 



boundary was studied in |41|, and we refer the reader to this work for the explicit boundary 



S-matrices. By interpreting the poles of the amplitudes in |^ as boundary bound states, 
we find a spectrum of masses that look like ( [4.11|) but with a slightly different range of 
parameters. Considering then scattering of breathers off a boundary with a bound state does 



not give rise to any new poles beside ( [4.3| ) and ( [4.10| ), with in the latter case an extended 



range of values of (for simplicity we do not give the relevant boundary S-matrices here). 
Therefore the complete boundary bootstrap is closed in principle. 



4.2.2 Integral representations of various S-matrices 



For comparison with results obtained from regularizations of the sine-Gordon model it is 
useful to write integral representations of the boundary ^-matrices ( [4. 21) , (|4.7| ) and 
using the well-known formula 



logr(z) 



°° dx 



X 



Z-1 + 



Rez > 0. 



(4.14) 



^To compute this solution explicitly requires using a bulk five-soliton configuration p^ , an expression 
which is very cumbersome. 
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Suppose first that l<2^/7r<A + l and denote 

TT 



(4.15) 



where the square brackets mean the integer part of the number. For such values of ^ there 
are n^: + l poles (|4.3| ) in the physical strip, i.e. the spectrum of excitations contains boundary 
bound states. Correspondingly, there is a finite number of F- functions in (|4.2|), (^.7|), ( [4.8| ) 
whose arguments have negative real part so that formula ([4.14|) is not applicable. Treating 
such F-functions separately, we obtain the following results: 

^{6) 1 _ 2A /•+°° (2Xex\ 



d . 

— fog 



(4.16) 



X 



sinh(2^/7r - 2n^ - 2)x ^ sinh(A - 2^/it)x' 



sinh X 



2 sinh X cosh Xx 



.d 



2A [+^ , /2A^a;\ 
— / dx cos 

TT J-oo \ TT J 



(4.17) 



sinh( A — 2^/ tt)x — 2 cosh x sinh( A + 1 + 2n — 2^/ 7r)x 



2A 



+~ , /2A^x\ 

dx cos 

71 J-oo \ ir J 



2 sinh X cosh Xx 

sinh(2n* + 2 - 2^/Tr)x 



sinh a; 



(4.18) 



+ 



sinh(A — 2^/'k)x — 2cosha;sinh(A + 1 + 2?t, — 2^/tt)x 
2 sinh X cosh Ax 



In the derivation of analogous representation for P there are no subtleties because the 
"dangerous" F-functions cancel. We get 



.d 



Ro{0) 



2X 

TT 



+°° , /2Aex\ sinh(A - 2^/7r)x 

dxcos — - — . 

-oo \ 71 J 2 sinh X cosh Ax 



(4.19) 



In the region < 2,^/7r < 1 where there are no poles and no boundary bound states in the 
spectrum, formula ( [4.19|) is valid, too. The expression for P"^ can be obtained from ( [4. 161) 
by setting formally ub = + 1 = 0, which gives 



.d 



Ro{0) 



2X 



+°° /2Aex\ sinh(A + 2^/7r)x 

UiOC COS I I ~ ~ ^ • 

7r J-oo \ vr / 2 sinh X cosh Ax 



(4.20) 



Note that if 2^/7r > 1, the integral in ( |4.2CI|) diverges. Finally, we complete this list by the 
following two expressions: 



X 



.d 

sinh(f 



PtM 
Ro{0) 

' 2?2* — 2)x 



.d 



Pli(^) 



2X 

TT 



sinhx 



RoiO) 

2 cosh X sinh (2S + 2N 
2 sinh X cosh Ax 



, /2A^x\ 
dx cos 

-oo \ 71 J 



X - Dx 



(4.21) 



For the integral representation of Rq see |]39[ . 
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4.3 Exact solution of the regularized boundary sine- 
Gordon model 



4.3.1 The XXZ chain with boundary magnetic field 

The XXZ model in a boundary magnetic field 
ri-i 



n 



7 



27r sin 7 



E + ^MVi + ^i^^^Ui - 1)) + hicrt - 1) + h\al - 1) 



Li=l 



(4.22) 



was discussed in |^4[ and in |^5[, where its eigenstates were constructed using the Bethe 
ansatz. As usual, these eigenstates T-C\n) = E\n) are linear combinations of the states with 
n down spins, located at Xi, on the chain: 

1^) = ^/^"''(^l' ■■■,Xn)\Xl, ...,Xn). 



Consider for simplicity the case n = 1. The wave-function f^^\x) reads 



(4.23) 



sinh |(i7 



a 



sinh i (27 



a] 



sin 7 sinh ^ (a + ^7-/^') 



sinh 1(^7 — a) sin |(7 + 7-ff') 
where we defined the new variables as in [45]: A = — C0S7, k = f{a,'y), 

h — i sin 7 



(a 



-a). 



'yH = /(z7, —i \n(h — A)) = —7 — z In 



h + i sin 7 



(and similarly for H'), and 



/(a, 6) 



-i In 



sinh ^{ib 



sinh i(i6 + a) 



(4.24) 



(4.25) 



When varies from to -|-cxd, 'jH increases monotonically from — tt — 7 to —7 according to 
( [4. 241) if we take the main branch of the logarithm. 

Denote h^h = 1 — COS7. This "threshold" value of h corresponds to 'jH = —tt; its meaning 
will become clear below. When h varies from —00 to 0, 'jH increases monotonically from 
—7 to TT — 7. For the purposes of the present work we confine our attention to the region 
h,h' > and choose 7 G (0, |). Other regions in the parameter space can be obtained using 
the discrete symmetries of the Hamiltonian ( [4.22| ): — >■ —a^ on each site or on the odd 
sites only. The parameter k in ([4.23|) is not arbitrary, but satisfies the Bethe equation |45 



e''' + h - A){e''' + h' - A) 
(e-*^ + h- A)(e-*'= + h' - A) 



or 



sinh i(a 



sinh |(a; + ij) 



2L 



sinh 2 (a 



i'')H) sinh ■ 



sinh |(a + i'jH) sinh |(a + i'jH') 



(4.26) 



(4.27) 
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Note that the wave-function ( |4.23| ) depends on H imphcitely through the solution of the 
Bethe equation ( ^.271 ) a{H,H'). Besides, one can multiply the amplitude ( [4.23| ) by any 
overall scalar factor depending on a, L, H and H'. The Bethe equations in the sector of 
arbitrary n > 1 can be found in [HH . 



4.3.2 The Bethe equations for the inhomogeneous 
XXZ chain 

The real object of interest for us is actually the inhomogeneous six-vertex model with bound- 
ary magnetic field on an open strip. The inhomogeneous six-vertex model is obtained by 
giving an alternating imaginary part ±iA to the spectral parameter on alternating vertices 
of the six- vertex model |4^. It was argued in generalizing known results for the 
that this model on an open strip provides in the scaling limit A, L — * oo, 
a lattice regularization of (|4.1|), with (3"^ = 87 and a value of (fo at the 



periodic case |^ 



lattice spacing 

boundary related to the magnetic field. The reader can find more details on the model in the 
references; it is actually closely related to the XXZ chain we discussed above. In particular, 
the wave function can be expressed in terms of the roots aj of the Bethe equations 0: 



sinh ^{aj -f A — ^7) sinh — A — ^7) ^ sinh ^{aj — i'jH) sinh ^{aj — i'-fH' 
sinh |(a;j + A + ^7) sinh — A + 27) sinh ^{aj + i'jH) sinh ^{aj + i'jH' 



= n 



2 

sinh l{aj — 



■ i'jH) sinh |(a;j — i'-fH' 



2^7) sinh + am — 2^7) 



sinh ^{aj — a 



■m 



m 

227) sinh \ {oij + a 



m 



+ 2^7) ■ 



(4.28) 



By construction of the Bethe-ansatz wave function, Rettj > 0. Note that the solutions of 
( [4.28|) aj = 0, in should be excluded because the wave function vanishes identically in this 
case. The analysis of solutions of eq. ( [4.28| ) is very similar to the case of the XXZ chain 
in a boundary magnetic field. We consider the regime < 7 < 7r/2, which falls into the 
attractive regime < < 47r in the sine-Gordon model ( [4.1D . We set hereafter 7 = vr/t and 
for technical simplicity restrict t to be positive integer. In the limit L — 00 this constraint 
implies that in the bulk only the strings of length from 1 to t — 1 are allowed, along with the 
anti-strings (see chapter 1). 

Taking the logarithm of eq. (f4.28 ), one obtains: 



L [f{aj + A, 7) + /(«, - A, 7)] + fia„ iH) + /(a,-, ^H') 
= 27ilj + J2 [/("i - ^rn, 27) + f{aj + am, 27)] , 



(4.29) 



where Ij is an integer. We also recall the formula for the eigen-energy associated with the 
roots aj 0, 



E 



TT — 7 



TT 



E[/K + A,7) + /'(«,-A,7)]. 



(4.30) 
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4.3.3 Thirring model with boundary 

Like the bulk sine-Gordon model is a bosonized version of the bulk massive Thirring model, 
one can expect that the boundary sine-Gordon model is a bosonized version of the Thirring 
model with certain boundary conditions. The quickest way to identify this boundary Thirring 
model is to use the Bethe ansatz equations ( |4.28| ). Write the most general f/(l)-invariant 
boundary interaction 



+ E «u-^M-(0) + Y: <,^t^jiL). (4.31) 







The entries of the 2x2 matrices A = {aij},A' = {a'ij} can be determined up to one 
arbitrary parameter by the hermicity of and the consistency of the boundary conditions 
(detA = 0). For the left boundary, the matrix A looks like 

and the boundary condition reads V'i(O) = — e*'^'?/'2(0) (similarly for the right boundary). 

To find the eigenvectors of the Hamiltonian ( [4. 311 ), Ht"^ = E^^, one can use the same 
wave-functions as for the bulk Thirring model ( |1 . 7|) , and modify them by analogy with the 
example of XXZ chain in a boundary magnetic field This way one gets the equations 
for the set of rapidities aj : 



^2imoL sinh aj 



cosh |(aj + icj)) cosh ^{aj + 
cosh |(aj — icj)) cosh |(aj — icj)') 



X 



-j-j- sinh l{aj - am- 2^7) sinh |(q;j + am - 2^7) 
sinh l{aj - am + 227) sinh l{aj + am + 2^7) ' 

where 7 is related to go in the usual way [l^. These equations look quite similar to ( [4.28| ). 
The mapping can be made complete by taking in ( 4.28| ) the limit A ^ 00 with the identifi- 
cation niQ = 4e~^ sin 7. 

The derivation of these equations is rather cumbersome, therefore to illustrate the pro- 
cedure we comment on the simplest case of one-particle sector, which is nevertheless suf- 
ficient to obtain the form of the boundary terms in (|4.33|) . We make an ansatz \1/ = 
lo d,yx^{y)'ipx{y)\0), where A is the spinor index, xiv) is the wave-function and |0) is the 
unphysical vacuum annihilated by ipx. 

The equation if^^ = E'if reduces to 

d 

-^(^3^X + moaix + AxS{x) + A'xS{x - L) = Ex- (4.34) 
where cxj are the Pauli matrices. We look for the solution of (|4.34|) in the form 

) = a{a) ( "^T/f ) e^™°^^^°'^" - a(-a) ( f'/, ) e^^o--'^" (4.35) 
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Substituting it into (|4.34|) we get E = mo cosh a and, besides, two boundary conditions to 



be solved. The first one, at x = 0, determines the form of the factor a{a) = cosh |(a; — 
while the second one aX x = L gives rise to the Bethe equation 

g2*moLsinha _ COsh |(q; + i0) COsh |(a + i^Q 

cosh \ {oi — i(f)) cosh ^{a — icf)') ' 

which determines a. Comparing the Bethe equation ( [4.28| ) with ( |4.33| ) and using the relation 
between ^ and H obtained below in section 4.5 we find the relation between the boundary 
parameters and ip^ in the Hamiltonians (|4.31|) and (|4.1| ) respectively: 



Thus, the integrable boundary condition for the f/(l)-invariant boundary Thirring model 
equivalent to ( |4.1|) reads: 

^2(0) = -e*/3V8-i/3^o^^(o). (4.36) 



It would be interesting to obtain the result ( [4.36|) directly from the Hamiltonian ( [4.1|) 
using an extension of the Coleman- Mandelst am bosonization technique to the case with 
boundary. However, to our knowledge such an extension has not been developed yet, except 



for the free point |48 . 



4.4 Solutions of the Bethe ansatz equations with bound- 
ary terms 

As is well known in the case of the bulk Thirring model or equivalently the periodic XXZ 
chain, the bound states are associated with various types of solutions of the Bethe ansatz 
equations involving in general complex roots (see chapter 1). By analogy, we expect the 
boundary bound states to correspond to new solutions made possible by the boundary terms. 

Consider first the example of the XXZ chain as given in section 4.3.1. Since our goal is 
to study purely boundary effects, we will look for the solutions of the Bethe equations that 
give rise to a wave-function localized at x = or x = L and exponentially decreasing away 
from the boundary. 

The states described by such wave-functions will be referred to as the "boundary bound 
states" below. For this, one should have a purely imaginary in ( [4.23| ). We consider here 
the limit of L large, when the left and the right boundaries can be treated independently 
and the overlap of the corresponding wave-functions is negligibly small (for the physical 
applications it is necessary to take the scaling limit anyway). In the limit L ^ cxd it is easy 
to check that there are two such solutions to ( |4.27|) : a = ioq = —i'-fH + ie{L, H, H') and 



a = ia'o = —i'~)H' + ie'{L, H, H'), where e ~ exp(— 2kL) and we defined k > as 

sin|(-7if-7) 



sin \{-iH + 7) 



51 



(similar relations are assumed for e',K'). Solution gives a wave-function ( [4.23 ) localized 



aX X = L: f^^^{x) ~ e '^''^^ Solution ao gives a wave-function localized at x = 0, 
f^^\x) ~ e~'^^, provided we renormalize the wave-function ([4.23|) : 



sinh 2 ('-^ ~ i'-fH) sinh 2 ('^ "I" ^7-^) 



1/2 

(4.37) 



In the special case H = H' there is only one proper solution a = iao = —i'yH + ie{L,H) 
with e ~ exp(— kL). The wave-function ( [4. 231 ) behaves as the superposition of the "left" and 
the "right" boundary bound states, f^^^ ~ ^^-kx _^ ^-K{L-x)y jsjote that the boundary bound 
state appears in the above example only when the boundary magnetic field is large enough:^ 
namely, h > hth- This follows from the fact that a should be such that < aQ < it. 

Now, consider the equations for the inhomogeneous model ( 4.28|) . The basic boundary 



1-string solution to ( [4.28| ) is a = iao = —i^H + ie, provided that < ao < ^r. This solution 
is possible due to an argument very similar to the one used in the bulk: as L — ^ 00, the two 
first terms of ( [4.28| ) decrease exponentially fast, while the third increases exponentially fast, 
and e ~ exp(— 2kL) with 

_^ _ sinh^ A + sin^ ^ 
^ - sinh2A + sin2£^^- ^^'^^^ 

Recall that for the bulk problem when there is no boundary term, the right hand side of 
( [4.28|) would have to decrease exponentially, forcing the existence of a "partner" root at 
a — 2i'y. 

One can construct similarly boundary n-strings which consist of the points iao, 
iao+2i'y, iao + 2i{n — l)'y (see figure 4.2). By convention n = means there is no boundary 
string, that is all complex solutions are in the usual bulk strings. The possible values of n 
are restricted by the fact that the upper point of the complex should be below in: max(?T,) = 



-aol 



L 2^ J + 1) where the square bracket denotes the integer part. To show that the boundary n- 
string is indeed a solution to ( [4.28| ), we introduce infinitesimal corrections Ei to the positions 
of the points of complex [Q. Taking the modulus of both sides of ( [4.28 ) with aj = iao + 2ik'y 



and multiplying equations for = 0, n — 1 we obtain exp{— 2L(ki + ^2 + ••• + K,n)} ~ ^1, 
where ei denotes the correction to the point iao. The behavior of the remaining Sk follows 
from 61 by recursion. For example, for the 2-string 62 is given by l^i — £2! ~ exp(— 2LK2). 
Note that associated with each boundary n-string there is also the solution to ( [4. 281) 



obtained by complex conjugation of all a's. The existence of such a "mirror image" is the 



2 



More generally, the criterion of existence of boundary bound state solutions allows us to determine 



threshold fields for any anisotropy A. For this, let us examine (4.26). The parameter k is defined modulo 
27r, therefore we restrict it to lie withing k G (0, 27r). Two possibilities k = ia and k — t: + ia, where a > 0, 
lead to two different threshold fields, determined by the fact that the denominator in ( 4 .26| ) should vanish: 



and the regions where boundary bound states could in principle exist are h < A — I, h > A + 1 (one has 
to be carefuU here and check that these solutions of BE indeed correspond to the stable states of the model) . 
When A > 1, there are two different threshold fields, in agreement with the results of Jimbo et al In 

the region of interest, |A| < f, there is only one threshold field h[j^ . 
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Ima 

71 Dirac sea 



,,iaQ + 6^7 
,,io;o + 2^7 

Rea 



Figure 4.2: The first type of boundary string. In tlie ground state tlie boundary string of 
maximum allowed length is filled. 



consequence of the symmetry of equations ( |4.28| ) and it is of no importance to physics. In the 



bulk case, it is easy to show that all solutions are invariant under complex conjugation, 
but this result does not hold here. In fact, a solution which has both the boundary n-string 
and its mirror image would lead to a vanishing wave-function. 

Additional boundary strings can be obtained by adding the roots iag below iao so that 
ias = iao ~ 2is7, with s = 1, 2, (see figure 4.3). Together with the boundary n-string 
above ao, they form the complex which we call boundary [n, A^)-string. To analyze the 
existence of such complexes as the solutions of ( |4.2(j| ) we introduce as before the infinitesimal 



corrections Es to the roots ag, where now s = n,n — 1, 1, — 1, — 2, — A^. Then, the 
equations ( ^.281 ) with aj = iug tell us that the range of A^ should be 

^ < AT < (4.39) 
27 27 ^ ^ 

In other words, the inequality ( [4.39|) states that the lowest root of the boundary string should 
be below the axis Ima = and above the axis Ima = — ivr. Another constraint follows if 
we multiply the equations (|4.28|) for all the roots of boundary (n, A^)-string. This gives 



exp(— 2LX] f^s) = ^i- So, one should have X) > 0. The latter sum can be easily evaluated 
if one uses the expression ([4.38|) simplified in the limit A 00: k = 4e~^ sin7 sin oq. The 
constraint obtained in such a way forces the number of roots above Ima = axis in the 
boundary string to be greater than the number of roots below Ima = 0. 

We have not been able to find any reasonable additional solution to the Bethe ansatz 
equations. The two sets of boundary strings we have encountered appear to be in one to one 
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Figure 4.3: The second type of boundary string. 
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correspondence with the boundary bound states identified in section 4.2 using the bootstrap 
approach. To clarify this identification we now compute related masses and S-matrices. 



4.5 S-matrices and bound state properties from the 
exact solution 

4.5.1 Bare and physical Bethe ansatz equations 

The hare Bethe equations follow from taking the derivative of ( [4.29|) . Defining 2L{pk+ p'l)da 
to be the number of roots in the interval da, one obtains coupled integral equations for the 
densities of strings pi, ...,pt-i and anti-strings pa- 

1 1 

2vr(p, + p^) = -ip'^ + /'*p, + |]/>p, + ^(«„ + ^/f)+^75^)) (4.40) 
where * denotes convolution: 

f*g{a)= / dpf{a-P)g{P). 



These densities are originally defined for a > 0, but the equations allow us to define pk{—a) = 
Pk{(y) in order to rewrite the integrals to go from — oo to oo. If we totally neglect the boundary 
terms (terms ~ L~^) in ( [4.40 ), we will end up with the same equations as for the periodic 



inhomogeneous six- vertex model |4^. The various kernels and sources in ( [4.40|) are defined 
as follows: 

Pa{a) = /(^TT + a + A, 7) + f{in + a - A, 7), 
Pfc(a) = J2 fi'^i + A, 7) + fiai - A, 7), 

where the sum in the last expression is taken over the rapidities of the bulk k-string roots 
centered on a. 

The kernels fki are the phase shifts of bulk k-string on bulk /-string obtained by summing 
( [4.25|) over the rapidities of string roots. The boundary terms are: 



Ua{a) = -2/'(2a, 27) - /'(« + ztt, ^H) - f'{a + tir, 'jH') - 27r5(«), 
Ukia) = E[2/'(2«., 27) + /'(«., iH) + f'ia,, jH')] - 27r<5(a), 
(the sum above is over the roots of bulk k-string centered on a), 

•''^ («) = E /(*^ + 27) + /(^vr + a + a„ 27), 
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and ftj denotes the rapidities of the roots in the boundary n-string. 



where a, denotes the roots in the boundary n-string, while aj denotes the roots in the bulk 
k-string centered on a. The parameters uj,uj' are equal to 1 or 0, depending on whether the 
boundary string is present or not. In our ferromagnetic case the ground state of the periodic 
inhomogeneous XXZ chain is filled with anti-strings. 

The physical Bethe equations are obtained |]5T| , |52| by eliminating the "non-physical" 
density pa from the right-hand side of ( [4.40|) . This is done simply by solving for pa in the 
last equation in (|4.4CI|) and substituting it into the others. The result is 



Mpa+Pa) 



1 1 f 



J ak 



7 * 2vrp^ 



E (fL + 

1=1 \ 



J akJ a 



al 



277-/7 

1 2np'^ f 



271 -f 
* Pl + ■^Un,n';k, 



f 

J al 



2 271 - f 271- f 



; * 27rp^ + ^ — ^ * 27rp, 



(4.41) 



2L 



where 



Un,n';a 
Un.n'ik 



Ua + UjfP + Uj'f:^, 



27r 



ci{R) 



271- f 

- ^-Jfn± - ^75? - fak * f/„,„';a/27r, 



(4.42) 

(4.43) 



and different products (ratios) of kernels are defined through their Fourier transforms. 



4.5.2 The mass spectrum of boundary bound states 

We assume at first that the ground state is built by filling up the Dirac sea with anti-strings, 
as in the case of the periodic XXZ chain. We will see below that this is not always true. 
The presence of the boundary strings in the Bethe equations deforms the distribution of 
roots and modifies the density of the Dirac sea pa by a term 6pa/2L of order L~^. With the 
boundary n-string, the Bethe equation for the density of the Dirac sea particles pa is 

1 1 /■+°° 1 

2P'a(«) - ^«a(«) = f'ia - f3)pa{mP " 2^P~a{a) + i^f». (4-44) 



where was defined above. Subtracting from ( |4.44| ) the equation for the density of the 
Dirac sea alone, 

1 1 /'+°° 

2p1(«) - ^^-(") = j ^ - P)paWP - 27rp,(a), (4.45) 



56 



one obtains the equation for 6pa- 

n+OO 







(4.46) 



where 



5 Pa = 2L{pa - Pa)- 

The solution to ( 4.46|) can be written in terms of the Fourier transform 



as follows: 



5paik) = J dae''"^6paia) 



2vr - f'{k) 

For the boundary n-string ia^ + 2i'ys, s = 0, 1, n — 1, we obtain 

4 cosh7/c sinh wyk cosh(ao + 'jn — 'y)k 



Spa{k) 



where we used 



sinh Tik 

2 cosh7/c sinhn7A; cosh(a;o + 'jn — 'y)k 
sinh 7A; cosh(7r — '~f)k ' 

sinh(7r — 27) 



(4.47) 

(4.48) 
(4.49) 



f\k) = 271 



sinh Tck 

Expressions (|4.48| ), ( |4.49|) are valid for the n-strings with n 
n-string with n 



1,2,..., 



t+H 



For the longest 



t+H 



+ 1 = + 1 the Fourier transforms /4, Spa differ from ( [4.48|) , ( [4.49| ): 
2 sinh(7r — 27)/;; cosh(Q;o + 27n=i, — n)k 



2ti 
2tx 



sinh Tik 

4cosh7A;sinhn*7A;cosh(Q;o + 7n* — '~f)k 
sinh irk ' 



(4.50) 



Spa{k) 



sinh(7r — 27)^; cosh(a;o + 27n^, — TT)k 

sinh 7/c cosh(7r — 7) A; 
2 cosh 7fc sinh 77,^,7/;; cosh(Q;o + 7?T,* — ■y)k 
sinh7fccosh(7r — '~f)k 



(4.51) 



The conserved U{1) charge in the boundary XXZ chain is the total projection of the spin on 
the z-axis. In the thermodynamic limit the charge of the boundary n-string with respect to 
the vacuum is determined by 



/•+00 r+00 ]^ r+00 

Qn = n+ 2Lpada- 2Lpada = n+- 5pada = n + -5pa{0). (4.52) 
JO Jo 2 J-00 2 
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Using ( [4.491 ), we obtain for the n-string Qn = 0, and for the longest boundary string Eq. 
( [4.51| ) yields = 7r/27. Similarly, the mass of the boundary strings in the thermodynamic 
limit according to ( [4.30|) is given by 



r+oo r+oo ][ r+oo 

mn = K+ / 2Lpahada - / 2Lpahada = K + - / hjpada, 

JO JO I J-oo 

where the expression for ha is 

? TT — 7 ^. , 2 sinh7fc COS Afc 

ha{k) = Pa = -2(7r - 7)- 



(4.53) 



TT 



sinh nk 



and the soliton mass |^| 

We obtain in the limit A ^ oo 



m = 2e 2(t-7) , 



rrir. 



m 



sin ^{2n-l-H)+ sin ^ {H + 1] 
2 a 2 a 



IT 



= m sm — (H 
2\ ^ 



V. 



(4.54) 
(4.55) 



Since the parameter H varies in the interval —A — 1 < if < — 1, the mass of the longest 
string m=K ( [4.55| ) is always negative, while the other boundary strings have positive masses 
( [4.54|) . This means that the vacuum we have been working with is an unstable one in the 
region —t<H<—l{h> hth)- To cure the situation we define a new correct ground state 
by attributing the longest boundary string to the Dirac sea. The boundary excitations are 
obtained by succeessive removing of particles from the top of the longest boundary string. 
The charge and mass of such excitations with respect to the correct ground state are given 
by 



Qn = - 



27' 



TT 

m cos —{X + 1 + H — 2n) , 
2a 



n 



0,1, 



(4.56) 



Note that the number of excitations ( |4.56| ) is equal to the number of roots in the longest 
boundary string, + 1. The charge of such boundary excitations is equal to the charge of 
the hole in the Dirac sea. We identify a hole with a sine-Gordon soliton, and the boundary 
excitations described above, with the boundary bound states |/5„) ( [4.4|) . Their masses and 
charge ( [4.5(j| ) and the counting coincide provided that 



t + H 



2i 

TT 



(4.57) 



and the lattice charge Q is properly normalized. This expression is in fact valid for all values 
of /i > 0. 

In the above discussion we considered the boundary bound states related to one of the 
boundaries (say, the left one). In principle, one should include into the ground state the 



longest boundary string ia'^ + 2^7/, / = 0, 1, 



t+H' 



corresponding to the right boundary 



as well. The energy of the excitations due to both boundary strings is a superposition 
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of energies of the form (|4.5ti|). When H = H', these two boundary strings overlap and the 



usual Bethe wave-function vanishes. However on physical grounds we do not expect anything 
special to happen when the boundaries are identical. So, in such a case one should use as a 
wave function a properly renormalized version of the limit H H' of the usual Bethe wave 
function. 

When the magnetic field varies, the above picture indicates a qualitative change in the 
structure of the ground state at values H = —t, — t + 2, — t + 4, .... At these values, the 
mass of the bound state with the highest mass approaches the soliton mass and it becomes 
unstable. As discussed in |^3[ and for the Ising case, this decay corresponds to large 



boundary fluctuations that propagate deeply into the bulk. 

The mass of the boundary {n, A^)-string with respect to the correct vacuum can be 
calculated analogously. The result is: 

e vr e 2n + l ^ ,e , ^ 
mn,N = mcos(- - — ) + mcos(- 2A~^'^ - mcos(- H — — vr), (4.58) 



where we used ( [4.57|) to express H in terms of ^. This result seems rather confusing, because 



the above mass does not correspond in general to one of the bound state masses found in 
the bootstrap apporach. It can be considered as a sum of such masses, hinting that the 
(n, N) string describes actually coexisting bound states, but the calculation of correspond- 
ing boundary S-matrix does not allow such an interpretation. We are forced (but see the 
conclusion) to consider that only the (n, A^)-strings with n = + 1 occur, that is the phys- 
ical excitations are built by adding roots to the ground state configuration below iao. The 
charge and energy of such excitations with respect to the correct vacuum is given by 

Q7V = — - — = 0, mAf = mcos(- - — ) -mcos(- + tt), (4.59) 



These coincide with the charge and mass of the boundary bound states |5„=o,Af) (|4.11|) . The 
range of ( |4.39| ) agrees with the range of corresponding parameter in ([4.10| ). 

4.5.3 Boundary S-matrices 

It remains to check that the boundary S-matrices obtained above by the bootstrap approach 
coincide with those of the lattice model. To extract the boundary S-matrices from the 
Bethe equations we will follow the discussion of [03. Briefly, the idea of the method is the 



following. The physical excitations of the lattice model in the limit A ^ oo can be thought 
of as relativistic quasi-particles with rapidities 9i. The integrability implies that the set {9i} 
is conserved. Moreover, if the scattering matrices are diagonal, each particle preserves its 
rapidity. Assuming that this is the case, the quantization of a gas of A/" quasi-particles on 
an interval of lenght L results in the integral equations for the set of allowed rapidities 



p 1 

27r(pb + ) = nib cosh 9 + ^Lpbc* Pc + 77766, (4-60) 

c=l 
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where subscript stands for the type of particle, and 



. d 



(4.61) 



dO 



dd 



Equations ([4.60|) should be compared with the physical BE ([4.41|) , which gives bulk and 
boundary S-matrices. We will confine our attention to the boundary S-matrices only, keeping 
track of those terms in ( f4.42| ), ( [4.43|) , ( [4.61| ), which depend on the boundary magnetic field 
(the field- independent terms contribute to Rq and their agreement has been shown in ||39(| ). 
The discussion for the left boundary is completely parallel to that of the right one. Also, 
it is sufficient to consider only h =soliton and h =anti-soliton in ( [4.60|) . We identify a hole 
in the anti-string distribution in ( [4. 411 ) with a soliton in ( [4. 601 ), and {t — l)-string with an 
anti-soliton. Below we give explicit expressions only for the kernels in ( [4.41| ) which we need 
for our analysis. The other expressions are listed in [^. 

Suppose first that h < hth {—t — I < H < —t). This corresponds to the case without 
boundary excitations in the spectrum, ^ < 7r/2. Choose u; = a;' = in ([4.41|) . Then 



sinh(27r + -fH)k 



+ 



sinh nk 

(we omitted the if-independent terms and if'-dependent ones) 



m = 2n 



sinh(27r + -fH)k 
2 sinh 7/c cosh(7r — '~f)k 



+ .... 



Using (|4.61| ) we compare this expression with ([4.20|) (recall that the rapidity a should be 
renormalized a ^ 9 = ta/2X) and find complete agreement under the identification ( [4.57| ). 
Similarly, one can use 



Ut-l 



-27T 



sinh(7r + '~fH)ksmh.{7T — '~f)k 



ML) 



-2ti- 



sinh vrA; sinh 7 A; 
sinh(2 + H)-ik 



+ ... 



+ . . . 



2 sinh7/c cosh('7r — 7) A; 

to compare Vt-x with ( ^19D and obt ain agreement as well. 

Next, suppose that h > hth {—t < H < —1). To obtain the boundary S-matrices for 
scattering on the ground state \0) b set u = u' = 1 and choose the boundary string to be the 
longest string, n = n^, -|- 1 in ( |4.41| ). Then, using ( [4. 50] ), fn,+i-t-i — ~f'* ^'^d 



u 



(L) 



sinh 'jHk 

27T . / , + . . . 



sinh nk 



u 



(L) 



^^ sinh(ff + 2[^])7fc ^ 
sinh ■jk 



(4.62) 
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we obtain 



tAL) 



sinh '-yHk 



27r 2 smh7fc cosh(7r — 7)fc 



+ 



sinh(7r — 2'~^)k cosh.{H + t — 2n^)'~^k 
sinh 7 A; cosh (vr — 7)/c 



2 cosh 7A; sinh n^'^k cosh(iJ — n^, + 1)7^ 



sinh 7^ cosh (tt — ^)k 



+ 



(L) 

n,+l;a 



27r 



sinh(i^ + 2[i^])7A; 



+ 



sinh 'jk 

which agrees with ( 4.16|) , ( [4. IS ) under the identification ( [4.57| ). Note that the last relation, 
which follows directly from ( [4.42| ), (|4.43| ) and (|4.62|) , is valid also for Un;a and Un-t-i with 
any n. In the same manner one can calculate the boundary S-matrices for scattering on 
the boundary n-strings and check that they indeed coincide with (|4.17| ), ( [4.18|) under the 
condition ( [4.57|) . For this, one needs to take u = uj' = lin ( |4.41| ) and use ( |4.48| ), = —fn- 

Finally one can compute also the boundary S-matrix for the scattering on the {n^, + 1, A^)- 
strings, again in agreement with the bootstrap results. 



4.6 Conclusion 

The question of boundary bound states even in the simple Dirichlet case appears rather 
difficult: using the XXZ lattice regularization or equivalently the Thirring model, we have 
only been able to recover the and \Sn=o,N) boundary bound states. A way out is to 
consider solutions of the Bethe ansatz equations made of an {n, N) string superposed with 
the rz* + 1 string that describes the ground state. This is not allowed in principle in the model 
we consider because the Bethe wave function vanishes when two roots are equal. However, 
putting formally such a solution in the equations gives the masses of the \Sn,N) states and 
the S-matrix also agrees with the bootstrap results! But the meaning of this is not clear to 
us. 



4.7 Remarks 



Finally let us mention that one can calculate the ground state energy in the thermodynamic 
limit by solving the equation ( |4.44| ) for the ground state density and using ( |4.3CI| ): 



E. 



gr 



+00 



2Lpa{a)h{a)da. 



As a result, we get the combination Egr = E^uik + Eboundary, where Ebuik is the well-known 



sine-Gordon ground state energy p4 

Ebulk 



tan 



777 
2(7r-7) 



61 



and Ehoundary IS the contribution of the boundary terms (A — > oo): 



m 

'boundary 



sm 



(H+2)77r 
2(7r-7) 



+ 1 + cot 



.2 



sm 



2{7r-7) 



4(7r - 7) 



We see that the ground state energy of the boundary sine-Gordon model is a smooth function 
of the boundary magnetic field for the whole range of h in the XXZ regularization, hence of 
ifQ. The changes in ground state structure do not affect E, as is expected since in such a 
unitary model there is no (one dimensional) boundary transition. 

The finite size corrections to the ground state energy themshelves (the genuine Casimir 
effect) can be computed using the technique developed in [Q. We will show how such a 
calculation can be carried out in the next chapter. It is also interesting to consider the 
inhomogeneous 6-vertex model with an imaginary boundary magnetic field ensuring com- 
mutation with Uqsl{2) ||55|. This should lead to a solution of minimal models with integrable 
boundary conditions ||56| . 
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Chapter 5 



Boundary energy in integrable 
quantum field theories 



The main purpose of this chapter is to study the ground state energy of 1+1 integrable 
relativistic quantum field theories with boundaries. This involves several questions. One 
is the energy associated with a boundary for an infinite system, the other is the way the 
energy of the theory on an interval varies with its length - the "genuine" Casimir effect. 



The elegant method of Destri and de Vega |Q for the periodic systems leads directly to the 



expression for the ground state energy from which the infinite size contribution and the finite 
size correction can be easily extracted. The heart of the DDV method is a non-linear integral 
equation (|5.45|) being derived from the Bethe equations. We generalize the Destri-de-Vega 



method to the systems with boundaries and apply it to compute the ground state energy for 



the boundary sine-Gordon model. This chapter is a part of a more complete work |]56 



5.1 TBA for a QFT defined on cylinder 

For a Quantum Field Theory defined on a torus, the standard way to compute its ground 
state energy is through the Thermodynamic Bethe Ansatz p[. If the theory is defined on a 
circle of circumference i?, one switches to a modular transformed point of view where now 
the theory is defined on a circle of very long circumference L and at temperature T = 1/R. 
The free energy F of the theory in the "i? channel" can be computed using TBA, and it 
is simply related to the ground state energy E^{R) of the theory in the "L channel" by 
F = —TLE^[R). The limit R —>■ corresponds to the UV limit and is described by a 
conformal field theory. It is known that E^{R) behaves as E^{R) ~ nc/GR in this limit p7| , 
where c is the central charge of the underlying conformal field theory [Q.Q 

Consider now a quantum field theory on a cylinder of finite length R and circumference 
L, with some boundary conditions (a, b) at the ends of the cylinder. 

^ More precisely, one has 
where (A, A) are the scahng dimensions of the ground state. 
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R-Channel 



L-Channel 



R 




B 



R 




B, 



Figure 5.1: Cylindrical geometry in the R and L channels. 



Depending on quantization scheme, we have two possible ways to compute the partition 
function of the system (see figure 5.1). The first possibility consists in choosing as direction 
of time the horizontal axis and therefore the partition function will be expressed as 



■lab 



Tr e 



-LH„, 



(5.1) 



where Hat is the Hamiltonian relative to the system with boundary conditions (a, 6). In 
the second method the time evolution takes place along the vertical axis and therefore the 
partition function is given by the matrix element of the time evolution operator between the 
boundary states, i.e. 

Z,,= {B,\e-^'' \B,) = e^^, (5.2) 

where now H is the Hamiltonian of the bulk system. The ground state energy Eab{R) is, by 
definition, the leading term arising in the large L limit of the first expression, eq. (^.1|), i.e. 



Zab ~ e J 



(5.3) 



However, in view of the equivalence of the two quantization schemes, we can compute this 
quantity by looking at the large L limit of the second expression, eq. (|5.2| ). In the large R 
limit, R oo, the partition function (|5.2|) becomes 



(5,|0)e-^^"(0|i?,) 



(5.4) 



where |0) is the ground state of Hamiltonian H on the circle. It is the only state that 
contributes in the limit i? — oo; other states do not propagate along the cylinder. The scalar 
products Qa = {0\Ba) and = (0|-Bb) are called the ground state degeneracy. The quantity 
S{0) = lia{gagb), corresponding to the zero-temperature entropy, is one other universal term 
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that appears in the expansion of the free energy For periodic boundary conditions one 
must have gaQb = 1- 

Since eq. (|5.2|) employs the boundary states of the model, let us shortly recall their basic 
properties (for more detail see the original reference [|I|). In the QFT description of the 
model, the information on boundary conditions is encoded into a boundary state | B), which 
for infinite length L reads M: 



B) = gexp 



2^ 



K{9)A^{-9)A^{9) 



0). 



(5.5) 



Here g is an overall normalization related to the boundary entropy. For simplicity, we ignore 
possible additional contributions to the boundary state from zero momentum particles. From 
the point of view of QFT, the boundary state can be therefore regarded as a particular state 
of the Hilbert space of the bulk theory, made of a superposition of pairs of particles of equal 
and opposite momentum ("Cooper pairs"). All information relative to a particular boundary 
condition is encoded into the function K{6) which can be seen as the elementary amplitude 
to create a virtual pair of particles. 

In this chapter we address the question of computing the ground state energy of the 
sine-Gordon model with Dirichlet boundary conditions. Instead of deriving a TEA in the R 



channel, we adopt the beautiful approach of Destri and de Vega [Q. This requires working 
for a while with the lattice theory, here chosen to be the XXZ model with boundary magnetic 
fields. 



5.2 TBA for the inhomogeneous XXZ model with bound- 
ary fields 

like in the previous chapter, we start from the inhomogeneous 6-vertex model with boundary 
fields h gularization of the boundary sine-Gordon model with Dirichlet boundary 

conditions, with the only difference that now we consider the antiferromagnetic regime. 

In the inhomogeneous antiferromagnetic 6-vertex model with anisotropy parameter 7, one 
gives an alternating imaginary part ±zA to the spectral parameter on alternating vertices 



i6| , P7[1 . The scaling limit is given by taking A 00, N —> 00, and the lattice spacing 
A ^ 0, such that R = NA remains finite. In the bulk, this provides a regularization of the 
sine-Gordon model with Lagrangian 



R 



dx 







(5.6) 



where /i oc exp(— constA), = 8{n — 7), and the field is fixed at x = and x = R 
(Dirichlet boundary conditions) to a value that is simply related to h (see previous chapter, 
eq. (p3)). 

The wave function of the inhomogeneous six-vertex model can be expressed in terms of 
a set of "roots" aj, where j = 1 . . .n. They must be solutions of the set of equations ( |4.29| ): 

[/(«, + A, 7) + fia, - A, 7)] + 2/(a„ ^H) = 
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27r^j + [/("i ~ ^rn, 27) + f{aj + am, 27)] , (5.7) 

m=l,m^j 

where Ij is an integer and all aj are positive. The function / is defined as 

f{a,b) = 2tan^^ ^cot-tanha^ 

and 

= -/(27,-iln(/i + cos7)). (5.8) 

7 

By construction of the Bethe-ansatz wave function, aj > 0. Even though there is a solution 
of (|5.7|) with one vanishing root for any and n, we emphasize that aj = is not allowed 
because the wave function vanishes identically in this case. Observe that equations ( |5.7] ) are 
formally satisfied as well by the opposite of the roots, — a^ . Often in what follows we shall 
consider that the roots take both signs in order to rewrite equations in a way which is similar 
to the bulk case. 

For simplicity, we restrict to the case 7 = f where t is an integer, and restrict to the 



choice e = — 1 In the sine-Gordon model, this falls in the repulsive regime. We make 
the standard assumption that all the solutions of interest are collections of "fc-strings" for 
= 1, 2 . . . t — 1 and antistrings a. A fc-string is a group of aj in the pattern a^^^ — i7i{k — 
1), a^^-* — iTT(k — 3), ... , a*-'^^ + i7c(k — 1) where a'^^^ is real. The antistring has aj = a^""^ + iir, 
where is real. 

The thermodynamic limit is obtained by sending — 00. In this case, we can define 
densities of the different kinds of solutions. The number of allowed solutions of ( p.7| ) of type 
k in the interval {a, a + da) is 2N{pk{a) + p^{a))da, where pk is the density of "filled" 
solutions (those which appear in the sum in the right-hand-side of ( |5.7] ) ) and p^ is the 
density of "holes" (unfilled solutions). The densities pa and are defined likewise for the 
antistring. The "bare" Bethe ansatz equations follow from taking the derivative of (|5.?1 ). 
For 7 = vr/t they can be written in the form: 

*-i . 1 
2vr(pfe + p^) = ak{a) - (pk,t~i * Pa + J2 'i^^i * Pi + l^'^k 

1=1 

Mpa + P'a) = Mpt-l+Ptl) + l^{Ua-Ut_^) (5.9) 

These densities are originally defined for a > 0, but the equations allow us to define pk{—Oi) = 
Pk{a) in order to rewrite the integrals to go from —00 to 00. The kernels in these equations 
are defined most easily in terms of their Fourier transforms 

/•OO fjfy f POO 

f{x) = / — e*'^*"/V(a), /(«) = - / e-*"*"/V(a;)rfx. (5.10) 

J-00 277 TT J~oo 



One has 



3 coshx sinh(t — A;)x sinh /x 

(Pkiix) = Oab-2 — —7 , (5.11) 

smh X smh tx 
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for k > I with 



ak 



^kh and 

sinh(t — k)x 



sinh(t-2A;)x/2 



cos htX j 71 

sinn tx 

sinh(t — sinh kx 

2 -|- 

sinh X sinh tx sinh tx/2 

sinhifx sinh(t — 2)x/2 



(5.12) 



sinh tx 



sinh tx/2 



with in particular ai (a) = I /(a + A, 7) + /(a — A, 7) , 0ii(a) = — /(a, 27). The bound- 
ary manifests itself in the first term in Uk, notice that even for /i = 0, it still modifies the 
equations. A few technicalities account for the other terms (these are relevant here because 
we are interested in subleading boundary effects). The second term in Uk arises from the 
fact that the sum in ( |5.7| ) does not include the term m = j; the integration over densities 
includes such a contribution and so it must be subtracted off by hand. The third term in Uk 
arises because p and are defined for allowed solutions, while as already explained, a = is 
not allowed because it does not give a valid wavefunction. Since it is a valid solution of 



but is not included in the densities, we must subtract an explicit ■^6{a) (corresponding to 
1/2N in Fourier space). Explicitly, one has 



Ml = 2f{a, -fH) + 2/(2a, 27) - 27iS{a). 
For compactness we rewrite ( ^.9] ) as 



27ia^''^ {pk + pI) = ak{a) + ^^i * pi 

I 

where a'^^'^ = — 1 for the antistring. 

The energy reads, with proper hamiltonian normalization, 

j^latt 



2N 



Uk, 



(5.13) 
(5.14) 



2N 



1 f°° 

'1^ I ak{a)pk{a)da. 

t J J — CO 



(5.15) 



It is easy to write the thermodynamic Bethe ansatz for this model. One finds that the 
TBA equations, since they are obtained by a variational method, do not depend on boundary 
terms, and read as usual 



■-ak{a) = Tin (1 + e^^) - T V * In (l + e 

f ATT V 



where 



Aki{a) =27ia^''^Ski5{a)-^ki- 
The free energy does depend on the boundary term and reads 



(5.16) 



(5.17) 



tny: 

k 

T 



or^'=^a,(a)ln 1 + e-^M — 

00 ^ ^ /TT 



da 



-00 ^ 



da 
2^' 



(5.18) 
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In the above formulas the temperature T corresponds in the two-dimensional point of view 
to having a cylinder of radius L = 1/T. We can deduce from this result the ground state 
energy. Indeed recall that the ground state is obtained hj pk = 0, k ^ 1 and pj* = so 

^latt ^ ai(a)|e^| / ui|e^| — , (5.19) 



2tt 2 j-oo 27r 



(5.20) 



where from ( 5.16|) we have 

1 cos (Atx/vr) 

Ci = . 

t cosh X 

Replacing and using ( |5.12| ) we find 

^ /Atx\ sinh(t-l)x 

Etl = — / cosM — ^ . ; 

TT j-oo \ Tx J coshxsmhrx 
r^iatt ^ 1 z-"" cos(Atx/7r) / sinh(t - ij')x sinh(t - 2)x/2 \ 
^^'y 2ttJ^oo cosh a; 1^ sinhte sinhtx/2 ; ^ ' 

where we used the formula 

a{a)b{a)da = 2t / a{x)h{—x)dx. 



In the continuum limit A — > oo the energy contains various terms. We keep only the finite 
part which is obtained by closing the above integrals in the upper half plane and selecting 
the pole at X = i|, leading to 

Ebuik = R— cot — 




, sin(t - m-K 2 tn \ , , 

Ebdrv = 2 — ^ cot 1 , 5.22 

^ ^ * sint7r/2 4 / ' ^ ' 



where m is the soliton mass, 



m 



All these results trivially generalize to the case of two different boundary fields by splitting 
the H dependent terms into the sum of an H and an H' term. The bulk result agrees with 
what is obtained by other methods. As in the bulk case, when t is even, there are additional 
logarithmic terms. The boundary entropy is actually the same in the UV and IR limits. 



5.3 The Destri De Vega equations for the 
boundary sine-Gordon model 

We now would like to compute the complete Casimir effect in a theory with boundary. TBA 
in the R channel is pretty intricate in the non diagonal case. We adopt an alternative method 
elaborated by Destri and De Vega (DDV) in the periodic case. 
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5.3.1 The DDV equations with boundary conditions 

Consider eq. ( |5.7] ) which we rewrite as 

2Np{aj) + Pbdryiaj) + ^ (f){aj - am) + (f>{aj + am) = (5.24) 

am>0 

where the sum runs over all roots (including m = j) and we introduced the notations 

= ^[/(a + A,7) + /(a-A,7)], Pnia) = f{a,-fH), 0(«) = 0n(«), (5.25) 

Pbdryia) = 2pHia) - 0(2a). (5.26) 

The ground state is obtained by filling the real positive solutions, aj = excepted. This 
corresponds to the choice nj = 1,2,.... Recall that if aj is solution of (|5.24|) with some nj, 
so is formally —aj (with —fij). Given the set of roots {aj > 0} representing the ground 
state, one can construct the counting function as follows: 

f{a)=2iNp{a)+ipbdry{a)+i 0(a - "m) + 0(a + "m), (5-27) 

am>0 

Define then 

F(a) = e^(°). (5.28) 

We have Y{aj) = Y{—aj) = 1 for every root aj from the ground state, as well as Y{0) = 1. 
Therefore we can rewrite ( |5.27| ) as 

r Y((y') do' 

fia) = 2tNp{a) + tp^ryia) - ^0(a) - / 0(« - a')-^-^ — , (5.29) 

Jc 1 — y [a ) Zir 

where i(f){a) at the right-hand side takes care of the unwanted contribution of the pole a' = 
and the contour C consists of two parts as shown in figure 5.2, Ci above and C2 below the 
real axis. 

Like in the bulk case 15^, simple manipulations allow us to rewrite this non-linear integral 
equation as 

f{a) = 2iNP{a) +iPMry{a) + [ ^{a - a')\n (l - e^^""'^) da' 

+ / $(«-a')ln(l-e~^("'^)c/a'. (5.30) 



In equation (|5.3CI|) one has 

= / ^-rp-**^"/" 

2 sinh(t — l)x cosh x 



*(«) = - A r cixe--/.^^^^^^^^^^, (5.31) 
^ ' 2tt^J-oo 2sinhft - Dxcoshx' ^ ^ 

together with 



, , f+'=° , e"^*^"/'' - 1 cos(Atx/7r) 

P(«) = dx : ^ ' ' 5.32 

J —IX 2 cosh X 
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C1 



C2 



Figure 5.2: Contours of integration in the DDV metliod. 



and 



bdry 



-hoo g- 

dx— 

oo 



-IX 



^ smh.{t — 2)x/2 coshtx/2 ^ 



sinh(t — H)x 
sinh(t — l)xcoshx 

sinh(t — 2)x 



sinh(t — l)x cosh x 2 sinh(t — l)x cosh x 



(5.33) 



Obtaining ( ^.30|) requires some care with the definition of logarithms. One proceeds as 
follows. Before integrating by parts in the integral over C2 one factors out -^hiY{x): 



Yix] 



-^in[i-r-i] + -^iny(x). 



1 - Y{x) 

Then both integrals over Ci and C2 can be taken by parts, resulting in 



/(«) - 



+°° • da' 

(f){a - a')f{a') — = 2iNp{a) + ipbdry{a) - i(f){a) 

-00 ZiT 

+00 . /Vq,' 

a - a' - zO) ln[l - Y{a' + iO)] — 



+ 



00 

-foo 



27r 



{a-a' + iO) ln[l - - iO)] — , 

2n 



(surface terms from Ci and C2 cancel against each other provided N, A are finite). To make 
the source term vanish at infinity we take the derivative of both sides in the latter equation, 

after which it can be Fourier-transformed and "dressed" by the factor (1 — 0)^^. Finally, 
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one goes back in the rapidity space and integrates the equation, using /(O) = 0, to obtain 
The energy of the ground state configuration can be expressed as 

E'"'' = -7 EpK) = -7E/K -A,7)+/(-«.-A,7) 

where fla,^) = /7(a). One finds after exactly the same manipulations as above, 
1 • , . , 1 /■+°° „ „ da 



E 

t 



A(A) - 7 / - A + *0) ln[l - Y{a + zO)]— (5.35) 



Substituting ( p.30| ) instead of /(a) in the last term of the latter, we obtain 



1 r+co I 

j^latt ^ j^latt^ ^ j^la^t^ _ / _ A + ,0) ln[l - Y {y + ^0)]- 

i J~co 

+ ' r°^s{y-A-zO)\n[l-Y~\y-zO)]^, (5.36) 
t J-oo Zn 



where we defined s{y) by 



kdk ^_^^y _ tanh(ty/2) 



s(y) = i / e-^'y = , (5.37) 

' 2cosh7fc 4 cosh(t?//2) ^ ^ 



— oo 



The last two terms in (|5.36|) represent finite-size corrections to the ground state energy, while 



j^iatt . j^iatt _ 1 / ^ Z^" dk 2Np{k)+p,,^y{k)-<P{k) 

t tJ-oo 27r 27T-(j){k) 

1 • 2N /■+°^ • 

= f^{A)- f^{a-A)p{a)da, (5.38) 

t t J —oo 

where the function p{a) defined so satisfies the following equation: 

27ip{a) = p{a) + * p{a) + ^ {pbdry{a) - 0(a)) , (5.39) 

which can be checked by solving this linear equation in Fourier space. Introduce pi{a) = 
p{a) - 5{a)/2N. Then 

ON /•+00 . 

<Ji + Eify = -— / Ua - A)p,ia)da (5.40) 



and, by virtue of (|5.39| ), pi satisfies the equation 



2npi{a) = p{a) + * Pi(a) + ^ {Pbdry{a) - 2vr(5(a)) . (5.41) 



and E\^^^ coincide with the quantities computed in the previous section. 



Hence p\ is the density of the ground state configuration ( p.9| ) (see also ( [5.13| )) and E\^^l^, 

:platt 
-'bdry 



71 



5.3.2 The continuum DDV equations 

Having checked the correct values of bulk and boundary energies, we now let the cutoff 
A — s> oo according to A = | log(2/mA) with the size of the system R = NA and the physical 
mass m = ■^e"*'^/^ fixed, and work only with the renormalized theory. In that limit one has, 
evaluating all integrals in Fourier transforms and keeping the leading terms, 

s{a + A) + s{a - A) = ~—msmh.9, NP{9) = mRsiiih9, (5.42) 

where we set 

= y . (5.43) 

Recall indeed that when one studies the excitations of the model, a relativistic dispersion 
relation is obtained provided the rapidity in the relativistic theory and the "bare rapidity" of 
the Bethe excitations are related by ( 5.4ij| ). We redefine implicitely all functions to depend 



on 9 from now on. The energy therefore reads now, 

d9 



1 r 

E = Ebuik + Ebdry ^ 2 Jc ^ ^^^^ 6' In (l - 



+ - f msmh9\n(l- e'f^'^A —, (5.44) 
2 Jc2 ^ ' 2m 

where / is solution of the integral equation 

j[9) = 2imRsmh9 + iPMry{9)+ [ <^{9 - 9') In (l - e^^^'^) d9' 

+ [ ^9-9')\n(l-e-f^^'Ad9', (5.45) 

JC9. 



and 

J-oo 21^^ cosh X smh(t — \)x 



which can be identified with 



where 5*++ is the soliton-soliton S matrix element (see eq. (|1.37|) ) 



5.3.3 The Casimir effect 

We define the effective central charge by the formula 



E = Ebulk + Ebdry — (5.48) 

Z4ri 
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From (|5.44|) we find 



CeffimR) = ( f msmh9ln(l-e^^^Ade+ f msinh^ln (l - e^^^^^) j (5.49) 



To study the ultraviolet behaviour of the above expression, let us use f{9) = f{6) and rewrite 
as 



f{e) = 2imRsmh9 + iPbdryiO) - 2i / - 9') Imln 

and (|5.49| ) as 



[mR) 



12mR 



1 - e 



ofie+io) 



f{e'+io) 



(5.50) 



(5.51) 



dOsmhO Imln 

It might be useful to remind the reader of the form of the corresponding bulk equations p4 

QmR 



CeffimR) 



sinh^ Imln 



1 + e 



fie+io) 



d9. 



with / satisfying 

f{9) = imR sinh 9 + iu-2i I $(0 - 9')lm In 



1 + e 



f{e'+io) 



d9\ 



(5.52) 



(5.53) 



and uj is the twist of the 6- vertex model. Note the deep similarity between these two systems; 
the factors of 2 can obviously be absorbed in a redefinition of mR and the minus sign in 
the arguments of logarithms in a redefinition of /, so the only essential difference is that the 
twist angle uo is replaced by Pbdry It is well known that the twist corresponds to a soliton 
fugacity ||3^ so we see that the boundary acts by some effective, rapidity dependent fugacity. 

In the limit when i? — >■ 0, only the region |^| large contributes to Cg//. Let us focus on 
the limit ^ 3> 1, the results for negative 9 following by symmetry. Then one finds f ^ fx, 
where 



fKi9) = iniRe^ + iPbdryioo) - 2% \ ^{9 - 9') Imln 1 - e^^(^'+^°) d9' , 
together with 



6 

CxifTiR) = ^ I mi?e ''imln 



1-e 



fK{e+io) 



d9. 



H J-00 

It is now useful to recall some well known results about dilogarithms p2l. Define 



L(x) = 



du 



ln(l + u) \nu 



u 



l + u 



Assume 



/oo 
dy G{x - y) Imln [1 + F{y + iO)] 
-00 



(5.54) 
(5.55) 

(5.56) 
(5.57) 



73 



with G an even function. Then one has 



Im / dx^'ix) In [1 + F{x + le)] = -Re {L[F(-oo)] - L[F{oo)]} 

+hm {(j){oo) In [1 + F(oo)] - 0(-oo) In [1 + F(-oo)]} , (5.58) 

(where we did not write the iO part of some arguments for simphcity). Set F = e^^~'^'^ and 
denote PMry{.oo) = a. Then, according to (|5.54| ): (j) = mRe^ + a — n, G = — $. We have 
(j){—oo) = a — TT, (f){oo) = oo. One has also F{oo + iO) = 0, and from this and ( p.54| ) one can 
get the value of F at — oo: 



/i^(-oo) 



a — 2 Im In 

t - 2 



t-1 



Imln 



So, if e^J^(-°°+*°) = e'^ one may use Imln(l ± e'' 



ln(±e*^) to find 



2i 



t — 1 TT 

exp <j 2z — - — a + 2i— ^ . 



(5.59) 



(5.60) 



From ( |5.58| ) it follows that 



Im 



niRe^ ln(l - e^^^'(^+*°V^ = - ReL(-e^") - - (tx - vr) 



In the region ^ 1 we have f ^ Jai and similar calculations yield 



Im 



mRe ^ ln(l — e 



fA{e+io) 



)de = -i ReL(-e-^'^) + ^ (tt - a) 



t-1 TT 

(T — vr H — 

t t 



t-1 TT 

a + 71 



Collecting both ^ ^ 1 and 6 <^ 1 contributions we obtain 



6 ri 



6 



TT 

~6 



L(-e2'^")) + L(-e 
t- 1 



t - 1 



(T - TT 



(a - tt)^ 



From (|5.33| ) we get 
Finally, from this we find 



Pbdryioo) = a = 27T - -'- 



7iH + H' 



2 t-1 



Cjjv = 1 - I 1 



t \ 2{t-l)) ' 
In the case with no boundary field, H = H' = t — 1 so cuv = 1 as expected. 



(5.61) 



(5.62) 



(5.63) 
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Remark. So far we tacitly assumed that < H < t — 1. In general, from relation 
it follows that when h > 0, H varies between — t — 1 and —1, while when h < 0, — 1 < 
H < t — 1. To generalize our results, we should use the most general form of p//: 

P^(k) = r pH(a)e''"'da = 2n sign(H f''^^^''~ "t"^^^ , -tt < < tt, (5.64) 
j-oo smh 71 k 

where we defined ljh = \'yH\. For — 27r < ■jH < —n set uh = 2it + •jH and sign(if) = 1 in 
( [5 .641 ). Then ( p.62| ) generalizes to 

2 TT — 7 

if H and H' are both positive or —2t < {H, H') < —t, and 

TT An - ujh - uj'h 
a = Zn 

2 TT — 7 

if they are both negative, but greater than —t. In the case when < H < t — 1 and 
-t < H' < -1 (that is, h < 0,h' > 0) we get: 

n uj'jj - ujh - 2'y 

O" = TT H . 

2 TT — 7 

So, when cu^ — cj// = 27 we have a = 71 and C[/y = 1, as in the free case. The condition 
uj'fj — ujh = 27 is equivalent to h = —h', as could be seen from (|5.8| ). That c = 1 when the 



two surface field are real and opposite is well known from lattice studies pO 



5.4 Remarks 



A particularly interesting case is when the XXZ chain or the inhomogeneous 6-vertex model 
commutes with the quantum group Uqsl{2). In that case h = —h' = 2isin7 and the net 
result is that all H dependent terms simply disappear from the equations, so, in particular 



bdry 



m f t7T 

— cot — 

2 V 4 



(5.65) 



At the N = 2 supersymmetric point, t = 3, the boundary energy vanishes, a result well 
expected from supersymmetric considerations. More generally, it vanishes if t = 4n + 3, n 
an integer. Notice that the bulk energy vanishes for t an odd number (as a consequence of 
the generalized fractional N = 2 supersymmetries studied in [|6H] ). 

In the quantum group symmetric case |^ one has H + H' = 2t so from (| 



1 - 



6 



t(t-l)' 

the expected result for the restricted sine-Gordon model [^, 



(5.66) 
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Chapter 6 

Surface excitations and surface energy 
of the antiferromagnetic XXZ chain 
by the Bethe ansatz approach 



We study an open XXZ chain in the regime A > 1 with a boundary magnetic field h and 
discuss some of its pecuhar features due to the presence of boundary. In the Bethe ansatz 
formahsm, boundary bound states are represented by the "boundary strings" as described 
in chapter 4. We find that for certain values of h the ground state wave function contains 
boundary strings, and from this infer the existence of two "critical" fields in agreement with 
49(| . An expression for the vacuum surface energy in the thermodynamic limit is derived 



and found to be an analytic function of h. We argue that boundary excitations appear only 
in pairs with "bulk" excitations or with boundary excitations at the other end of the chain. 
The case where the magnetic fields at the left and the right boundaries are antiparallel has 
non-trivial differences with the case of the parallel fields. The Ising (A = oo) and isotropic 
(A = 1) limits are discussed thoroughly and found helpful for the intuitive understanding of 
the behavior of the boundary XXZ chain at arbitrary A. This section is based on the work 
6l. 



6.1 Introduction 



In this chapter we study the XXZ chain with even number of spins L in a boundary magnetic 
field, 



(6.1) 



in the regime A > 1, 
boundaries lESl. At hi 



hi > 0, h2 < 0, focusing on the effects peculiar to systems with 

= h2 = this model describes one-dimensional antiferromagnet with 

non-magnetic impurities, accessible experimentally. We exploit the Bethe ansatz solution 
for this model, first derived in ||^, together with the well-known results for the periodic 
chain [^, We find new "boundary string" solutions to the Bethe equations, similar to 
the boundary strings existing in the | A| < 1 regime |^0[. For certain values of the boundary 
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magnetic field the ground state configuration contains boundary 1-strings. Boundary exci- 
tations are obtained by removing (or adding, depending on the sign of h) boundary strings 
from the ground state wave-function. Their energy was first obtained in |^ by the algebraic 
approach. 

A pecuhar feature of the Bethe ansatz solution of the periodic chain is that the excitations 
(holes in the Dirac sea) appear only in pairs [Q. We argue that similarly the boundary 
excitations can appear only in pairs with bulk excitations or with boundary excitations at 
the other end of the spin chain. There is no such restriction in the soultion of the semi-infinite 
chain by the algebraic approach 

Using the Bethe ansatz solution we calculate the surface energy (see e.g., [|6^] ): 

EsurfiL,A,h)=Egr-E''g„ (6.2) 

in the thermodynamic limit L = oo. Here Egr is the ground state energy of ( |6.1| ) and E^^ is 
that of the periodic chain. We give an interpretation of our results in the limits A oo and 
A ^ 1, corresponding to the ID Ising and XXX models respectively. Finally, we comment 
on the structure of the ground state when the boundary magnetic fields are parallel. 



6.2 The Bethe ansatz equations 



Let us first set up the Bethe ansatz (BA) notations and list the relevant results about the 
XXZ chain ||4^, In the eigenstates of (|6.1|) were constructed for arbitrary A. As 
usual in the BA picture, the n-magnon eigenstates |n), satisfying T-C\n) = E\n), are linear 
combinations of the states with n spins down, located at sites xi, x„: 



\n) 



The wave-function 



/(Xi, ...,Xn) =J2^pMPi, ■■■,Pn)e 



i{piXl + ...+p„Xn) 



(6.3) 



contains n parameters pj G (0, vr) which are subject to quantization conditions, called Bethe 
equations (BE): 



^2iLpj 



e'P^ +hi- A e'P^ + hi - A 



n 



(6.4) 



1 + (/ii - A)e'P^ 1 + (/i2 - A)e'P^ 

The summation in ( |6.3|) is over all permutations and negations of pj. The energy and spin 
of the n-magnon state are given by p5| : 

1 ^ L 

E = -[{L-l)A + hi + h2]+2Y,{cosp,-A), S, = --n. (6.5) 

^ 3=1 ^ 



It is convenient to rewrite BE using the following mappings: 

A = cosh 7 > 1, 7 > 0, 



(6.6) 
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p 



-i In 



cosh ^{ia + 7) 
cosh ^{ia — 7) ' 



(6.7) 



(our definition of p{a) differs from that of ||2y, ^ by the shift a a + n and it was chosen 
in such a way that p{a) be an odd function that maps —n < a < n to —n < p < n), 



h 
h 



^"^^^+ sinhi(l + i/) 
coshm-H) 



7 



sinh(7) ■coth^(i7 + l) 



7 



sinh(7) ■ tanh — {H + 1 



hum < \h\ < 00, (6.8) 
\h\ < hum- (6.9) 



The latter two mappings are defined on H E {—00, 00) and are necessary to cover the region 
—00 < h < 00, with positive h corresponding to H E (— l,oo). The value hum = ^(00) = 
sinh7 lies between two critical fields h^J^\ hf) defined as follows [HOl: 



/i« = A-l, M2) = A + 1. 



(6.10) 



Both critical fields correspond to = 0, and the gap h'^^) < h < hf) corresponds to 
< < 00. In these notations eq. ( |6.4| ) becomes: 



cosh i {iaj + 7) 
cosh \ {iaj — 7) 



2L 



B{aj,H^)B{a,,H2) 

sinh \ {iaj — iam + 27) sinh \ {ioij + idm + 27) 



n 



sinh \ {ioij — io-m — 27) sinh \ {ioij + idm — 27) ' 



(6.11) 



where 



B(a,H) 
B{a,H) 



cosh ^{ia + •yH) 
cosh ^{ia — •yH) ' 
sinh ^{ia + yH) 
sinh \ {ia — yH) ' 



^Zim <\h\ <00, 
\h\ < hum, 



(6.12) 
(6.13) 



are called boundary terms. The energy eq. (|6.5| ) takes the form: 



E = - [{L-1) cosh 7 + hi + h2] - 2 sinh7 ^p'(aj) 



P(aj 



sinh 7 



cosh 7 + cos a 



(6.14) 



In the thermodynamic limit L ^ 00 the real roots aj of BE form a dense distribution in 
the open interval (0, n) with density p(a), dl = 2L{p + ph)da being the number of roots in 
the interval da. The logarithm of eq. ( |6.11|) is: 

11 " 
2Lp{aj) + -\nB{aj,Hi) + -\nB{aj,H2)+^{2aj) = Y,^{aj-ai)+(f){aj+ai)+27rlj, (6.15) 
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where Ij form an increasing sequence of positive integers, and 

/ N , sinh ^(27 + ia) , , , , , 

Taking the derivative of eq. ( |6.15| ) and defining p for negative a by p{a) = p{—a), we obtain 
+ = J%'ia-P)p{P)dp + 2n{p{a)+p,{a)), (6.17) 

with 

« = -^ p) rr ! - ^ r.) rr i + 20' 2a - 27r5 a - 27r5 a - vr . 6.18 
^ B{a,Hi) B{a,H2) 

The presence of delta-functions in (|6.18| ) is due to the fact that aj = and aj = it are 
always solutions to ( |6.11| ), which should be excluded, since they make the wave-function 
(|6.3|) vanish identically. 



6.3 Solution for the bulk part 



In eq. ( |6.17|) the "boundary terms" are down by the factor 1/2L. Neglecting p[^j.y 
setting ph = 0, we obtain the equation for the ground state density of the periodic XX Z 



chain ||6^. Solving it by the Fourier expansion 

/(«) = E /(Oe^'", /(O = TT / 

/ ^ ^TT J-7T 



(6.19) 



and using ( |6.14] ), we recover the result for the ground state energy of the periodic chain pB 



E 



2Trpperin) 
LA 



p [n) 



l + d'in) 



gr 



2Lsinh7 / pper{(y)p {(y)da 



-27|n| 



LA 



p [n) 



;-l)"e 



n -7|n| 



L sinh 7 E 



-7|n| 



cosh 7n 



(6.20) 



(6.21) 



The spin of the ground state is Sz = L/2 — L J^^ pperda = 0, which is the well-known result 
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An elementary "bulk" excitation above the vacuum in the model (|6.1| ) is a hole in the 
distribution of Ij, but only a pair of holes can occur for the periodic chain, as argued in |]n . 



Thus physical excitations contain an even number of holes. The energy of the hole with 
rapidity 9 can be easily computed: 



e,(e) = sinh7 E >0' 

n^oo cosh7n 



(6.22) 



79 



and the spin with respect to the vacuum is S'^ = 1/2. (Our result, eq. ( |6.22D , differs from 
the conventional one by the shift 9 + tt, but the dispersion relation is unchanged by 
rapidity reparametrization.) 

Analogous arguments can be applied to analyze "bulk" string solutions with complex 
values of a. Although there exists an infinite hierarchy of complex strings of arbitrary 
length, and quartets, their energy vanishes with respect to the vacuum pO|. 



6.4 Boundary excitations 

So far we discussed the bulk excitations, which are essentially the same as in the periodic 
chain. Let us now turn to the new solutions of eq. ( |6.11|) , boundary strings. The analysis 



is close to that of section 4.4. Boundary excitations have their wave-function ( |6.3| ) localized 
at the left or right ends of the chain, and in the limit L ^ oo the two boundaries may 
be considered separately. Let us study first the left boundary, /ii > 0. The fundamental 
boundary 1-string consists of one root located at = —i'jHi for < hi < h^}^ ^ and at 
ao = TT — i'~^Hi for h^^) < hi < oo (in both cases — 1 < f^i < 0). It is a solution of 
BE due to the mutual cancellation of the decreasing modulus of the first term in ( |6.11|) 



and the increasing modulus of the second term B{a,Hi) as L ^ cxd and a ao- When 
h^J^^ < hi < h^^\ no such solution exists. Introduction of such a string into the vacuum with 
the density of roots p{a) defined from 

+ = r - l3)pil3)d/3 + 27rp(a), (6.23) 



2^j^jr odry \ 

leads to the redistribution of roots by 6p = 2L(p — p) satisfying the equation 

= r 0'(a - p)6p{P)dp + 0'(a - Oq) + 0'(a + "o) + 27r5p. (6.24) 

J —IT 

Prom the latter we find 

, 2cosnaoe-2'>l'^l 

2^^P(^) = 1 + • (6.25) 

For the energy and spin of the boundary 1-string with respect to this vacuum we obtain: 

£6 = — 2 sinh(7)p'(Q;o) — sinh7 / 5 pp (a) da = — sinh.'y , (6.26) 

J~n „_ „ cosho'n 



S, = -- . 
2 



cosh 7n 

1 



We see that the energy ( |6.26| ) is negative, so the correct ground state wave-function 
should contain the boundary 1-string root when the value of hi is not in the gap h^^} < 
hi < h^^\ The ground state density p in this case satisfies the equation 



+ -«o) -</>'(" + ao))= / _0'(a-/?)p(/?)(i/? + 27rp(«). (6.27) 
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The boundary excitation is obtained by removing from vacuum ( |6.27| ) the root ao, which 
means that it has the energy —ib > and spin 1/2, equal to the spin of the bulk hole. Sub- 
stituting the value of ao into (|6.26|) , we get the boundary excitation energy, which precisely 
agrees with the one obtained in 

eb{h) = smhj ' -1<H^<0, (6.28) 

cosh7n 

with K = 1 if /ii < /i^J.) and k = 2 ii hi > hf) . 

From ( |6.22| ) and ( |6.28| ) we see that for hi < h^^) the energy of the boundary excitation 



is smaller than the bottom of the energy band of bulk excitations, and becomes equal to it 
at hi = /i^J.) (see Figure 6.1). So in this regime we may interpret the boundary excitation as 
the bound state of the kink, which gets unbound at hi = h^J^\ For hi > h^^) the energy of 
the boundary bound state is bigger than the top of the energy band. Therefore it is stable, 
in spite of its huge energy. 

Besides the fundamental boundary 1-string, there exists an infinite set of "long" boundary 
strings, consisting of roots ao — 2zA;7, ao — 2z(/c — 1)7, ao + 2m7 with n. A; > 0. We will call 
such solution an (n,k) boundary string (thus the fundamental boundary string considered 
above is the (0,0) string). One can use the same arguments as given in section 4.4 to show 
that the (n,k) string is a solution of BE when its "center of mass" has positive imaginary 
part and the lowest root ao — 2i/c7 lies below the real axis. Thus, sufficiently long boundary 
string solutions exist even in the region h'^}^ < hi < hf^ . However, a direct calculation shows 
that their energy vanishes with respect to the vacuum, so they represent charged vacua. Q 
(Analogous phenomenon occurs for the "long" strings in the bulk [^|: if the imaginary part 



of a lies outside the strip —27 < Ima < 27, the root a gives no contribution to the energy.) 
For Q < hi < h^^) and h^^) < hi < 00 the (n,0) strings also represent charged vacua, while 
(n,k) strings with k > 1 have the same energy (|6.28|) as the boundary bound state found 
above, and hence represent charged boundary excitations. Q 

Consider now the right boundary, h2 < Q {H2 < — 1). Now the fundamental boundary 
1-string solution ao = —i'~fH2 exists for any value of h2 in the interval —hum < h2 < (resp. 
ao = vr — i'yH2 for /12 < —hum)- Explicit calculation shows that it has non- vanishing energy 
only if —2<H2< —1, which corresponds to —hl}J < /i2 < (resp. /i2 < — /i^^-*). For such 
values of /i2 the energy of the 1-string with respect to the vacuum ( |6.23| ) is positive and 



equal to £&(— ^12) (see eq. (|6.28|) ), and its spin is = —1/2. In some sense the pictures are 



dual for the positive and negative h cases: there exist two states when \h\ is not between 
\h^cr\ and \h^^^\, one with boundary 1-string and one without. One of them is the ground 
state and another is the excited state at the boundary, and these states exchange their roles 
when the sign of h changes. The analysis of long boundary strings is very similar to that at 

^As an example, consider the boundary (l,0)-string consisting of the roots ao + 217, ao- It exists if 
— l<Hi< 1, although the (0,0)-string exists only if —l<Hi< 0. The (l,0)-string has charge 5*^ ~ — 1 
and vanishing energy with respect to the vacuum. 

^E.g., (1,1) string with roots ao + 2ij, ao, ao — has 5*2 = —3/2 and energy given by ( |6.28 ) with respect 
to the physical vacuum. 
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the left boundary, and therefore will be omitted. The net result is again that long boundary 
strings represent charged vacua or charged boundary excitations. 

In all examples shown above, the charge of boundary excitations turned out to be half- 
integer. One can easily check that this is true for all boundary strings representing charged 
excitations. Since the charge of physical excitations is obviously restricted to be an integer 
(see ( |6.5| )), we conclude that a boundary excitation can appear only paired with the bulk 
excitation of half-integer charge (i.e. containing an odd number of holes), or with a boundary 
excitation at the other end of the chain. We give a qualitative interpretation of this fact 
below. 



6.5 The surface energy 

To compute the vacuum surface energy, eq. (|6.2|), of the model (|0|), one should use eq. 
( |6.14|) in the limit L = oo with the root density determined from eqs. (|6.23|) or ( |6.27| ) and 
the boundary terms (|6.12|) or (|6.1ij|), depending on the value of h. Define for convenience 



giA) = - + 2sinh7^ 



oo g— 2n7 



n=l 



cosh 2n7 



(6.29) 



We consider separately the following intervals for positive hi and negative /12: 



1) 1^1,2! < hl^^. The ground state contains one boundary 1-string, corresponding to hi. The 
spin of the ground state can be found to be Sz = 0. Using eqs. ( |6.13| ), ( |6.18| ) and 
(|6.27|) , and subtracting the bulk contribution ( |6.21D we get 



E. 



surf 



1 00 p-jHin _ p-fH2n 

-{hi + h,)-g{A)~smh^J2{-ir 



n=l 



cosh 7n 



(6.30) 



2) |/ii,2| > ^cr''- The ground state contains one boundary 1-string and has Sz 
eqs. (EJT) and (|07|) it follows: 



0. From 



E. 



surf 



-{hi + h2)-giA)-smh^Y. 



00 ^-^Hin _ ^jH2n 



n=l 



cosh 7n 



(6.31) 



^) ^cr"* < 1^1,2! < hiim- The ground state has no boundary strings and its spin is zero. From 
( |6.23|) and ( |6.13| ) one obtains the same expression as in case 1. 



4) hum < \hi,2\ < h^l From ( [Ol and (|6A^ ) one obtains the same expression as in case 
2. The ground state has the same structure as in case 3. 

A qualitative plot of the surface energy as a function of h {h = hi = —h2) is given in Fig- 
ure 6.2. The apparent difference between ( |6.30| ) and ( |6.31D is an artefact of our parametriza- 



tion of h in terms of H. In fact, Egurf is an analytic function of h in the domain h G (0, cxd). 
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which can be seen after substituting if as a function of h according to ( |b.8D -( p3?D . In this 
sense the fields are not actually "critical." We find for hi = h2 = the value 

A fl °° p2n7 _ 1 oo (_l)n \ 

Note that one can obtain the boundary magnetization (af) |^ immediately from the formula 
for the surface energy ( p.30| )- (|6.31| ) by differentiating it in hi. 



6.6 The Ising A = oo and rational A = 1 limits 

In the extreme anisotropic limit A — > oo, /i ~ A of the XXZ chain ( |6.1|) one gets the 
one-dimensional Ising model: 




Aafa^+i + hia'i + /isfxlj , (6.33) 
In this limit from ( |6.8D -( |OD one has 



h^A± e-^", (6.34) 

and the gap between h^J^^ and h^^^ dissapears, so for any h there exists a boundary bound 
state. The energy of the "bulk" hole ( |6.22 ) becomes 6'-independent and equal to A, since 



only n = term contributes to the sum when 7 ^ oo. The energy of the boundary bound 
state (|6.28| ) becomes = A ± e~'^^^ = hi. This suggests the following interpretation in 



terms of the Ising chain. In the Ising ground state the i-th. spin has the value (—1)*. Local 
bulk excitation of the smallest energy 2A can be obtained by flipping one spin (the first 
and last spins excepted). The arising two surfaces (domain walls) separating the flipped 
spin from its right and left neighbours are called kinks and carry the energy A each. Kink 
corresponds to a hole in the Bethe ansatz picture, and kinks obviously appear only in pairs, 
which demonstrates that holes can exist only in pairs, too. The charge of the one-spin-fiipped 
state is equal to one, in agreement with the charge of two holes in BA. In addition to charge 
1 excitation, one has charge excitation of the same energy obtained by flipping any even 
number of spins in a row. In the BA this corresponds to the "2 holes and 2-string" state. 
In the Ising model the left (right) boundary bound state is obtained by flipping the first 
(last) spin. Such a state has the energy hi + A above the vacuum energy, where hi is the 
contribution of the boundary term in ( |6.33| ) and A is the energy of the kink created due to 
the boundary-bulk interaction. Thus flipping the boundary spin actually gives a combination 
of the boundary excitation and the bulk kink. Still another possibility is to flip all spins, 
creating two boundary bound states, one at each boundary. This explains why, in the BA 
picture, a boundary excitation can exist only if paired with a hole in the Dirac sea or with 
another boundary excitation. The vacuum surface energy ( |6.2[ ) of the Ising chain in the 
thermodynamic limit is {A — hi + h2)/2. The A/2 contribution here is the bulk interaction 
energy that we lost when we disconnected the periodic chain, and ±/ii^2/2 is the contribution 
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of each of the boundary terms. Taking the hmit 7 — (X) in eqs. ( |b.3(JD - ( |t).31 ), we obtain the 
expected resuh Egurf ^ {A — hi + /i2)/2. 

In the isotropic (rational) hmit A ^ 1 {i.e., 7-^0) one gets the XXX chain in a 
boundary magnetic field, which was discussed in the BA framework in for < hi^2 < 2. 
From (16. 81) -(16. 91) one has in this limit 



h 



l + H 



(6.35) 



There is only one critical field her = 2, which is the limit of h^\ Passing from summation 
to integration in eq. (|6.31|), we obtain for < /ii < 00, < — /i2 < 00: 



E. 



surf 
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h. 
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dx- 



-{-^-l)x (-r--l)x 



+ e- 



dx- 



coshx 
(— 



+ e- 



coshx 



(6.36) 



where the second line was obtained from the first one after a simple manipulation. This 
agrees with the results of For hi = h2 = one has from (|6.36| ) Esurf = (vr — l)/2 — In 2. 



6.7 The case of parallel magnetic fields 

Another aspect is the structure of the ground state in the regime hi^2 > 0. Assuming that, 
for example, for hi^2 > hf^ the ground state contains both left and right boundary 1-strings 
to minimize the energy, we end up after a short calculation with a half-integer spin of the 
vacuum, which signals that such a state cannot, in fact, be the vacuum. Hence, the ground 
state must have a more intricate structure. Appealing to the Ising limit 7 00, one sees 
that for hi 2 > A the ground state must have both boundary spins directed opposite to 
the magnetic field, and therefore contain a kink in the bulk (recall that L is even). On 
the other hand, for /ii 2 < A the lowest energy configuration is such that the boundary 
spins are antiparallel, which means that the physical vacuum contains what was called a 
boundary excitation at one of the ends. This suggests that for finite A the correct ground 
state wave-function of the Hamiltonian ( |6.1|) should contain a bulk hole with the minimal 
possible energy (i.e. the kink with zero rapidity ^ = 0) and both boundary 1-strings when 
^1,2 > . Such a state has spin zero. Changing the rapidity of this stationary kink away 
from zero, one obtains in such a way an excited state whose energy can be arbitrarily close 
to the vacuum one, which means that there is a new gapless branch in the spectrum. Q 
Similarly, when h^J^^ < hi^2 < ^cr^ ^^e ground state to have the integer charge it should 
also contain a kink in the bulk. When /ii 2 < h^J^^ the physical vacuum contains only one of 
the two boundary 1-strings and no stationary kink in the bulk (when hi = h2 there are two 
possibilities to have either left or right boundary 1-string in the vacuum, corresponding to 

^In the Ising limit 7 ^ cx3 the energy of the kink is independent of rapidity, and therefore this branch 
degenerates to the vacuum. 
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the obvious two-fold degeneracy of the Ising ground state in this case). Such a state has a 
smaller energy for /ii 2 < h^^r than the one with a hole in the bulk and two boundary strings, 
while for /ii 2 > h'^^) the state with the bulk hole is energetically preferable, since in this case 
Eh > Sh (see Fig. 12 and |^). This situation is in some sense analogous to the case of the 
periodic antiferromagnetic XX Z chain with odd L, where the ground state contains a kink. 
According to the above discussion the surface energy in the case hi^2 > hum is: 



1 / 00 ^-yHin _|_ g- 

Esurf = 7^ihi + h2) - giA) + EhiO) - sinh7 1 + ^ r 

2 \^ cosh7?T, 



(6.37) 



In the rational (7 0) limit eh{0) vanishes and eq. (|6.37]) becomes 



E, 



surf 



n 
2 



,-(7^7-1)- 



dx- 



+ e- 



coshx 



(6.38) 



This expression agrees with the one obtained in |^2|. Note that the authors of |^2| obtained 
eq. ( |6.38| ) under the assumption that < hi^2 < her, while our derivation shows that this 
result is valid for < hi^2 < cxo. In the Ising limit eq. (|6.37|) gives the correct result Egurf = 
(3A — hi — h2)/2. Observe that for the XXX chain the following equality holds (see ( |6.36 ) 



and (|6.38|) ): i?s«r/(^i5 ^2) = -Ssur/(^i5 — ^2)- This is the consequence of the decomposition 
Esurf = f{hi) + f{h2) +const, which takes place in the limit L = 00 when two boundaries are 
independent, and the obvious property of the semi-infinite chain f{—h) = f{h). The same 
statements are true for the surface energy of XXZ chain apart from the 6h{0) contribution 
(see (p3)). 



6.8 Remarks 



We would like to mention also that within the BA technique it is possible to calculate also the 
boundary S-matrix for the scattering of kinks (represented by holes in the Dirac sea) in the 
ground state of the Hamiltonian ( |6.1|) or in the excited boundary state. Such a calculation 
has been performed in H^] for the boundary XXX chain and in chapter 4 above for the 
boundary sine-Gordon model. For the boundary XXZ chain these S-matrices have been 
obtained by Jimbo et al. 



49[| by the algebraic approach. 
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h 



Figure 6.1: Solid line: a schematic plot of the energy of the boundary excitation, £b{h), 
as a function of the boundary magnetic field h. Shaded area: the energy band of the bulk 
excitations. 
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Chapter 7 



Calculation of correlation functions 
for the problems with impurities 



We show how to compute analytically time and space dependent correlations in one-dimensional 
quantum integrable systems with an impurity. Our approach is based on a description of 
these systems in terms of massless scattering of quasiparticles [^. Correlators follow then 



from matrix elements of local operators between multiparticle states - the massless form- 
factors. Although, in general an infinite sum of these form-factors has to be considered, we 
find that for the current, spin and energy operators only a few (two or three) are neces- 
sary to obtain an accuracy of more than 1%. Our results are valid for arbitrary impurity 
strength, in contrary to the perturbative expansions in the coupling constants. As an exam- 
ple, we compute the frequency dependent condunctance, at zero temperature, in a Luttinger 
liquid with an impurity, and also discuss the succeptibility in the Kondo model and the 
time-dependent properties of the two-state problem with dissipation. This chapter is based 
on 



7.1 Introduction 

In this chapter we present the technique to calculate current-current correlation functions in 
the quantum field theories of the form 

1 PCO 1 

H=- dx[87igU' + —{d,^f]+B, (7.1) 
Z Jo ong 

where S is a problem dependent boundary interaction and the fields are defined on the 
positive half-line. The method we adopt is based on the form-factor formalism |]71| and 
provides us with the series expansions of correlators. One has to insert intermediate n-particle 
states in the correlator (0|j^jy|0) = J2n{^\jiJ.\^){''^\ju\0)] then the well-known technique for 
integrable models [^] gives the exact values for the above scalar products between the 
vacuum and n-particle states (the form-factors). The integral representation for each term 
of the series is available. Moreover, the observed rapid convergence at any scale from 
small to large distances allows us to truncate the series after a few terms (typically two 
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or three) to obtain a 1% accuracy. There are few novelties that deserve our attention: first, 
we have to deal with the massless form- factors, following the pioneering work . Second, 
we work in a half-plane geometry. An instructive example, although somewhat simple and 
abstract, is an Ising model in a boundary magnetic field ||7i| . 

The role, played by integrability in our approach is two-folded. First, it gives us the exact 
expressions for the form-factors and other necessary quantities. Second, it provides us with 
a basis of particles (intermediate states) which allow us to truncate the series expansion of 
correlators due to the rapid convergence. Thus, the integrability insures the relatively small 
contribution of the multi-particle terms. 

To pave our path through the calculations, we consider some simple model of the free 
massless 1-dimensional scalar field with the boundary interaction of the form B = Mb cosh (p. 
This model can be viewed as the massless limit of the integrable sinh-Gordon model with 
boundary interaction: 



S = dxdt 



-Mb dtcosh[^]U=o (7.2) 



1 T77^ 

{d,<j>r - — cosh{g<j>) 

2 

The advantage of this approach, as opposed to starting with the free massless boson quantized 
in the plane- wave basis, is the following: in the massless limit, obtained by scaling the 
energy and momentum of particles along with the boundary mass Mb to infinity, we get the 
convenient basis of massless particle states which are particular combinations of plane waves 
that scatter diagonally off the boundary. Corresponding classical solutions are presented in 
3T[| . Working with these massless particles at first sight adds some complexity, but it is paid 



off by the final simple and manageable results, while in the plane wave picture one has to 
do infinite summations and the final result is difficult to extract. 

We apply our technique to find the correlation functions in three models of condensed 
matter physics: the Kondo model, the spin-boson model (two-level system with dissipa- 
tion), and the Luttinger liquid with impurity (realized as the quantum Hall liquid with 
constriction). A common feature to all these models is that they can be reduced to a model 
described by massless excitations in the bulk interacting with an impurity at the boundary. 
The boundary interaction B in the first and second models is B = \ (S+e^'^^^^^'^ + S^e^^'^^^^^'^), 
and the third model has B = M cos 0(0). The absence of a mass gap leads to a power law 
behaviour for the current correlators in both the ultra-violet and the infra-red regime. The 
cross-over between these two regimes is non-trivial because of the renormalisation group fiow 
induced by the impurity. For each model that we study, correlation functions can be related 
to the measurable quantities in the model dependent way. For example, for the Luttinger 
liquid, the Kubo formula gives the AC conductivity: 



CO roo 



a^,{uj)=J^ dxj^ {0\UO)j^{t + ix)\0)e-''^'dt (7.3) 



The perturbative analysis of shows that the conductivity scaling function, depending 
on the impurity strength, fiows from the insulating at T = to the perfect conductance 
at T = cx) in the repulsive regime. The alternative exact methods, which do not employ 



the knowledge of correlation functions, have been developed in fT^ to find the static DC 
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conductivity at non-zero temperature and voltage. These methods do not allow, however, the 
determination of the AC conductivity. The experimental data as well as the numerical 
Monte-Carlo simulations |79|, Q for the Luttinger liquids with impurity are available. 



7.2 Integrable models in condensed matter physics 

In this section we introduce and review three prominent models of condensed matter physics: 
the Kondo model, the spin-boson model (two- level system with dissipation), and the Lut- 
tinger liquid with impurity. The results of application of our technique for computing cor- 
relators to these models are partially shown here and in part further. When bosonized, all 
three models look alike in the bulk (free massless boson), but have different boundary inter- 
actions B: A (S'+e*'^(°)/2 _^ S'_e-^<^W/2) for the first and second, and Mcos(/)(0) for the third 
one. 

7.2.1 Kondo model 

The Kondo model describes the free bulk electrons interacting with an impurity via the spin: 

Hk = J2 <k)ctcka + JmSrmp (7.4) 
k,(j 

Here, Simp = is the impurity spin, a are Pauli matrices, and s(0) is the spin induced by 
the electrons at the point of impurity, 

^"(0) = ^E^^(0)^-'^-'(0)' V^-(O) = L-^'^Y.Cka. (7.5) 

era' k 

A useful generalization of (|7.4|) is an anisotropic Kondo model ||81|| : 

Hk = J2^(^)4aCka + J±[(T+ij^1p'^{0) + a-lp^ tp l{0)] + J,a, ^4aCk'a, (7.6) 

k,a fc,fe',(7=± 

which reduces to ( |7.4|) if = J±. 

We need also the bosonized version of (|7.6| ). For it, two bosonic fields are necessary: one 
associated with charge and one with spin. The charge- density field decouples completely and 
only the spin-density has interaction at the boundary. The Hamiltonian for the spin-density 
field is of the form 

1 1 

Hk = - dx[87igU^ + — (9,.0)2] + A (S+e'-^^o^/^ ^ ^_e-^^(o)/2)_ (7,7) 

2 JO Sng 

The coupling constant A is related to J±. The Sz term in the Hamiltonian has been eliminated 
by a unitary rotation, but dependence on is implicit in g. The case g = 1 corresponds to 
the isotropic Kondo model. 

The quantities of interest in the Kondo model are the spin-spin correlation functions for 
the induced electron spin in the bulk (screening cloud problem |8^) and for the impurity 
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spin on the boundary, (>S'f^p(t)S'f^p(0)). We consider here mostly the latter ones, since they 
are relevant for the dissipative two-level problem. But the screening cloud problem can be 
approached by the form-factors technique, too, and for partial results we refer to |]70|. The 
dynamic susceptibility, or the response function is given by 

x"M = i/rfte^-*([5^(t),5^(0)]). (7.8) 



In it was shown how to express the spin correlators in terms of the bosonic field correlators 
of (|7.7|). In particular, for the response function we haveQ 

1 1 



Im [Q{—iuj, 13b) — Q{—iuj, —oo) 



where 



/■oo 

{d-,ct>{x, y')d,<P{x, y"))x = / dEg{E, 13b) exp [-2Ex - iE{y' - y")] , 

JO 



(7.9) 



(7.10) 



and the boundary temperature 13b is related to A. The foregoing technique (see section 7.4) 
allows us to compute the quantity Q{E,f3B) in the form of series expansion. Every term 
of the series is known analytically. In practice, it is enough to keep only few first terms of 
these series to obtain a very good accuracy of 1% or more. E.g., for the Toulouse limit, 
corresponding to g = |, the series get truncated and we obtain the exact result 



x"i^) 



2 Tb^ 
— — Im 



X + iTB){uj - X + iTB) 



ATI 




1 -^ ^ 
+ -z— tan —— 

Tb Tb 



For the less simple g 

6x"iuj)^'^ 



7r2 cj2 + 

1/3 case the leading contribution to the response function reads 



2 j2 



■Tie 



tanh lo^ 



tanh(|log(^; 



+ 



(7.11) 



where n ~ 3.14, d ^ 0.1414. Similar computations give rise to the results in figure 7.1 where 
we plotted x"{^)/^ for the values g = 3/5, 1/2, 1/3, 1/4. Notice the emergence of the peak 
for g < 1/3. It is remarkable that this peak appears at g = 1/3 and not at = 1/2 as was 
expected from other means of calculations. 

Using eq. ( [7.9| ) and the general form of Q{E,Pb) given below, we can find the exact 
static succeptibility Xo which is given by the Shiba's relation |8F 

X"(^) 



lim 



UJ 



Omitting the details, we list the result: 



Xo 



(7.12) 



(7.13) 



B 



(but see ||7^ for the complete account of this derivation). 
^In the most of this chapter it is tacitly assumed that h = 1. 



91 




92 



7.2.2 Quantum systems with dissipation 

The effect of dissipation on quantum tunnelling has been addressed in the pioneering work 
of Caldeira and Leggett Suppose first one has a system with friction, or dissipation, 

described by the equation of motion 

dV 

Mq + r]q + — = F,,t{t), (7.14) 
oq 

with the potential V admitting a quasi- stationary state that can tunnel into the continuum. 
Then, speaking loosely, the question one asks is whether there is a qualitative change in the 
behavior of such a system for 77 = and 1] 07 

It is an interesting problem in general how to introduce friction in the formalism of 
quantum mechanics. In the following phenomenological model is suggested, in analogy 
with ([7l4l ): 

C<p + R-'^+ — = (7.15) 

It is usually refered as RSJ (resistively shunted junction) model. The nature of the resistance 
Rn is not clear, however. It is introduced in the model phenomenologically. 



As the analysis of shows, the net effect of the friction is that the WKB exponent for 



the tunnelling rate is increased by dissipation and tunnelling rate decreases: 

P ~ exp f ) (7.16) 



huj, 



eff, 



As a consequence of interactions, the characteristic frequency of classical system uq (defined 
by the product LC) is renormahzed to uj^ff- 

As for experimental applications, the theory of MM describes a trapped flux in SQUID 



(superconducting interference device). The phase of the flux, 0, is a quantum degree of 
freedom governed by the potential 

u{<i>) = ^ '^*°'°H^J' ^ ^ ^ 

where (px is an external flux, L is a self- inductance of the device, Ic is a critical current and 
"^•0 = /i/2e is a flux quantum. The kinetic energy is given by = |C0^ with C a capacitance. 
Another experimental realization can be found in the single Josephson junction biased by 
a fixed external current Iq. The quantum degree of freedom now is the phase difference of 
Cooper pair across the junction, ip, and the potential is 

U = — — y?-— — cosy?. (7.18) 
Zn zn 



Another fundamental system has been studied in [^]. Here, instead of tunnelling into 
the continuum, a double well problem is being considered with the tunnelling between two 
separated wells. Since the particle can be localized in either of two wells, the system is refered 
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to as a two-level system. More generally, it could be not necessarily a continuum degree of 
freedom with two spacially separated states, but also an isospin degree of freedom with two 
states, e.g. a strangeness quantum number in particle physics. Examples and applications 
include the strangeness oscillations of a neutral K-meson beam, or an inversion resonance of 
the NH3 molecules. 

The isolated two-level system can be modeled by a Hamiltonian 

1 1 

Ho = --Aoa^ + -ea, (7.19) 

where Aq is a matrix element for tunneling, and £ is a "detuning" parameter called bias (e.g. 
in the double- well problem it is the difference of the ground state energies of the wells). It 
is obvious that the model (|7.19|) is the effective Hamiltonian in the quasiclassical description 



of the double-well problem. 

The dynamics of the two-level system with the Hamiltonian ( |7.19|) can be easily solved. 



It is well-known that in the absense of bias, e = 0, system shows coherent behavior with the 
probability distribution 

P{t)=PR-PL = cos{Aot), (7.20) 

where Pl {Pr) is the probabilty of finding particle in the left (right) well, and we assumed 
that initially the particle was localized in the right well. In the presence of non-zero bias the 
oscillatory behavior ( [7.20| ) is destroyed. 

Another factor that can potentially destroy the coherence is an interaction with the 
environment with sufficiently strong coupling. Environment couples through the term in 
( [7. 191 ) and can be modeled the by following spin-boson Hamiltonian 



The environment is modeled by the set of N harmonic oscillators (phonons) described by the 
last term in ( [7.21[ ), and, of course, we are interested in the thermodynamic limit 00. 
With the non-zero coupling to environment, go? one state of the two-state system becomes 
more preferable than the other, depending on the "mood" of the environment. At T = 0, 
classically, all the oscillators are at = and have no effect on the two-state system. 
However, quantum-mechanically the problem becomes non-trivial even at T = due to the 
quantum fluctuations. The question of theoretical interest can be formulated as follows: Can 
there exist oscillatory behavior (|7.20|) in a macroscopic system (|7.21|) , or the coherence will 



be destroyed by the interaction with environment? 

The time-dependent quantities of interest which are useful in the analysis of this problem 
are conveniently encoded in the following: 

P{t) = (a,), (7.22) 
C{t) = \{[aM.<y^m)■ (7.23) 
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Here, P{t) measures how the average of spin varies with time, provided that at t = it 
was in a certain state, while C(t) describes the probabihty to be in a state (Tzit) given that 
the system was in a state a^iO) at t = 0. The operators are understood in the Heisenberg 
representation az{t) = e-'^^shtfj^^iH^it ^ 

It turns out that the effect of phonon bath on the two-level system is very non-trivial 
and depends on the form of the spectral function of the environment, 

^M = fE;^^(^-^0 (7-24) 

It is usually assumed that J{uj) = rju^ exp(— cu/c^c) where Uc is a cut-off frequency. According 
to the analysis of |]85[, we have, at T = 0, the weakly damped coherent oscillations for s > 1, 



and the complete localization for s < 1. For s = 1 the analysis becomes more complicated 
and the result depends on the value of the dimensionless coupling constant a, 

a = ft (7.25) 

The s = 1 case is usually refered to as the ohmic regime. For a > 1 and s = 1 it is believed 
that the system is completely localized while for a < 1 the situation is not clear yet. It 
is known that there exists a critical point at the value of a equal or about etc = | where a 
phase transition in the ground state occurs. For a less than this critical value the damped 
oscillations are observed, while above it there is an incoherent relaxation. 

The value a = |, called the Toulouse limit, is an exactly solvable point. At this point 
the model can be mapped onto the following model: 

Ht = J2 ^4ck + VY.{d+Ck + ctd) (7.26) 

k k 

where d^,d create and annihilate a localized state (corresponding to the spin degree of 
freedom in the spin-boson model), while c^, Ck are the creation and annihilation operators for 
the fermions in continuum (corresponding to the bath).0 The Toulouse model is a particular 



case, f/ = 0, of a more general exactly solvable resonance-level model |8^, which can be 
mapped onto the spin-boson model : 

Hrl = Ht + Uj2ictck' - l)id^d - 1) (7.27) 

k,k' 

The correspondence is given by the simple relations 

(7+ = d+, a_ =d, az = d+d - ^ (7.28) 



and V is directly proportional to the hopping Aq, while U is related to a as follows P? 



(7.29) 



TT 



The bias e corresponds to the energy of locaUzed state. 
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The exact solution of the Toulouse model gives the following behavior at a = |: 

P{t) = exp(-t/r). (7.30) 

A particularly useful fact for us is that the long-time behavior, t ^ uj^^, of the spin- 
boson model in the ohmic regime with ^ < a < 1 is the same as the long-time behavior of 
the anisotropic Kondo model f\ ( |7.6| ). This result can be established rigorously by means of 



the bosonization HS^j- The harmonic oscillators in the spin-boson model play the role of the 
spin-density excitations in the Kondo case, and the impurity spin {Simp = |) corresponds to 



the discrete degree of freedom in ( [7.21|) . We have 



a ~ (7.31) 

The point a = 1 maps into = 0, in such a way separating the ferro and antiferromagnetic 
regimes in the Kondo model. Note that the mapping works only for J±, Jz small enough, i.e. 
a close to 1. The regions where a is far from 1 can be approached by the renormalization 
group analysis . By continuity, the results should hold at least up to the strong coupling 
fixed point, a = |, where the exact solution ( [7. 301) holds. Note that the analysis of the 



previous section, in particular figure 7.X suggests that the incoherent relaxation takes place 
in fact up to a = 1/3. This is supported by a recent RG numerical study |]89[| . 

7.2.3 Quantum Hall liquid with constriction 

Transport in one- dimensional interacting electron systems in the presence of impurities is an 
instructive problem with a broad rangle of experimental applications. Important ingredients, 
defining a proper physical model, are the Coulomb correlations in the vicinity of a tunneling 
point, as well as the electron interactions in the "feeding leads". Namely, when the leads 
are one-dimensional, they ought to be described by a Luttinger liquid, rather than the 
Fermi liquid since the latter is de-stabilized by interactions. In real ID wires the impurities 
away from the point contact will complicate matters, tending to localize electrons. This 
localization makes it difficult to realize ID Luttinger liquids in experiment. 

However, systems are available J75| , [T^ which are free of undesirable localization - where 
the leads feeding a point contact are 2D fractional quantum Hall (QH) fluids. In the ap- 
propriate regime of the QH liquid the incoming current will be carried by edge states. Wen 
has demonstrated that the gapless edge excitations of a QH system are chiral Luttinger 



liquids. Let us recall briefly the logic behind Wen's theory. The long-length-scale physics of 
the bulk 2D QH state is described by the massive 2+1 dimensional Chern-Simons theory 

Sbuik = T— / CLf.d^axS^uxd'^xdT. (7.32) 



'Only the low-lying excitations are relevant in this limit. 
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The coefficient v is uniquely fixed by the quantized Hall conductivity. ^ In the presence of 
a boundary, say at ?/ = 0, one can integrate out the bulk degrees of freedom. The resulting 
ID action for the edge field reads 

Sfs.ge = / dxdr{d^<^R){idr + vd^)^R, (7.33) 

where v is the velocity of edge excitation (we set f = 1 hereafter). The boson field was 
introduced as aj = djip to solve an incompressibility constraint on the electron density, 
eijdittj = 0. Similarly, one writes for another edge 

Sedge = J dxdr^d^Lp l) {-idr + vd^)ipL- (7.34) 

The charge density along the edge is given by p{x) = dx^p/27i, and the momentum operator 
conjugate to is 11 = p/u. Adding an extra electron to the edge is equivalent to creating a 
soliton in ip with electron creation operator being 

^(x) ~ exppvri J"" U{x')dx'] = e^^^^^/^ (7.35) 

A quasiparticle of fractional charge ue is created by e*'''*^^^ (speaking loosely, an electron is 
composed of z/~^ quasiparticles). 

Due to the chirality, backscattering is only possible when opposite edges of the sample 
are close together, i.e. at the point contact. Thus, localization in such leads is absent. The 
analog of impurity that causes backscattering is a narrow constriction which brings left and 
right edges close enough for Laughlin quasiparticles to tunnel, as illustrated in figure 7.2. 
This is achieved by applying a gate voltage Vg across the narrow region in the Hall bar. 

Suppose that tunnehng takes place at x = 0. The edges are no longer independent, but 
rather coupled via the tunneling term in the total action: 

S = Sjfdgei^L) + S^^g^i^n) + j drVimpi^L, ^r)x=o- (7.36) 
The most general form of the tunneling term is 

oo 

V = ^^^i?(0) + h.c. = t;„e''"['^^(°)-'^«(°)l + c.c. (7.37) 



m=l 



where Vm are tunneling amplitues. The term with m = 1 corresponds to quasiparticle 
tunneling, while m = term - to the electron tunneling. In what follows we choose 
u = 1/3. As was argued in [|7^, for this value of u only the term m = 1 in ( |7.37] ) is 



relevant, corresponding to the transfer of the Q = e/3 Laughlin quasiparticle. For generic 
filling fraction, many types of quasiparticles contribute as relevant charge transfer. So, the 



*For simplicity, we choose ly ^ to be an odd integer to have only one branch of edge state. 
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Figure 7.2: Quantum Hall bar with tunneling between the edge states. 

model of interest has the following Hamiltonian (in a slightly different notations for further 
convinience) : 

1 1 

H=- dx[8nun^ + ^{d^^f] + X5ix) cos{ipL - ^r), (7.38) 
Z J-oo ovrz/ 

where the L and R components depend on x,t a.s ipL{x + t),(pR{x — t). One can tune the 
tunneling amplitude A by varying the gate voltage Vq- The Hall voltage V between the two 
edges of the QH liquid can be introduced in the model by adding the term —e{QL — QR)V/2h 
to the Hamiltonian. This substitution has an effect of placing the charge carriers injected 
into the left and right edges at different chemical potentials. As can be easily seen, by shifting 
the fields (Pl,r we can bring the voltage dependence in the boundary term: cos{ipL — fR+^^t). 
We will deal with the = and T = case only. 

As discussed in ||7^, in order to map (|7.38|) to a boundary problem on a half-line, it is 



convenient to proceed in two steps. First introduce : 

0"(x + t) = -^[fL{x,t)+yDR{-X,t)] 

r{x + t) = -^[^L{x,t)-^n{-x,t)] (7.39) 

which are both left moving. The even and odd charges are related to the charges of the 
original left- and right-moving edges by 

AQ = Ql-Qr = V2Q°, Ql + Qr = V2Q'. (7.40) 

The backscattering current is related to Q°, whereas the total charge on both edges, Q'^, is 
conserved even in the presence of impurity. It is clear that the interaction term does not 



98 



affect the even field, which therefore remains free. As for the odd term, it can be mapped 
onto a boundary problem as follows. Define : 



x,t) = \/2(f)°{x + t), X > 0, 



(pR{x,t) = V2(j)°{-x + t), x>0. (7.41) 
The odd hamiltonian then reads : 

1 roo 1 1 

H=-l Mn-'f + — (5.0°)2] + A5(a;)cos-</)°, (7.42) 

and in the following we will write = 0° and g instead of u. Thus, for this problem, 
B = X cos |0(x = 0, t). 

The quantity of interest in this case is the AC conductance at vanishing temperature. A 
standard way of representing it is through the Kubo formula ( [7.3|) . It is usually easiest to 
calculate ( |7.3| ) using Matsubara formalism |^T|. First, one defines current- current correlator 
in the Matrubara formalism : 

'i 

UJ 

n„^(x,r) = -{T,Ux,T)jf,ix,0)) (7.43) 

n„^(x,zcuM) = f dTe'^^'^Il^p{x,T) 
Jo 

Then, one substitutes iuM ^ uj + i6 into Uap and sends 6—^0. Eventually, we get for the 
real part of conductivity tensor : 

Rea«/3 = --ImU^piu). (7.44) 

UJ 

The imaginary part of the spectral function in the latter formula can be expanded as follows: 
- ImU^pioo) = 7r(l - e-'^")e'^^ ^ e-^^"(n|j„|m)(m|j^|n)<5(u; + E„ - E^). (7.45) 

n,m 

In the limit of zero temperature, /3 —>■ oo, the sum over n contains only one term - the 
ground state, so one has : 

-lmU^f,{u) = 7iY.{0\ja\m){m\jp\0)S{uj - Em). (7.46) 

m 

In the above formulas ja,/3 = j is the physical current in the unfolded system, j = dt{(pL—'^R)- 
Without impurity, the AC conductance of the Luttinger liquid is frequency independent, 
G = g. When adding the impurity, it becomes G = |+AG. After some simple manipulations 
using the folding, one finds : 



AG{ujm) = / dxdx' / dye 



[(a,0(x, y)d,4{x', 0)) + {d,(j){x, y)dA{x\ 0))] , (7.47) 



where z = x + iy. 
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7.3 Boundary sinh-Gordon model 



We discuss in this section one other boundary integrable model, the boundary sinh-Gordon 
model with B = M cosh 0(0). It is related to the boundary sine-Gordon model by the analytic 
continuation in the coupling constant g, but we do not make use of this fact here. Instead, 
we use it as a toy model for studying the correlation functions. It is fairly straighforward 
then how to tackle more complicated models. The advantage of the sinh-Gordon model is its 
technical simplicity: the particle spectrum of this model consists of only one scalar particle 
- a sinh-Gordon boson. The bulk sinh-Gordon model is well-studied in the literature. 



7.3.1 The boundary reflection coefficient 

We need to know the boundary reflection coefficient for sinh-Gordon model: 

Mb f , fg(l){x = oy 



SshG 



dxdt 



l{df,4>Y - ^ cosh{g4>) 
2 g^ 



9' 



dt cosh 



(7.48) 



For this we use the fact that the action (|7.48|) is related to that of sine-Gordon model 



5, 



SG 



dxdt 



Tfl 



+ 



B 



(3' 



dt cos 



'i3(t){x = oy 



(7.49) 



by the analytic continuation in coupling constant g = i(3. Recall that the sinh-Gordon model 
has only one particle with neutral charge. A useful observation is that the bulk scattering 
matrix for this particle. 



S{d,B) 



tanh 



tB- 



tanh |( 



+ zf)' 



B 



2g^ 



171 + g' 



(7.50) 



can be obtained from the scattering matrix of the lightest breather of the sine-Gordon model 
1^ by the mentioned above analytic continuation in coupling constant. We assume that the 
same is true for the boundary reflection matrices. Partial conflrmation of this can be found 
in |9^, where this statement was proven in the semi-classical limit. So, we will use the 



result obtained in for the reflection coefficient of the lightest breather in the boundary 
sine-Gordon model. The general solution of Ghoshal has two free boundary parameters pl| : 



Rmass{0) =U{r^,d)U{i§,9)R,{e) 



(7.51) 



One can argue that for the boundary term of the form Mcos(/50/2) one should set ?7 = 0, 
the other parameter d being related to Mb- So, the reflection coefficient for the boundary 
sinh-Gordon model with boundary interaction cosh(5f0/2) obtained from (|7.51| ) by the 
analytic continuation reads: 



Rmass (^) 



tanh i| 



+ 



2 



tanh ■ 



- f ) cosh \ 

2 2 



) cosh I ( 



- HLB)coshi( 



I cosh 



1 , 



— y) cosh h 



\ _j_ ittB ' 



cosh i ( 



inB iiT \ 

^ ^ (7-52) 

' 4 2 
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The massless limit of ( [7.52D is obtained by sending 9 and d to infinity. We have in the 
massles limit 

1 ?7r 

R{e) = -i^nh-{e-eB--) (7.53) 



where 6 is now massless rapidity related to the momentum as in (|7.64|) . This is as well the 
massless limit of the sine-Gordon reflection matrix. The boundary parameter 9b is related 
to the boundary mass Mb in the way unknown to us. Note that in terms of momentum 
expression (|7.53|) looks like 

R{p) = —— 

P + IPB 

The definition of K{6) in the massive case 

-^massi,^^ Rmass ( "2" 

becomes in the massless limit 

K{e) = r{^^ + 0) (7.54) 

This follows from 

ZTT 

K{e) = lim K^assiOo + 9)= lim Rmassi— -eo-e) = 

= lim s{2e)Rmass{^ + eo + e) = r{- + e) 

where we used massive crossing-unitarity relation |I| above. 
7.3.2 Sinh-Gordon form-factors 

The form-facotrs for the massive Sinh-Gordon model were found in [Q. Let us list here 
some of them that we need. For the field itself, the form-factors between the ground state 
and n-particle states are: 



2n+l 



(/?l,...,/?2n+l) = (O|0(O)|/5i,...,/32„+l) 



i^2.+iQ2n+i(e^\ e^--) n J'^ (7-55) 



where H„ are normalization constants 



4sin^ 



V2_ 

Q2n+i are symmetric polynomials in the variables Pi = e'^*, 

Qi = l, Q3 = exp(/?i+/?2 + /33), 
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Figure 7.3: Plot of \Fmin{P, B)]"^ for two different couplings. 



(5)^^(5)^(5) 2 (5) 



ci = 2 cos 

2 



By parity, only form-factors with odd number of particles are allowed. We used the standard 
basis of elementary symmetric functions af"'' in variables pj in the formulas above, and 



F 



D exp 
D = exp 



(ixsinh^sinh(l - y)|sinh| . ^x{i'K - (3) 



sm 



X 



sinh X 

dx sinh sinh(l - f )| sinh | 



27r 



a; 



sinh X 



The plot of -Fmm is given in Figure 7.3. 

Since the sinh-Gordon form-factors ( |7.55| ) for the field operator are scalars with respect 
to the Lorentz transformations, they are invariant under the simultaneous shift of all the 
rapidities (3i — > (3i + (3q, where (3q is some constant. Taking the massless limit implies making 
the substitution (3i = 6i + (3q with (3^ ^ oo for the right-moving particles, and making the 
substitution [3i = —Oi — (3q for the left-moving particles. Here 6i are the rapidities of massless 
particles. So, the massless form-factors derived by this procedure are: 



(010^(0)1^1,..., 



(O|0i(O)|^l,...,e2n+l)j 



where we used the fact that F„ 



F„ 



( Gil} Gl f 



g&i _|_ g&j 



i{9i — 9j) 



i<3 



-6) = Fmin{(^)- The normalization factor Nq can be fixed 
by comparison with the massless free boson two-point correlation function (|7.57|) . Inserting 
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into the vacuum two-point correlation function of free fields the full set of normalized in- 
termediate states, and using the form-factors above, we obtain the contribution from each 
n-particle state in the form 

with some constant coefficients a„, 

= / %.yT' + ■ ■ ■ + e»-)%-<-"--'")|F„(9,. . . . .9„)P. 

The normalization factor A^^q can then be obtained by summing the series J2 o-n and requiring 
the result to be 1. One can compute a„ numerically and observe that they decrease very fast 
with n, and so can be truncated after a few terms. In particular, using one, three and five 
particle form- factors we found the approximate value of Nq: 

No = 1.005 B = 1.0 

No = 1.0002 B = 0.1 

Note that the rate of convergence of the series (^n gives us a hint about the convergence 
of analogous series in the case with the interaction at the boundary (see below). The rapid 
convergence holds in the massive case as well, and has a physical explanation P^ . 

7.4 Calculation of correlation functions 

In the massless limit the Green's function of current-current type on the half-line can be 
factorized as follows: 

G{xi,U,g,MB) = (O|9,0(xi,ti)9,0(x2,t2)|O) = 

= Go{xi- X2,ti-t2,g) + Gi{xi+ X2,ti-t2,g,MB) (7.56) 

(this can be easily seen from the form-factor approach below). The first term Go does not 
depend on boundary coupling and is in fact equal to the current-current correlation function 
in the translationally invariant system on the line. It is the second term Gi where the 
breakdown of translational invariance manifests itself explicitly and which carries all the 
dependence on the boundary coupling. 

In general Mb is the only dimensional coupling that enters the Green's function. There 
are two values of Mb where the theory is scale-invariant. One of them, Mb = 0, corresponds 
to the free boundary condition dx(t>\x=o = 0, while another one, Mb = oo, corresponds to 
the fixed boundary condition (j)\x=o = 0. For these values of Mb the Green's functions are 
known exactly: 

G\mb=o = Goixi - X2) + Go{xi + X2) 
^1^-/3=00 = ^0(3^1 — X2) — Go{xi + X2) 
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where ^ ^ 

Go(x,t) = ^^^^ (7-57) 

Thus, the boundary mass Mb induces the boundary RG flow from short to large distances, 
which vary Gi between —Gq and +Go- 

We will work on the half-plane which geometry is shown in Figure 7.4. Based on the 
euclidean duality, there are two alternative ways to introduce the Hamiltonian picture [0. 
First, one can take x to be euclidean time. In this case the equal time section is an infinite 
line X =const, y G (—00,00). Hence the associated space of states is the same as in the 
bulk theory. The boundary at a; = appears as the "time boundary," or initial condition at 
X = which is described by the boundary state \B) (a particular state from the bulk Hilbert 
space). The correlators are expressed as 

where Oi{x,y) are the Heisenberg local field operators 

O,(x,i/)=e-^'^O,(0,y)e^^ (7.59) 

and %; means x-ordering. 

Alternatively, one can take the direction along the boundary to be the time. In this 
case boundary appears as a boundary in space, and the Hilbert space of states is associated 
with the semi-infinite line y =const, x G [0, 00). The correlation functions of any local fields 
Oi{x,y) can be computed in this picture as the matrix elements 



/n t \ n ( \\ ^-^\Ty[Ol{.xuyl)...ON{.XN,yN)W) 
{Oi[xi,yi)...ON[xN,yN)) = j^^^ , (7.60) 

where ||0) is the ground state of the boundary Hamiltonian, Oi{x,y) are understood as the 
corresponding Heisenberg operators 

0,{x,y) = e-^^^O,(x,0)e^^^ (7.61) 

and Ty means y-ordering. 

The equality of expressions ( [7.58|) and (|7.60|) can be understood as a definition of the 



boundary state, which is chosen such as to provide the equivalence of correlators. 



7.4.1 Boundary-in-time representation 

First suppose that a; > is imaginary time coordinate and y is space coordinate, z = x + iy. 
It means that time translation is performed by the operator T = exp{xH), with H being 
the bulk Hamiltonian. 

We wish to compute the following matrix element: 

(s|a,0(xi,2/i)9,0(x2, 1/2)10) = {B\d,Mzi)d-zM^2m 

+ {B\d,Mz-,)dMz2m (7.62) 
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Figure 7.4: The geometry of space-time. 



where (/> is a massless free field 
state of the sinh-Gordon model 



^l{z) + (pR^z) and \B) denotes the massless boundary 



oo „ 

\B) = |o) + E / 



^_-|^ J— oo<6*i<...<9jv<oo 



.1=1 



AL{9M)...AL{9^)An{9,)...AR{9Mm (7.63) 



Because \B) has chirality zero, products of the fields of the same chirality do not contribute 
to the right hand side of eq. ( [7.62| ). Parameter ^ is a massless rapidity in terms of which 
the energy and momentum of particles read: ^ 



Er = -Pr 



lie 



E 



R 



Pr = /^e^ 



(7.64) 



and Al ji{9) are the left and right moving particle creation operators, the particles in our 
case being the massless sinh-Gordon bosons. The fact that the sinh-Gordon particle bounces 
elastically off the boundary allows us to have a much simpler expression for the boundary 
state ( [7. 631) than that for the plane waves. Substituting the boundary state ( [7.63|) into ( [7.62|) 
and using the fact that left (right) moving field acts only on left (right) moving particles, we 
obtain the following expansion in terms of the form-factors: 



B\d^(j){xi,yi)d^(f){x2,y2)\0) = (O|9^0(xi, yi)9^0(x2, ?/2)|0) 



J— oo<Si<...<9jv<oo ^71" 
2 

-{xi+X2+iyi-iy2 



X 



' N 



(010^(0)1^1, ...,e,v)L(O|0,j(O)|ei,...,^^^) 

HEL^n + c.c.) 



R 



^Here /x is an arbitrary mass scale. 
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Go{Zi - Z2) -J2 ^2n+l (7.65) 



n=0 



The product of two massless form-factors for the left-moving and the right-moving field above 
is in fact equal to \F]sr{6i, 6n)\'^, the modulo square of the bulk massive form-factor (with 
the only difference that 6 is now the massless rapidity). 

The first term in ( [7.65| ) is Go, the Green's function of the free massless fields on the plane: 

The explicit expressions for the first few terms /2n+i are: 

I, = J^ ^e^^ / -(.i+.,-i,i+i,,)^e« ^ \ tll^ (7.66) 

2N0J-OC 277 V J 2 ^ ^ 



Changing variables, p = exp{6), it can be rewritten as 

AttNq Jo ^ + Tb 



''-£k L PrfP^(^~'"'**"" + «). Ts^e" (7,67) 



and 



,^HI rf[f!^t,„hi^ (7.68) 



3!iVo J ^00 f-J^ 271 2 



i<j 



2(1 + cosh(^i - Oj)) 




Plot of these two contributions to the two-point correlation function for the points OA 
(xi =1/1 = 1/2 = 0) is shown in figures 7.5 and 7.6, and few values of are given in Table 
7.1. It turns out that the series J^^n converge fast because each term is by the factor 
of hundred smaller than In-2- We checked it up to n = 5. Each integral converges for 
any finite value of (xi -|- X2), but is divergent for xi = X2 = 0. Therefore, we need to do an 
additional work to find the correlation function between two points on the boundary, OB 
(Fig. 7.4). 

Let us find equivalent expressions for the integrals /„ which would be finite for xi = 
X2 = 0- For this, we note that the contour of integration can be rotated to go from to 
ioo in the complex p-plane, eq. ( [7.67] ). This rotation is equivalent to the shift of countour 
of integration in ^-plane up by i7r/2, with having the contour to pass below the poles of 
K{6) (figure 7.7). For the complex conjugated term one has to rotate countour clockwise in 
p-plane (shift by —i7i/2 in 6'-plane). So, we obtain 

n?^(^0 Ee'M |F„(^^i,...,^„)re-(^-"-^-HE:U<=«0+c.c (7.69) 
-00 [.^-^ Ztt J \.^-^ ) 

where we used eq. ( |7.54| ) and yi — j/2 > is assumed. 
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0.02 - 



0.01 




< y C "0 "0 "OOP ^ 



-0.01 



-0.02 



-0.03 I 1 1 1 1 1 

2 4 6 8 10 

Figure 7.5: Plot of the one-particle contribution to correlator between two points (0, 0) and 
{x, 0) as a function of x. 



0.0002 



0.0001 




-0.0002 I 1 1 1 1 1- 

2 4 6 8 10 



Figure 7.6: Plot of the three-particle contribution to correlator for two different couplings; 
the curve for B=0.1 is shown scaled by the overall factor 10 compared to the true curve. 
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X 


h 


1.1 


-4.65436256 ■ 10'^ 


2.1 


-2.02534138 ■ 10'^ 


3.1 


-1.11571005 ■ 10-6 


4.1 


-7.03902032 ■ 10"^ 


5.1 


-4.81526684- 10"^ 



Table 7.1: Five-particle contribution for i? = 1. 




Figure 7.7: The contour of integration in the complex 9 plane. The vertical intervals of the 
contour are assumed to be sent to infinity. 



Notice that this expression could be also obtained if we had started from the dual picture 
where time goes along the boundary and space is a half-line and inserted a full basis of 
asymptotic states of the form |^i, 6'„)i -|- R{6i)...R{0n)\0m ■■■^Oi) r between the fields in 
the correlator. This dual picture is to be discussed in detail in the next section. 

For the points OB in Fig. 7.4 (xi = X2 = Q,y2 = 0) expression ( |7.69| ) reads: 



°° d6i...d6n 



D (27r)" 
E.g. for n = 1 we have 



.1=1 



i=l 



, . . . , 9j^ 



|2g-?/iM 



(EL,' 



vt=l 



(7.70) 



fx 



2'KNr 



/oo „ 
-oo 



7.4.2 Boundary-in-space representation 

One can in principle compute the same correlation correlation functions adopting a different 
point of view on the space-time geometry (figure 7.4). In euclidean field theory the role of 
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space and time can be interchanged. Now, in the alternative representation, x > is a space 
coordinate, and y is an imaginary time coordinate. The theory is invariant under the global 
time shift, which is performed by the operator T = exp{yH). 

However, the theory, being defined on a semi-infinite line, has rather intricate structure 
of the ground state, and the form-factors are not known to us a priori. In particular, the 
ground state is not invariant under the space translations. One can speculate that the 
ground state for a theory on a half-line ||0) can be obtained from the bulk vacuum |0) by 
the action of a "boundary creation operator" B: ||0) = -B|0). Speaking loosely, boundary is 
an infinitely heavy inpenetrable particle sitting at the origin. The asymptotic particle states 
are not pure L or R moving, but rather the superpositions of those. E.g., one-particle states 
are 

\\9) = \9)l + R{9)\9)r, 9>0 (7.71) 

To find the correlation function {0\\dx4>{xi, yi)dx(p{x2, y^) 1 10) we conjecture that the following 
equality holds: 

{0\\dx(j){xi,yi)dx(f){x2,y2)m = }^ --— , (7.72) 

with the particular ansatz for the intermediate states \ \n): 

\\9u ...,9n) = \9i,..., 9:^)l + R{9^) ■ ■ ■ R{9n)\9n, . . . , 9i)n (7.73) 

(in the bra states, (?t,||, the factors R{9) become complex conjugated). Substituting ( [7. 731) 
into ( |7.72| ), combining the complex conjugated terms in the series and using the bulk form- 
factors, we indeed get the correct expression for the correlator ( [7.69|) . In the course of this 
calculation one has to use the unitarity of the reflection matrix, R{9)R{9) = 1, as well as 
the property that the vacuum average of the quantum field vanishes, (0) = 0. 

Note that in principle it is possible to formulate the set of general equations for the 
form-factors of integrable models on a half-line, (O||0(O)|n), in analogy with the equations 



of [^]. E.g., it is obvious that the following should hold: 

mm\o) = R{0m\m\-o)- 



7.4.3 The renormalization group analysis 

Let us study the behavior of integrals /„ under the change of scale z e^z. Such a rescaling 
can be compensated by the change 9b ^ 9b + X and by the overall normalization factor 
Z{X) = e^'^ to have the integrals (and hence the correlator) unchanged. Repeating this RG 
transformation, we will flow to the UV or IR fixed points 9b = ±C)0 (depending whether 
A > or A < 0). For such values of 9b the hyperbolic tangent factors in the integrands are 
equal to ±1, and the integrals are proportional to ±l/z^. On the plots Figure 7.5 and Figure 
7.6 one can see two regimes: /xz <^ 1 and fiz ^ 1 when the functions behave as ±l/z'^. The 
non-trivial behaviour at the intermediate scales is due to the presence of boundary, which 
introduces a scale /ze^^ corresponding approximately to the position of the deep. Shifting 
9b corresponds to the motion of the deep to the right or left on Figures 7.5 and 7.6, untill it 
will go away completely and one of the regimes will dominate over all scales. 
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7.4.4 The use of Kubo's formula 



To use the Kubo formula ( |7.3| ), we adopt the first point of view where the boundary is taken 
into account through the introduction of the boundary state \B). Write : 



{d,(j){x,y)d,,(j){x',0)) 



dEg{E) exp [-E{x + x') - iEy] . 



(7.74) 



One obtains g{E) simply by fixing the energy to a particular value in ( [7.65| ). When this is 
done, the remaining integrations occur on a finite domain for each of the individual particle 
energies since J^fZt'^ = E, and there is no problem of convergence anymore. One then 
gets : 



QiE) 



oo rlnE 

E 



n=0 



2n 



E^ 



3 (27r)2"+i(2n + l)!E-e^i 
K{e^)...K{e2n)K[\n [E-e'^-. 

F2n+1 {ei...e2nM [E - c'^ - . . . - 



— e 



e2n 



(7.75) 



with the constraint e^' < E. The denominator might suggest some possible divergences; 
it is important however to realize that it vanishes if and only if the particle with rapidty 
(^2n+i has vanishing energy, in which case the form factor vanishes too. By using the dual 
picture, one writes : 

(9,-0(x, y)d,'(j)ix', y') >= / dET{E) exp [-tE{x + x') - E{y - y')] , (7.76) 

Jo 

The two expressions are in correspondence by the simple analytic continuation : 

g{E) = iJ^{iE). (7.77) 



Eq. ( [7.74|) is the only correlation contributing to AG for positive Matsubara frequencies, 



and 



4 um 



(7.78) 



Here we have used the fact that ojmL <^ 1, i.e. the system, although large, is much smaller 
than the wavelength associated with the (modulus of the) AC frequency. To go to real 
frequencies, we can simply substitute uj^i —>■ ~iuj in the K matrices in the integrals ( 7.65| ): 



^Giuj) = -^Imgi-iu) = -^ReJ^(cj). 

4:U Aui 



(7.79) 



7.4.5 The numerical work 



We computed the integrals for 3 and 5 particle contributions to correlation functions numer- 
ically using Monte-Carlo simulations. The domain of integration was the hypercube with 
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the length of the side equal 40. Because of exponential decay of the integrand, the con- 
tribution of the region outside the box is small withing our accuracy. Since the integrals 
diverge as the imaginary time approaches zero, we found it technically more difficult to find 
the reliable results for a; < 0.1. We took 5 ■ 10^ points to evaluate the integrand inside the 
box. We checked that the Monte-Carlo result is stable by increasing the number of points 
and estimated the relative error in 3 and 5 dimensional integrals to be about 1%, while the 
1-dimensional integrals were evaluated with 0.001% accuracy. The higher-particle contribu- 
tions can be computed in the same manner, but require considerable amount of computer 
time. Because we expect them to be very small, we found their evaluation unnecessary for 
the purpose of the present research. 



7.5 Boundary sine-Gordon correlators 

We wish to apply now the formalism developed in the previous section to the massless 
boundary sine-Gordon model, with B = M cos 0(0). The basic procedures and formulas of 
the previous section apply directly to the present section, the difference appearing in a few 
technicalities. 



7.5.1 Conductance at g = 1/3 

We are interested in the value of coupling constant = Svr/S in ( [7.49| ), which corresponds 
to the quantum Hall liquid regime u = 1/3. The particle spectrum in the sine-Gordon model 
depends on the coupling constant. At g = 1/3, we have three particles in the spectrum: a 
soliton, antisoliton and their bound state - breather. The breather is completely analogous 
to the scalar particle of the sinh-Gordon model, while the other two particles introduce 
technical differences. 

Let us describe what these differences are. In the boundary state, (|7.63|) , one has to intro- 
duce soliton and antisoliton creation operators A^~^^ and A^~\ and their boundary reflection 
matrices: 



\B) = \0)+Y^ I ^ ^ K^'^'{9i) . . . K^^^^ 

N=l 



7iV 



-co<Si<...<9]v<oo 

AT (6^) . ..Al^{e,)A'AiOi) . ..A'n^e^m (7.80) 
(the summation over the particle indices a, b is assumed) 

K^'^e) = Ri(— + e 



The breather reflection coefficient is still given by formula ( [r.53| ), while for the soliton and 
antisoliton massless reflection matrices are 



R^{e) = e-v^i?(^), (7.81) 
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Figure 7.8: Frequency dependent conductance at T = 0. 



R{9) 



1 



2n- 



r(l + 



2tt' 



2 cosh(^ - 



ITT ^ 
4 > 



(7.82) 



r(|-f)r(| + f) 

(for the arbitrary (3 expressions we refer to |^). Correspondingly, there are more terms 
in the series expansion of the correlators since there are more intermediate particles. The 
leading contribution, Ji, is still the one-particle contribution of the breather. The next- 
after-leading contribution comes from two particles in the intermediate state, soliton and 
antisoliton. Its magnitude is approximately 20% of the value of Ji. Note that soliton and 
antisoliton can appear only in pairs since the total topological charge of the intermediate 
states should vanish in order for it to contribute to the correlator. The next terms come 
from the three-particle states of the breather with soliton and antisoliton, and the three 
breathers. Their magnitude is 1% of the value of the leading contribution. The form-factors 
for different particles in the sine-Gordon model have been originally obtained by Smirnov in 
the massive case |j72|, and their massless limit is given in |70|. These expressions are much 
more cumbersome than that of sinh-Gordon model. 

The leading contribution to the conductance computed along the lines of the discussion 
above is given by 



AG{uj) 



(1) 



2 1 

— Retanh 
8 



1 / UJ \ 171 



(7.83) 



where /i ~ 3.14 and d ~ 0.1414. The contribution from the soliton-antisoliton state can be 



found in We plot the full function G{uj) in figure 7.8. 
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7.5.2 The free point, g = l/2 

In the free case (3'^ = Att one has simply: 



R^{9) = p{e) 
Rt{0) = Q{e) 



+ i 
i 



& + 1 

As for the form-factors, only the soliton-antisoliton form-factor is non-zero, 
Thus, we find for the conductance 



(7.84) 



1 — tan ^{ijj /Ti 



B) 



(7.85) 



This is in agreement with the solution of ||7^ and also 
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Chapter 8 
Conclusion 



This review is concerned with the study of 2D integrable models with boundary. In chapters 1 
through 6 the formahsm is developed for dealing with such models and particular examples 
are discussed. It involves a construction of solutions to the classical equations of motion 
for a model on a half-line, solution of the quantum XXZ chain in a boundary magnetic 
field, investigation of the boundary bound states (a phenomenon caused by the presence of 
boundaries), generalization of the Destri-deVega technique, which allows to find the exact 
ground state scaling energy, to the theories with boundaries. All this work was carried 
out mostly for the quantum sine-Gordon/Thirring model, which provides a good basis for 
theoretical investigations withing the framework of quantum integrable models. 

There is a room for the further development of the above formalism. For example, we 
often have chosen the boundary conditions to be a specific, Dirichlet boundary condition. 
Consideration of more general boundary conditions will most certainly be associated with 
the increased complexity, while it is not clear to us whether the outcome will prize us with 
a new interesting physics, or at least whether the most general boundary condition has any 
vital applications. An interesting direction of research is associated with the technique of 
the thermodynamic Bethe ansatz and the corresponding finite-difference equations p7| . 
Another interesting direction has to deal with the non-integrable deformations of integrable 
models [||. 

In chapter 7 we described some of the applications of boundary integrable models to 
condensed matter physics. In particular, we developed a technique to compute exactly 
time dependent properties of the Kondo problem and/or two-state problem of dissipative 
quantum mechanics, as well as the conductance of the one-dimesional wires with a tunneling 
through impurity. These problems have potential experimental applications and a room for 
further theoretical work. The latter includes the computation of the voltage and temperature 
dependent behavior of the conductance, quantum noise, as well as studying of the Coulomb 
blocade phenomenon. 
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